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ABSTRACT 
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An  echosounder  is  used  to  probe  various  atmospheric 
parameters.  An  acoustic  wave  is  transmitted  into  the 
atmosphere  and  information  deduced  from  the  backscattered 
energy . 

This  thesis  seeks  to  understand  the  range  limitations  of 
the  echosounder  and  to  explore  methods  to  quantify 
atmospheric  turbulence  parameters  at  a  given  range.  The 
propagation  of  the  acoustic  energy,  including  the  effects  of 
excess  attenuation,  are  modeled  to  predict  the  performance 
of  an  echosounder  when  various  parameters  are  changed.  The 


electronics  of  an  existing  echosounder  are  investigated  to 

ts 

understand  inherent  or  design  limitations.  ■ 
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I.  INTRODUCTION 


Many  atmospheric  turbulence-dependent  processes  take 
place  in  the  troposphere,  the  lowest  ten  to  twenty 
kilometers  above  the  earth's  surface.  Acoustic  energy 
interacts  with  irregularities  in  the  atmosphere  more 
strongly  than  electromagnetic  waves  and  is  potentially  a 
better  probe  for  these  irregularities.  Acoustic  echoes  have 
been  used  to  detect: 

1)  wind  speed  and  direction  profiles, 

2)  humidity  profiles, 

3)  strength  and  location  of  temperature  inversions, 

4)  temperature  inhomogeneities, 

5)  mechanical  turbulence. 

The  acoustic  sounder,  also  known  by  the  names 
echosounder,  echosonde,  sodar,  and  acoustic  radar,  transmits 
a  pulse  of  acoustic  energy  into  the  atmosphere.  The  various 
atmospheric  parameters  can  be  determined  based  on  the 
intensity  and  frequency  of  the  scattered  energy.  The  range 
to  the  scattering  center  is  determined  from  the  elapsed  time 
between  transmission  and  return  of  the  scattered  energy. 

In  this  thesis  we  are  attempting  to  improve  our 
understanding  of  the  fundamental  range  limitations  of  the 
acoustic  sounder.  In  general,  the  acoustic  sounder  suffers 
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from  multiple  scattering.  Turbulence  at  shorter  ranges 
alters  the  phase  front  of  the  propagating  pulse  and  reduces 
the  magnitude  of  the  return  from  longer  ranges. 

We  are  also  exploring  the  use  of  return  from  shorter 
ranges,  along  with  theory  for  the  degradation  of  the  energy 
along  the  path,  to  compensate  for  short  range  degradation  in 
a  boot-strap  fashion  and  quantify  the  return  from  longer 
ranges.  Acoustic  sounders  at  present  indicate  the  presence 
of  inhomogeneities  at  a  specific  altitude  but  the  magnitude 
of  the  inhomogeneities  may  have  errors  of  a  factor  of  four 
or  more. 

The  following  factors  affect  the  range  of  the 
echosounder : 

1)  Atmospheric: 

a)  Pressure, 

b)  Temperature, 

n 

c)  Temperature  Structure  Parameter  CT  , 

2 

d)  Velocity  Structure  Parameter  Cv  , 

e)  Water-vapor  Pressure,  and 

f)  Ambient  Noise. 

2)  Echosounder: 

a)  Antenna  Aperture  Factor, 

b)  Antenna  Diameter, 

c)  Efficiency, 

d)  Frequency, 


II.  BACKGROUND 


Acoustic  energy  propagates  in  the  atmosphere  as  a 
longitudinal  wave  of  pressure  variations.  The  energy  is 
scattered ,  attenuated ,  and  refracted.  Energy  which  is 
scattered  constitutes  the  returned  signal.  The  attenuation 
decreases  the  energy  ensonifying  a  given  volume  and 
decreases  the  returned  signal.  Energy  that  is  scattered 
through  small  angles  or  that  is  refracted  results  in  further 
degradation  or  excess  attenuation  of  the  signal. 

A .  SCATTERING 

The  echosounder  transmits  a  pulse  of  acoustic  energy 
that  is  scattered  by  temperature  and  velocity 
inhomogeneities.  Information  about  the  inhomogeneities  is 
then  based  on  elapsed  time  between  transmission  and  receipt 
of  the  return  signal,  the  strength  of  the  returned  signal, 
the  equations  for  scattering,  and  the  Doppler  shift. 
Variations  in  the  propagation  velocity  of  the  wave,  which 
are  a  function  of  the  temperature  and  velocity  variations  in 
the  air  mass,  produce  the  scattering,  in  addition  to 
scattering,  the  temperature  variations  cause  refraction 
while  the  velocity  turbulence  causes  both  a  shift  in  phase 
and  direction  of  propagation  [Ref.  l:p.  60]. 
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turbulence  and  uniform  gradients  of  temperature  or  velocity 
contribute  to  the  scattering.  The  gradient  must  change  in 
a  scale  size  that  is  comparable  to  or  smaller  than  the 
acoustic  wavelength  in  order  to  contribute  to  the 
scattering.  It  appears  that  the  effects  of  uniform 
gradients  are  limited  to  beam  bending  and  possible  specular 
reflection  for  the  acoustic  frequencies  of  interest  here. 
This  refractive  structure  of  the  troposphere  and 
stratosphere  would  be  better  probed  with  what  has  come  to  be 
called  infrasound,  sound  of  frequencies  below  twenty  Hertz. 
(Ref.  1 : p .  52] 

The  turbulent  atmospheric  temperature  and  velocity 

fluctuations,  at  high  wave  number,  follow  the  Kolmogorov 
K  5//3  power  spectral  density  law  (in  one  dimension).  The 
fluctuations  are  correlated  spatially  on  the  order  of  one 
centimeter  to  tens  of  meters.  The  nonzero  correlation 
lengths  and  the  declining  power  spectral  density  imply  a 
different  scattering  process  than  for  random  point 
scatterers,  even  though  the  turbulence  is  a  stochastic 
process  (Ref.  l:p  61]. 

An  expression  for  the  power  scattered  from  a  unit  volume 
per  unit  incident  flux  into  a  unit  solid  angle  is  [Ref.  2:p 
84] 

Of  =  s-Q-55  cos2©  (£y2  cos2(2)  +  .13  Ct2 ) ( sin(f ) )”11/3  , 
y'*  cz  2  Tz  1 

I  o 


where 


is  the  wavelength  of  the  transmitted  sound, 

T  is  the  average  atmospheric  temperature, 

0  is  the  angle  of  scattering  from  the  direction  of 

propagation, 

2 

Cv  is  the  velocity  structure  parameter, 

Cv2“  <VOO  j/yxtr)>2  , 

and  Ctz  is  the  temperature  structure  parameter, 

Ct2=  <T(x)  -^(x+r)>2  . 

This  scattering  equation  indicates  that  backscatter 
(O -  rr )  is  only  a  function  of  the  temperature  structure.  The 
velocity  at  any  point  is  the  vector  sum  of  the  phase 
velocity  of  the  sound  and  the  particle  velocity  of  the 
turbulence.  This  vector  sum  is  always  in  the  forward 
direction  because  the  particle  velocity  is  always  less  than 
the  speed  of  sound  c.  Therefore  the  backscatter:  is  only  a 
function  of  the  temperature  structure  [Ref.  l:p  60],  The 
vertical  turbulent  velocity  does  cause  a  Doppler  shift  of 
the  frequency  of  the  backscattered  radiation. 

Scattering  over  a  region  with  correlated  scattering 
centers  produces  constructive  and  destructive  interference. 
The  backscattered  waves  are  partially  coherent.  This 
results  in  much  greater  intensities  than  would  be  received 
from  incoherent  scattering  such  as  Rayleigh  scattering  [Ref. 


l:p  61].  Also  the  scattering  with  the  interference,  over 
the  nonzero  coherence  length,  acts  like  an  array  of 
scattering  centers  in  a  regular  crystal  lattice  of  spacing 
L.  The  Bragg  condition 
L»  ^  /sin(‘ft/2) , 

is  satisfied.  The  dominant  scattering  is  for  scale  sizes  of 
X/2.  [Ref.  lsp  61] 

B .  ATTENUATION 

The  atmosphere  absorbs  some  acoustic  energy  that 
propagates  through  the  atmosphere  and  reradiates  this  energy 
at  different  frequencies.  Quantitatively  the  power  lost 
(P^)  over  a  path  length  1  is  given  by 


where  is  the  attenuation  coefficient  in  Nepers  per  meter. 
Historically  the  attenuation  coefficient  has  been  considered 
to  be  the  sum  of  several  terms; 
c*;  cl=  classical  viscous  losses, 
r=  molecular  rotation  losses, 
cKvib*  molecular  vibrational  losses,  N2  and  o2* 

Classical  and  rotational  loses  are  negligible  below 
about  three  kilohertz,  the  region  of  interest  for 
echosounder  operation  [Ref.  l:p  54  and  Ref.  3:p.  18-2], 
consequently  vibration  of  N2  and  02  produce  most  of  the 


attenuation. 


Rotational  and  vibrational  losses  are  referred  to  as 
relaxation  processes.  The  acoustic  energy  excites  internal 
energy  modes  of  the  Nj  and  02  molecules.  The  rate  of 
collisions  with  water  vapor  determines  the  rate  and 
efficiency  of  conversion  of  the  energy  into  translational 
energy  (heat). 

The  phase  is  shifted  due  to  the  relaxation  processes. 
This  is  one  of  the  reasons  some  of  the  microwave  radar  pulse 
compression  techniques  cannot  be  used  in  echosounders .  A 
number  of  the  pulse  compression  techniques  rely  on  the  phase 
not  changing  during  propagation. 

The  dependence  of  the  attenuation  on  the  water  vapor 
pressure  is  believed  to  be  due  to  a  resonance  process 

between  the  lowest  vibrational  states  of  the  02  and  N2 
molecules  with  the  water  molecules.  Henderson  and  Hertfeld 
in  Reference  4  [p.986]  state  that  the  lowest  vibrational 
states  of  O2  and  water  vapor  are  only  thirty-nine 
wavenumbers  (56°K)  apart  at  1600  cm”*.  For  this  reason 
oxygen  was  thought  to  be  primarily  responsible  for  the 
humidity  dependent  absorption  of  sound.  Henderson  and 
Herzfeld  (Ref.  4],  and  many  others,  assumed  nitrogen  was  an 
inactive  dilutent  having  no  effect  on  the  relaxation 
processes  (Ref.  4:p.  987].  Unfortunately  the  theory  did  not 
agree  with  the  data  for  low  frequencies  and  high 
humidities . 


Theory  incorporating  the  relaxation  processes  involving 
nitrogen  and  water  at  low  frequencies  and  relative 
humidities  greater  than  twenty-five  percent  bring  the  theory 
into  agreement  with  the  data.  At  high  frequencies  and  low 
relative  humidities,  theory  and  data  match  well  with  oxygen 
making  the  main  contribution  and  nitrogen  acting  as  an  inert 
dilutent.  At  low  frequencies  and  high  humidities,  nitrogen 
seems  to  make  the  main  contribution  and  oxygen  acts  as  an 
inert  dilutent.  For  the  range  of  relative  humidities  in 
which  many  people  live  and  over  a  good  part  of  the  audible 
frequency  range  it  appears  nitrogen  is  the  major  contributer 
to  the  relaxation  processes  [Ref.  5:p.  165].  This  is  in 
contrast  to  the  previous  assumption  that  oxygen  was 
responsible  [Ref.  6:p.  604]. 

An  expression  for  attenuation  due  to  the  molecular 
absorption  is  given  by  [Ref.  7sp.  34], 

=  ^ax  £oi8»f  \2  +/2»(f/fm)2J\  2  1  1/2  , 

304. 8^  fm)  Kl  +<f/fm)  /  J 

in  dB/m, 

where 

f  is  the  frequency  of  the  transmitted  pulse. 


f  is  the  Napier  frequency. 


fm  =(10+6600  h+4400  h2) 


(  *\  (  .aia,  'l 0J 

\  10T?  )  \1 . 8T+492  / 


*  =.0078  fj  1.8  T+492V2 * 5  exp(7 .77(1-  519  A\ 

>  \  V  1.8T+492/; 


where 

T  is  temperature  in  degrees  Celcius, 
h  is  the  percent  mole  ratio  of  water  vapor, 
h=  100  (e/p) , 

e=  water  vapor  pressure  in  mb, 

P=  atmospheric  pressure  in  mb. 

To  convert  from  dB/m  to  Nepers/m  note  that 
10  log(I(x)/I(0)  )=  10  log (exp( -Xx ) , 

==>  4.34  X (Nepers/m )=  X(dB/m). 

The  Napier  frequency  is  the  frequency  of  maximum 

absorption  per  wavelength  and  “^max  is  the  attenuation  at 
the  Napier  frequency.  The  Napier  frequency  is  shifted  to 
higher  frequencies  by  even  small  amounts  of  water  vapor. 
Plots  of  attenuation  versus  water-vapor  pressure  in 
millibars  (Figure  1)  and  versus  relative  humidity  (Figure  2) 
are  shown  for  a  range  of  frequencies. 

For  a  given  temperature,  the  attenuation  has  a  maximum 
and  decreases  for  higher  or  lower  humidities.  With  higher 
temperatures,  the  maximum  attenuation  increases  and  the 
relative  humidity  at  which  that  maximum  occurs  decreases. 
Plots  of  attenuation  versus  water-vapor  pressure  in 
millibars  (Figure  3)  and  versus  relative  humidity  (Figure  4) 
are  shown  for  a  range  of  temperatures. 
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Pig.  3.  Attenuation  versus  Water-vapor  Pressure  (mb.) 
for  Temperatures  Around  20  Celsius. 


Fig.  4.  Attenuation  versus  Relative  Humidity  (!) 
for  Temperatures  around  20  Celsius. 


C .  EXCESS  ATTENUATION 

The  atmosphere  scatters  and  attenuates  the  acoustic 
energy.  In  the  previous  section  on  scattering,  we  were  only 
concerned  with  the  energy  scattered  from  a  given  volume. 
Excess  attenuation  accounts  for  scattering  and  energy  loss 
for  the  round  trip  in  the  atmosphere  to  and  from  a  given 
volume  of  interest.  Excess  attenuation  arises  because  the 
atmosphere  degrades  the  mutual  coherence  of  the  acoustic 
wave.  The  divergent  solid  angle  of  the  acoustic  wave  is 
larger  than  would  occur  for  a  coherent,  diffraction  limited 
wave. 

The  excess  attenuation,  Ze2,  accounts  for  this  energy 
lost  due  to  small-angle  scattering.  This  turbulent  beam 
broadening  reduces  the  on  axis  intensity.  Clifford  and 
Brown  in  Reference  8  [p.  1972]  develop  the  equation 
Ze  =  1/ ( 1+N )  for  N< 1 , 

Ze  =1.5/ (1+N)  for  N> 1 , 
where 


Dq=  antenna  diameter, 

y^Oe  is  the  atmospheric  acoustic  coherence  length 

=  1.46*k2  *Qf  R°  ds*Cne2(s)*(  (k^)5/3  +  [  s_]5/3)"3/5  . 

2R0  2R(j 

The  term  varies  from  a  value  of  one  for  no  excess 


attenuation  to  an  asymptote  of  zero,  implying  the  energy 
would  be  spread  over  a  2  tTsolid  angle. 

There  is  a  step  in  the  functional  dependence  of  the 
excess  attenuation  when  N=1  that  can  be  seen  in  the 
equations  above  and  in  the  plots.  There  clearly  is  no 
physical  discontinuity  but  rather  a  transition  between  the 
theoretical  dependence  between  two  asymptotic  regions. 

Figures  5,6,  and  7  are  plots  of  excess  attenuation 
versus  range  for  antenna  sizes  of  .5,  1,  and  1.5  meters, 
respectively,  for  several  frequencies.  As  can  be  seen  from 
these  plots  or  the  equations  above,  the  relationship  of  the 
antenna  size  to  the  coherence  length  has  a  significant 
effect  on  the  excess  attenuation.  A  larger  antenna  will 
not,  by  itself,  produce  greater  range.  The  antenna  size 
must  be  matched  to  the  transmission  frequency  in  terms  of 
antenna  design  guidelines  and  effects  upon  the  excess 


attenuation. 
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Pig.  6.  Range  versus  Excess  Attenuation 
£or  an  Antenna  Diameter  of  1.0  m. 


Pig.  7.  Range  versus  Excess  Attenuation 
for  an  Antenna  Diameter  of  1.5  m. 


2  I 


Ill  EXPERIMENT 


The  acoustic  sounder  is  an  excellent  tool  for  probing 
the  lower  troposphere.  It  suffers  from  several 
shortcomings : 

1)  The  range  has  typically  been  limited  to  a  few  hundred 
meters . 

2)  It  has  proven  difficult  to  quantify  the  measurements 
accurately  for  a  given  range. 

In  this  thesis,  we  approached  the  problem  from  several 
directions.  We  looked  at  the  hardware  to  understand  the 
signal  to  noise  limitations  and  to  understand  how  the 
magnitude  of  the  returned  signal  could  be  calibrated.  We 
developed  a  software  model  which  would  allow  one  to  estimate 
the  power  backscattered  from  the  atmosphere  at  a  given 
range,  based  on  profiles  of  atmospheric  characteristics  and 
input  parameters  for  the  acoustic  radar. 

A .  HARDWARE 

The  acoustic  sounder  we  were  working  with  was  the 
Aerovironment  Model  300.  It  consists  of  an  electronic 
module  which  generates  a  1600  Hertz  electrical  pulse.  The 
pulse  is  converted  to  acoustic  energy  by  a  transducer  which 
feeds  a  1.25  meter  parabolic  reflector.  Energy 


backscattered  by  the  atmosphere  is  then  received  by  the 
reflector  and  transducer.  The  electrical  signal  is  then 
filtered  and  amplified.  In  addition,  a  ramp  amplifier 
compensates  for  the  1/r  decreasing  signal  amplitude  with 
range  to  decrease  the  dynamic  range  requirements  of  the 
present  data  display,  a  strip  chart  recorder. 

We  replaced  the  various  integrated  circuit  amplifiers 
and  filters  with  more  current  designs  with  lower  noise.  We 
replaced  the  pre-amplifier  with  an  OPA  111  and  the  rest  with 
LF  356  BN  devices. 

We  traced  the  amplifier  and  filters  of  the  receiver 
board  to  understand  the  undocumented  choices  the 
manufacturer  had  made.  Figures  8  and  9  show  a  pre¬ 
amplifier,  high  and  low  pass  filters  and  two  stage 
amplification.  The  filters  are  bi-quad  configured  with 
notched  outputs.  Figure  10  shows  the  notching.  Figure  10 
represents  the  frequency  spectrum  output  of  the  receiver 
board  as  measured  with  the  HP  3561A  Frequency  Spectrum 
Analyser  with  random  noise  from  the  HP  3561A  providing  the 
input  signal  before  the  filters.  Figure  11  represents 
frequency  spectrum  of  the  receiver  board  with  the  random 

noise  across  the  transducer.  A  10^  ohm  resistor  was  used  to 
match  impedences  as  shown  in  Figure  12. 


Fig.  12.  Connection  of  Noise  Source  to  Transducer 
The  manufacturer  had  a  filter-oscillator  board  after  the 
receiver  board  with  a  switch  for  three  pass  bands;  narrow, 
medium,  and  wide.  With  the  same  input  as  in  Figure  10  we 
measured  the  frequency  output  for  the  three  settings  of  the 
filter-oscillator  board  (Figure  13).  As  can  be  seen,  there 
is  little  difference  between  the  three  settings. 

From  previous  measurements,  the  resistor  labeled  R25  in 
the  schematic  (Figure  8)  of  the  receiver  board  was  found  to 
load  down  the  input.  The  resistor  was  initially  25  kilo- 
ohms.  Figures  14,  15,  and  16  represent  the  signal  from  the 
input  transformer  with  R25  =25K,  100K,  and  oO  ohms 
respectively.  Random  noise  from  the  HP  3561A  was  input 
across  the  transducer  as  shown  in  Figure  12.  The  Q  of  the 
input  transformer  was  improved  by  increasing  the  resistance. 
The  lower  curve  in  each  Figure  is  with  no  noise  signal  input 
from  the  HP  3561A. 

It  appeared  the  Q  of  the  filters  could  also  be  increased 
by  adding  the  resistors  labeled  Rq  in  Figure  17.  The  upper 
curves  in  Figures  18,  19,  and  20  are  the  output  of  the 
receiver  board  with  each  Rg  =  c,°,  750,  and  560  ohms 
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Pig.  13.  Frequency  Spectrum  of  the  Pilter  Oscillator  Board 
of  the  Aerovironment  Model  300. 

(Input  before  the  filters  of  the  receiver  board) 
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Pig.  14.  Frequency  Spectrum  of  the  Input  Transformer 
(Input  across  transducer  with  R25*25k  j 
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Pig.  18.  Frequency  Spectrum  of  the  Receiver  Board  with  Rq*  <>^ 

(input  before  filters) 


Fig.  19.  Frequency  Spectrum  of  the  Receiver  Board  with 
rq-750  -a  (input  before  filters) 


respectively.  The  input  was  random  noise  from  the  HP  3561A 
after  the  pre-amplifier  but  before  the  filters  of  the 
receiver  board.  The  effect  of  decreasing  the  resistance 
was  an  increase  in  gain  but  the  pass-band  was  no  longer  flat 
on  top.  However,  for  the  range  of  vertical  wind  velocities, 
about  two  or  three  meters  per  second,  a  pass-band  of  fifty 
or  sixty  Hertz  is  adequate  to  allow  for  the  Doppler  shift. 

As  can  be  seen,  the  filters  could  be  tuned  for  a  narrow 
pass-band  which  would  have  better  shape.  We  would  gain 
amplification  and  still  have  a  pass-band  that  is  wide  enough 
to  meet  our  needs. 

The  lower  curves  in  Figures  18,  19,  and  20  are  with 
notching,  the  upper  without.  The  effect  of  the  notching  was 
a  sharper  cut  off  near  the  bandpass  frequency  (1600  Hz)  but 
a  loss  at  higher  and  lower  frequencies. 

B .  SOFTWARE 

The  signal  received  by  an  acoustic  sounder  from  a  given 
range  tells  us  that  scattering  centers  exist  but  we  can  say 
little  about  the  size  and  magnitude  of  the  scattering 
centers.  We  don't  know  how  much  energy  was  incident  on  the 
scattering  volume  nor  how  much  the  signal  was  degraded  on 
the  return  path.  If,  on  the  other  hand,  previous  signals 
are  used  to  determine  the  atmospheric  characteristics  then 


an  estimate  of  the  degradation  of  the  incident  power  and  the 
degradation  of  the  return  signal  can  be  made. 

The  software  programs  I  have  written  here  go  in  the 
opposite  direction.  Given  certain  atmospheric 
characteristics,  the  power  returned  is  estimated.  With  the 
modules  tested,  it  would  then  be  a  matter  of  turning  it 
around  to  take  actual  data  and  estimate  atmospheric 
conditions.  In  the  present  form  they  serve  our  purposes  in 
allowing  us  to  explore  the  effects  of  parameter  changes. 

All  the  programs  are  fundamentally  built  around  the 

echosonde  equation,  also  sometimes  referred  to  as  the  radar 

equation  in  meteorology  [Ref.  9  and  Ref.  7:p.  3]. 

Pr  =  Er*(Pt*Et)*(exp(-2XR)  )*  5'0(Rrf)*(G*r  )*(A^)*Ze2, 

2  R 

where 

Pr  is  the  power  returned  from  a  range  R, 

Pt  is  the  power  transmitted  at  frequency  f, 

Er  is  the  efficiency  of  conversion  of  acoustic  power 
to  electrical  power  by  the  transducer  , 

Et  is  the  efficiency  of  conversion  of  electrical  power 
to  acoustic  by  the  transducer, 
c  is  the  speed  of  sound  in  m/sec, 

T  is  the  pulse  length, 

A  is  the  area  of  the  antenna, 

R  is  the  range,  and 

G  is  the  effective-aperture  factor  of  the  antenna. 


Transducer  efficiencies  must  be  measured  for  each  driver 
and  typically  range  from  a  few  percent  [Ref.  7:p.  3]  to 
twenty  five  percent  iRef.  10:p.  11-10] . 

The  scattering  cross  section  per  unit  volume  *-s 

the  fraction  of  incident  power  backscattered  per  unit 
distance  into  a  unit  solid  angle  at  a  given  frequency.  From 
Reference  7  [p.  4]  and  Reference  2, 

<f0  =  .0039  *  k1/3  *  Qt2  , 

U  IT»  2 

T0 

where 

k  is  the  wavenumber  =  2  ^/wavelength, 

T0  is  the  local  mean  temperature  in  Kelvin, 

Ct*  is  the  temperature  structure  parameter. 

The  power  scattered  from  a  scattering  volume  is 
Pb(  I  )  =  (Pt*Et-Pb(  1-1 )  )*exp(-2»<R)*(c^/2)*Ze2*  CC.  , 

The  power  returned  to  the  antenna  is 
Pr=Pb*A*G*Er/R2, 

where  the  return  path  attenuation  was  already  included  in 
Ps.  The  excess  attenuation  Ze2  was  discussed  in  the 
background  section. 

The  dependence  with  height  of  the  velocity  structure 
parameter  Cv2  was  needed  for  the  calculation  of  the  excess 
attenuation.  Reference  2  and  Reference  11  [p.149]  give 
Cv2  ~  2  **=  ^3 , 

where  £  is  the  dissipation  rate  of  turbulent  kinetic  energy. 


Reference  11  [p.  154]  and  Reference  (p  194)  give 

<£  =  ( . 33m/s ) 3  *  (1  +  .07  *  R3/5)3/2, 

k  *  R 

for  a  stable  surface  layer.  Figure  21  is  a  plot  of  range 
versus  this  Cv  . 

The  programs  have  four  temperature  structure  parameter 
(CT2)  profiles.  The  operator  must  choose  one.  The  first  is 
based  on  data  as  presented  in  Reference  13  [p.  398]  for 
midday  clear  weather  above  the  Tularosa  Basin  desert  in  New 
Mexico.  The  second  choice  is  for  the  same  data  multiplied 
by  a  factor  two  to  approximate  looking  up  a  convective 
plume.  The  third  case  for  a  nocturnal  atmosphere  assumes  a 
dependence  proportional  to  the  negative  exponential  of  the 
range  as  presented  in  Reference  13  [p.  399]  with  tower  data 
from  the  same  reference  used  for  the  first  sixty-five 
meters.  The  fourth  case  assumes  a  dependence  proportional 
to  height  to  the  -4/3  and  a  surface  vertical  heat  flux  of 
.095  [Ref.  7:p.  7],  Also  case  four  allows  for  the  operator 
to  input  the  height  of  an  inversion  layer  with  CT  being 
proportional  to  height  to  the  -4/3  above  the  inversion 
layer.  Figures  22,  23,  24,  and  25  are  plots  of  these  four 
profiles. 

^  2  o 

Cv  and  Ct'5  were  then  used  to  calculate  the  acoustic 
refractive  index  factor  [Ref.  14:p.  119] 


TEMPERATURE  STRUCTURE  PARAMETER 


Fig.  22.  Range  versus  Temperature  Structure  Parameter  (C^,2) 

for  Profile  1 


Fig.  23.  Range  versus  Temperature  Structure  Parameter  (C2) 

for  Profile  2 
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Fig.  24.  Range  versus  Temperature  Structure  Parameter  (CT2) 

for  Profile  3 


TEMPERATURE  STRUCTURE  PARAMETER 

Fig.  25.  Range  versus  Temperature  Structure  Parameter  (CT2) 

for  Profile  4 


cr»e~  =  Ct_  +  Cvf 

4T02  ^ 

Figure  26  is  a  plot  of  Cne2  using  the  third  profile  for  Ct2. 

Figure  27  represents  an  input  flow  chart.  It  summarizes 
the  effect  of  each  variable  or  atmospheric  parameter. 

Figures  28  and  29  are  plots  of  the  range  versus  the  power 
returned  to  the  echosounder  for  profiles  3  and  4. 

The  segments  of  the  programs  were  tested  against 
existing  data  to  verify  proper  operation.  Based  on  input 
temperature,  atmospheric  pressure,  water-vapor  pressure,  and 
the  frequency  of  operation  of  the  acoustic  sounder  the 
programs  calculate  the  attenuation.  If  the  operator 
desires,  the  first  program  will  plot  the  attenuation  as  a 
function  of  absolute  water-vapor  pressure  and/or  relative 
humidity  for  frequencies  at  one-third  octaves  around  the 
input  frequency  and  then  again  for  temperatures  at  ten 
degree  Celcius  intervals  around  the  input  temperature. 

These  plots  were  used  to  check  the  attenuation  against  data 
[Ref.  15]  and  [Ref.  5]. 

Ambient  noise  levels  of  acoustic  sounders  are  found  to 
be  about  ten  to  forty  dB  above  the  theoretical  Johnson  noise 
limit  [Ref.  16:p.  19-4].  This  noise  level  determines  the 
maximum  range.  This  maximum  range  is  comparable  to  range 
capabilities  of  the  Aerovironment  System  300  when  the 
operating  parameters  of  the  Aerovironment  are  used  in  the 
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Fig.  27.  Input  Flow  Chart  for  Computer  Model 
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Trends,  as  input  parameters  where  changed,  were  also 
used  to  check  the  program  output.  As  the  frequency 
increased,  the  backscatter  of  the  turbulence  in  the 
atmosphere  decreased  slowly  and  the  absorption  of  the 
atmosphere  increased  rapidly  in  agreement  with  Reference  17 

[p.  1-10], 

The  calculation  for  excess  attenuation  was  difficult  to 
check.  At  ranges  of  about  400  to  500  meters  the  excess 
attenuation  of  a  typical  acoustic  radar  has  been  found  to  be 
about  .25,  in  agreement  with  the  program  calculation. 

PROGRAM  1 

The  first  program  takes  atmospheric  parameters  and  the 
echosounder  parameters  as  input  and  outputs  nine  plots.  The 
inputs,  by  the  operator  at  the  keyboard,  are: 

Atmospheric  parameters; 

Atmospheric  pressure  in  millibars. 

The  profile  of  Ct2  from  four  options, 

Temperature  in  degrees  Celsius, 

Water-vapor  pressure  in  millibars, 

Echosounder  data; 

Antenna  diameter  in  meters, 

Frequency  of  the  echosounder  in  Hertz, 

Power  transmitted  by  the  echosounder  in  Watts,  and 
Pulse  length  of  the  transmitted  acoustic  energy. 
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The  program  outputs  the  following  plots: 

1.  Attenuation  (1/m)  versus  water-vapor  pressure  (mb)  for 
five  frequencies  at  one-third  octaves  around  the  input 
frequency. 

2.  Attenuation  (1/m)  versus  relative  humidity  (%)  for  five 
frequencies  at  one-third  octaves  around  the  input 
frequency. 

3.  Attenuation  (1/m)  versus  water-vapor  pressure  (mb)  for 
five  temperatures  at  ten  degree  intervals  around  the 
input  temperature. 

4.  Attenuation  (1/m)  versus  relative  humidity  (%)  for  five 
temperatures  at  ten  degree  intervals  around  the  input 
temperature. 

5.  Range  (m)  versus  excess  attenuation. 

6.  Range  (m)  versus  the  temperature  structure  factor  Ct  . 

7.  Range  (m)  versus  the  velocity  structure  factor  Cv2. 

8.  Range  (m)  versus  the  acoustic  refractive  index 
structure  factor  Cne2. 

9.  Range  (m)  versus  the  power  backseat tered  to  the 
echosounder . 

All  of  the  programs  prompt  the  operator  for  inputs  and 
with  a  series  of  yes/no  questions  allows  the  operator  to  re¬ 
run  with  the  same  inputs  or  change  the  inputs.  Figure  30  is 
a  flow  chart  of  program  one  and  appendix  1  is  a  copy. 
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Figures  1,  2,  3,  4,  21,  22,  23,  24,  25,  26,  28,  and  29  are 
output  plots  of  program  one. 

The  second  program  takes  the  same  atmospheric  parameters 
used  in  the  first  program.  With  the  exception  of  the 
frequency,  it  takes  the  same  echosounder  parameters  also. 
This  program  outputs  a  plot  of  range  versus  frequency  and  a 
plot  of  range  versus  excess  attenuation  for  various 
frequencies.  Figure  31  is  a  flow  chart  of  the  program  and 
appendix  2  is  a  copy.  Some  output  plots  of  the  program  are 
included  in  the  conclusions  section. 

The  third  program  trkes  the  same  atmospheric  and 
echosounder  parameters  as  the  first.  The  program  ou touts  a 
plot  of  range  as  a  function  of  efficiency  of  the  transducer, 
assuming  the  same  efficiency  for  transmit  and  receive.  All 
the  other  programs  and  plots  in  this  thesis  assume 
efficiencies  of  25%,  which  is  on  the  high  side  of  typical 
performance.  Appendix  3  is  a  copy  and  an  output  plot  is 
presented  in  the  conclusions  section. 

The  fourth  program  has  the  same  inputs  as  the  first  with 
the  exception  of  antenna  area,  which  is  the  dependent 
variable  for  the  output  plot.  It  plots  the  range  as  a 
function  of  antenna  area  for  several  frequencies  and  is 
presented  and  discussed  in  the  conclusions.  The  model  does 
not  include  different  efficiencies  based  on  optimum  antenna 
design  for  different  frequencies.  The  output  reflects  the 


effect  of  different  antenna  diameters  on  the  excess 
attenuation  and  the  echosonde  equation.  I  am  uncertain  how 
much  these  drive  the  considerations  for  optimum  antenna 
design.  Appendix  4  is  a  copy  of  the  program. 

The  fifth  program  has  the  power  to  the  transducer  as  the 
dependent  variable  plotted  against  range  for  several 
frequencies.  Appendix  5  is  a  copy  and  the  output  plot  is 
presented  and  discussed  in  the  conclusions  section. 

The  sixth  program  has  the  ambient,  atmospheric 
background  noise  as  the  dependent  variable.  Appendix  6 
contains  a  copy  of  the  program  and  the  output  plot  is 
presented  and  discussed  in  the  conclusions  section. 

These  programs,  taken  together,  allow  a  parametric 
analysis  of  the  effect  of  different  parameter  changes. 


IV  CONCLUSIONS 


Reference  10  Ip.  II-l]  indicates  that  as  the  frequency 
increases : 

1)  Background  and  wind  noises  decrease  except  for  marked 
peaks  due  to  fans,  etc.  This  relationship  between 
noise  and  frequency  is  not  included  in  the  this  model. 

2)  The  reflectivity  of  the  turbulence  in  the  atmosphere 
decreases  slowly  and  the  absorption  of  the  atmosphere 
increases  rapidly.  The  model  agrees  with  this  for  low 
values  of  the  water-vapor  pressure.  The  effects  of 
water-vapor  pressure  on  the  attenuation  are  shown  in 
Figures  1  to  4 . 

3)  Wildlife  sounds  tend  to  increase  with  increasing 
frequency  and  dominate  the  background  noise  at  about 
3000  Hertz.  This  effect  is  not  included  in  the  model. 

4)  The  Doppler  shift  requires  the  receiver  bandwidth  to  be 
increased.  The  model  doesn't  take  this  into 
account.  The  bandwidth  is  used  to  calculate  the 
Johnson  noise  but  Johnson  noise  is  not  significant. 

The  noise  really  should  be  scaled  with  frequency  to 
depict  the  gain  in  range  accurately  for  lower 
frequencies . 


Transducer  efficiencies  vary  with  frequency.  This  is 
not  included  in  the  model.  An  increase  in  efficiency  will 
have  a  much  greater  impact  on  the  potential  range  than 
increases  in,  say,  power.  An  increase  in  efficiency  will 
help  both  with  transmission  and  return,  increasing  the 
transmitted  power  and  the  returned  electrical  signal 
strength.  Figures  32  is  a  plot  of  the  maximum  range  for 
efficiencies  of  .05  to  .5.  This  figure  is  the  output  of 
program  3.  The  range  increases  quickly  with  improvements  in 
efficiency.  The  efficiency  could  be  improved  by  using 
better  designed  horns,  such  as  catenoidal  or  exponential. 

Optimum  antenna  diameters  vary  with  frequency.  The 
output  of  program  4,  Figure  33,  demonstrates  this  effect. 

Not  included  are  the  effects  on  the  antenna  effective 
aperture  factor  G,  except  to  the  extent  it  may  be  effected 
by  the  excess  attenuation.  The  discontinuities  in  the 
curves  are  due  to  the  discontinuities  in  the  equations  for 
the  excess  attenuation.  For  a  given  frequency,  the  range 
increases  and  then,  as  antenna  diameters  increase  beyond  an 
optimum,  the  range  decreases.  The  decrease  for  larger 
antenna  diameters  is  the  result  of  the  excess  attenuation. 

One  other  consideration  as  to  choice  of  frequency  is 
that  the  resolution  is  increased  with  increasing  frequency. 

A  specific  need  for  the  echosounder  may  drive  this 
constraint  and  is  not  considered  in  this  model. 
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Fig.  32.  Range  versus  Transducer  Efficiency 
for  Profile  3,  Parameters  of  the  Aerovironment 
T  Model  300,  and  S.T.P. 
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Pig.  33.  Range  versus  Antenna  Diameter 
for  C_2  Profile  3,  Parameters  of  the  Aerovironment 
Model  300,  and  S.T.P. 


Figure  34  is  a  plot  of  range  as  a  function  of  frequency 
for  an  antenna  size  of  1.5  meters.  The  water-vapor  pressure 
input  was  ten  millibars.  Figure  35  is  for  the  same  antenna 
size  but  with  a  water-vapor  pressure  of  2  millibars.  The 
step  in  each  curve  is  not  a  physical  effect  but  rather  due 
to  the  step  of  1.5  in  the  theoretical  equation  for  the 
excess  attenuation.  Figure  1  showed  the  greater  value  of 
the  attenuation  at  around  2  mb.  water-vapor  pressure  for  the 
range  of  frequencies  we  are  dealing  with.  These  previous 
two  plots  demonstrate  the  effect.  Figures  34  and  35  are 
some  of  the  output  plots  of  program  2.  Both  figures 
demonstrate  the  greater  attenuation  of  the  acoustic  wave  as 
the  frequency  increases. 

Significant  increases  in  range  can  be  achieved  by 
judicious  choice  of  frequency.  Cons iderations  are  the 
water-vapor  pressure  (Figures  1  to  4),  the  excess 
attenuation  (Figures  5  to  7),  which  depends  on  the  antenna 
diameter  (Figure  33),  the  frequency  (Figure  33  to  35),  and 
the  anticipated  frequency  spectrum  of  the  background  noise. 

Figure  36  models  the  effect  of  increasing  the 
transmitted  power  to  increase  the  range.  This  is  the  output 
plot  of  program  5.  The  slope  of  the  curve  is  not  very  steep 
and  becomes  less  so  for  higher  frequencies.  A  considerable 
increase  in  power  is  required  to  double  the  range.  Also  the 
model  does  not  include  nonlinear  transducer  effects  which 


£  (METERS)  RANGE  (METERS 


1200 

1100 

1000 


600- 

500- 


FREQUENCY 


2500 


Fig.  34.  Range  versus  Echosounder  Frequency 
for  Profile  3,  Parameters  of  the  Aerovironment 
Model  300 ,  and  S.T.P. 
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Fig.  35.  Range  versus  Echosounder  Frequency 
for  C-,^  Profile  3#  Parameters  of  the  Aerovironment 
Model  300.  and  S.T.P.  except  a  Water-vapor  Pressure  of  2  mb 


arise  when  the  power  is  increased.  Reference  1  [p.  57] 
points  oat  that  at  some  point  the  acoustic  wave  becomes 
distorted,  which  implies,  from  Fourier  analysis,  there  is  a 
flow  of  energy  out  of  the  fundamental  frequency  into  the 
higher  harmonics.  Since  the  higher  frequencies  are 
attenuated  more  quickly  we  will  soon  reach  a  point  of 
saturation,  which  Brown  [Ref.  1  :p  57]  called  nonlinear 
saturation  to  distinguish  it  from  another  saturation  effect 
he  discusses. 

Figure  37  is  the  output  of  program  6.  This  shows  the 
effect  of  decreasing  the  ambient  background  noise  level. 

This  decrease  in  noise  might  be  achieved  by  using  digital 
processing  and  fast  Fourier  transforms  to  achieve  a  narrower 
bandwidth.  The  bandwidth,  and  therefore  the  noise,  can  be 
much  smaller. 

These  plots  demonstrate  which  parameter  changes  might 
best  improve  the  range  of  an  echosounder,  which  was  one  of 
our  goals  in  this  thesis.  The  modeled  performance  changes 
should  allow  for  intelligent  decisions  of  the  necessary 
parameters  for  the  expected  uses  of  the  echosounder. 

We  also  sought  to  explore  the  ability  to  use  the 
returned  signal  to  quantify  atmospheric  parameters 
accurately  at  a  given  range.  As  can  be  seen  from  Figure  27 
(the  input  flow  chart),  for  a  single  acoustic  radar  return 
you  could  calculate  CT^  based  on  an  assumed  profile  of  Cv^. 
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However  this  need  to  assume  a  Cvz  profile  could  be 

eliminated  with  the  use  of  two  or  more  echosounders . 

Reference  11  and  many  others  describe  techniques. 

Another  approach  would  be  to  measure  the  Doppler  width. 

The  Doppler  width  is  the  spread  of  frequency  around  the 

Doppler  shifted  frequency.  Epsilon  and  therefore  Cv2  can  be 

related  to  the  Doppler  width.  With  the  use  of  fast  Fourier 

2  2 

transforms  and  digital  processing  Cv  and  CT  could  be 
measured  simultaneously  with  one  echosounder. 

rising  the  returned  signal  from  lesser  ranges  the  energy 
incident  on  a  given  volume  could  be  estimated  and  the 
degradation  of  the  return  signal  could  be  estimated.  In 
this  way  the  computer  program  would  allow  one  to  essentially 
boot-strap  up  to  a  given  range  and  more  accurately  depict 
the  atmospheric  parameters  based  on  the  returned  signal. 
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APPENDIX  A 

COMPUTER  PROGRAM  1 

FULLER,  ROBERT  PROS  1 

10  SEP  85 


•••♦•••PURPOSE******* 

THIS  PROGRAM  WILL  TAKE  PARAMETERS  OF  AN  ECHOSOUNDER  AND 
ESTIMATE  THE  RAN6E.  THE  FOLLOWING  INPUTS  ARE  REQUIRED: 

ATMOSPHERIC  DATA 

1 1ATOMOSPHERIC  PRESSURE  in  millibars 

Z ) SELECT  ONE  OF  FOUR  ATMOSPHERIC  PROFILES  FOR  THE  TEMPERATURE 
STRUCTURE  PROFILE. 

a) FOR  ONE  PROFILE  THE  HEIGHT  OF  THE  INVERSION  LAYER  IS  ALSO  INPUT 

3)  TEMPERATURE  IN  DEGREES  CELSIUS 

4)  WATER  VAPOR  PRESSURE  IN  millibars 

ECHOSOUNDER  DATA 

5)  ANTENNA  DIAMETER  IN  METERS 

6) FREQUENCY  OF  ECHOSOUNDER  IN  Hz 

71P0WER  TRANSMITTED  BY  ECHOSOUNDER  IN  WATTS 
8)PULSE  LENGTH  OF  THE  TRANSMITTED  ACOUSTIC  ENERGY 

THE  PROGRAM  OUTPUTS  THE  FOLLOWING  PLOTS: 

1 ) ATTENUATION! 1/m. )  VERSUS  WATER-VAPOR  PRESSURE ( mb)  FOR  FIVE 
FREQUENCIES  AT  ONE-THIRD  OCTIVES  AROUND  THE  INPUT  FREQUENCY. 

Z) ATTENUATION  1/m. )  VERSUS  RELATIVE  HUMIDITY(X)  FOR  FIVE  FREQUENCIES 
AT  ONE- THIRD  OCTIVES  AROUND  THE  INPUT  FREQUENCY. 

3) ATTENUATION 1/n. )  VERSUS  WATER-VAPOR  PRCSSURE(mb)  FOR  FIVE 

TEMPERATURES  AT  TEN  DEGREE  INTERVALS  AROUND  THE  INPUT  TEMPERATURE. 
4>ATTENUATI0N( 1/m. )  VERSUS  RELATIVE  HUMIDITY(X>  FOR  FIVE  TEMPERATURES 
AT  TEN  DEGREE  INTERVALS  AROUND  THE  INPUT  TEMPERATURE. 

5)RANGE(m. )  VERSUS  EXCESS  ATTENUATION 
G)RANGE( m. )  VERSUS  TEMPERATURE  STRUCTURE  FACTOR. 

7>RANGE<  m. )  VERSUS  VELOCITY  STRUCTURE  FACTOR. 

8) RANGE<  m. )  VERSUS  ACOUSTIC  STRUCTURE  FACTOR. 

9 )  RANGE < m. >  VERSUS  BACKSCATTERED  TO  ECHOSOUNDER. 


VARIABLES 


Again  TEST  VALUE  IN  WHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  MAKE  ANOUTHER  RUN. 

Ano_change  TEST  VALUE  IN  WHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  CHANGE  THE  VALUE  OF  INPUT  DATA  BEFORE  MAKING 
ANOUTHER  RUN. 

Ant_area  ECHOSOUNDER  ANTENNA  AREA  CALCULATED  FROM  INPUT  OF 


a 

-V-TV7A  ■  .VV.V.V".* .  “  .*  '  ■  - ■»  V*  ?  ■■  ^  7 

! 

ANTENNA  DIAMETER. 

t  1 

Ant_dia« 

INPUT  OF  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS. 

!  i 

Aton_pres 

INPUT  OF  ATMOSPHERIC  PRESSURE  IN  nb. 

:•  i 

At  ten 

ATTENUATION  OF  ACOUSTIC  WAVE.  CALCULATED  IN 

•;  ! 

SUBPROGRAM  ATTENUATION. 

■I-  i 

Att_f req 

VARIABLE  USED  TO  DETERMINE  IF  OPERATOR 

!•  i 

i  ; 

! 

UANTS  TO  PLOT  ATTENUATION  VERSUS  UATER.VAPOR 

PRESSURE  FOR  VARIOUS  FREQUENCIES  AROUND  THE  INPUT 

FREQUENCY.  IF  SO  THEN  PLOT  IS  DONE  IN  SUBPROGRAM 

Att_f req. 

1 

Att_nax 

VARIABLE  IN  SUBPROGRAM  ATTENUATION.  IT  IS  THE 

j 

! 

ATTENUATION  AT  THE  FREQUENCY  OF  THE  MAXIMUM 

ATTENUATION  FOR  THE  INPUT  CONDITIONS. 

.  • 

Att_temp 

VARIABLE  USED  TO  DETERMINE  IF  OPERATOR  WANTS 

! 

! 

TO  PLOT  ATTENUATION  VERSUS  UATER.VAPOR  PRESSURE 

FOR  VARIOUS  TEMPERATURE  AROUND  THE  INPUT  TEMPERATURE. 

THE  PLOT  IS  DONE  IN  THE  SUBPROGRAM  Att_temp. 

Bn 

BANDWIDTH  OF  RECEIVER  NEEDED  TO  RECEIVE  SIGNALS 

! 

I 

DOPPLER  SHIFFED  BY  THREE  METER  PER  SECOND 

VERTICAL  VELOCITIES. 

c 

VELOCITY  OF  SOUND  CALCULATION  AT  ZERO  CELSIUS. 

®  ! 

Cl  &  cz 

VARIABLES  USED  TO  SCALE  INVERSION  HEIGHT  IN 

! 

TEMPERATURE  STRUCTURE  PROFILE  THREE. 

j 

C3 

CONSTANT  USED  IN  CALCULATION  FOR  CteZ  PROFILE  4. 

>•  ! 

CneZ< • ) 

ARRAY  OF  VALUES  OF  THE  ACOUSTIC  REFRACTIVE  INDEX 

PARAMETER.  CALCULATED  BASED  ON  SELECTION  OF 

i  ; 

! 

PROFILE  FOR  TEMPERATURE  STRUCTURE  PARAMETER 

AND  ASSUMED  PROFILE  FOR  THE  VELOCITY  STRUCTURE 

PARAMETER. 

! 

CteZ( • ) 

ARRAY  OF  VALUES  OF  THE  TEMPERATURE  STRUCTURE  PARAMETER. 

! 

THE  OPERATOR  SELECTS  A  PROFILE  FROM  SEVERAL  PROVIDED. 

•  ;  j 

CveZ( • ) 

ARRAY  OF  VALUES  OF  THE  VELOCITY  STRUCTURE  PARAMETER. 

i  ! 

VALUES  BASED  ON  CALCULATION  USING  ASSUMED  OIS5APATION 

RATE. 

Dx 

STEP  SIZE  FOR  X  AXIS  FOR  VARIOUS  PLOTS. 

j 

Dy 

STEP  SIZE  FOR  Y  AXIS  FOR  VARIOUS  PLOTS. 

j 

Epsilon 

DISSAPATION  RATE  USED  IN  CALCULATION  OF  CveZ. 

! 

Et 

TRANSMISSION  EFFICIENCY  OF  ECHOSOUNDER. 

I 

Er 

EFFICIENCY  OF  ECHOSOUNDER  WHEN  RECEIVING  ACOUSTIC 

B  1 

BACK SCATTER. 

r-  j 

Es 

SATURATION  VAPOR  PRESSURE  AT  GIVEN  TEMPERATURE. 

Exc_att< « ) 

EXCESS  "ATTENUATION”  AT  GIVEN  RANGE. 

F 

VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.  IS  THE 

1 

RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 

! 

Fnax 

FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 

i  i 

ATTENUATION. 

• 

Freq 

INPUT  FREQUENCY  OF  ECHOSOUNDER. 

! 

Freq_con 

USED  TO  INSURE  VALUE  OF  FREQUENCY  PASSED  TO  SUBPROGRAM 

■'  I 

Att_f req  WAS  NOT  CHANGED. 

»  . 

6 

ANTENNA  EFFECTIVE  APERATURE  FACTOR. 

! 
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H  VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

I  MAIN  INDEX  FOR  LOOP  TO  CHECK  RETURNED  POWER  UNTIL 

POWER  WAS  LESS  THAN  MINIMUM  DETECTABLE  SIGNAL. 

Interval  MAXIMUM  EFFECTIVE- SCATTERING  VOLUME  THICKNESS. 

Inver  .  HEIGHT  OF  INVERSION  LAYER. 

J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

K  WAVENUMBER 

L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

Mess_up  TEST  VALUE  TO  CHECK  RESPONSES  WHEN  OPERATOR  IS  INPUTTING 
RESPONSES. 

N  TEST  VALUE  IN  CALCULATION  OF  COHERENCE  LENTH  IN  CALCULATION 

OF  EXCESS  ATTENUATION. 

Neu_va  TEST  VALUE  IN  CHECKING  IF  OPERATOR  WISHES  TO  CHANGE 
A  VARIABLE  BEFORE  A  NEW  RUN. 

Noise  ASSUMED  MINIMUM  DETECTABLE  SIGNAL. 

Pou_back<  * )  POWER  BACK SCATTERED  FROM  GIVEN  RANGE. 

Pow_ret <  * )  POWER  BACKSCATTERED  TO  ECHOSOUNDER  FROM  GIVEN  RANGE. 

Pou_trans  INPUT  OF  POWER  SUPPLIED  TO  TRANSDUCER  OF  ECHOSOUNDER. 

Profile  OPERATOR  INPUT  OF  CteZ  PROFILE  FROM  AVAILABLE  PROFILES. 

Pstar  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 

CALCULATION. 

Pulse  TRANSMITTED  PULSE  LENGTH  IN  MILLISECONDS.  INPUT  BY 

OPERATOR. 

R  RANGE  FOR  INPUT  TO  INTEGRAL  OF  EXCESS  ATTENUATION. 

R1  VARIABLE  USED  IN  SCALING  INVERSION  HEIGHT  FOR  PROFILE  THREE. 

Rangel*)  ARRAY  OF  RANGE  VALUES. 

Remainder  REMAINDER  OF  MOOULO  FUNCTION  USED  TO  DECREES  THE 

NUMBER  OF  PASSES  THROUGH  THE  INTEGRATION  FOR  EXCESS 
ATTENUATION. 

Rge  VALUE  OF  RANGE  IN  INTEGRAL  FOR  EXCESS  ATTENUATION. 

Rho  CORRELATION  LENGTH  USED  IN  CALCULATION  OF  EXCESS 

ATTENUATION. 

Signal • )  FRACTION  OF  ENERGY  BACKSCATTERED  FROM  GIVEN  RANGE. 

Speed_sound  SPEED  OF  SOUND  IN  AIR  AT  INPUT  TEMPERATURE. 

Sumpou.back  SUM  OF  THE  BACKSCATTERED  ENERGY 

T  INPUT  TEMPERATURE  IN  DE6REES  KELVIN. 

Temp  INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

Tenp_con  USED  TO  INSURE  VALUE  OF  TEMPERATURE  PASSED  TO  SUBPROGRAM 
Att_tenp  WAS  NOT  CHANGED. 

Temp*  VARIABLE  STRING  USED  IN  FUNCTION  YES. 

Title*  STRING  PASSED  TO  SUBPROGRAM  PT  FOR  TITLE  OF  PLOT. 

Tstar  INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

IN  SUBPROGRAM  ATTENUATION. 

Var  USED  TO  SELECT  THE  VARIABLE  THE  OPERATOR  WISHED  TO 

CHANGE  BEFORE  MAKING  ANOTHER  RUN. 

Wat_pres  ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.  INPUT  BY 
OPERATOR. 

X  THIRD  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 

X*  STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 


I  RESPONSE  TO  YES  OR  NO  QUESTION. 

!  XlabelS  LABEL  ON  X  AXIS  PASSED  TO  SUBPROGRAM  Pt  FOR  PLOTTING. 

1  Xnax  VALUE  OF  LARGEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

I  Xnin  VALUE  OF  SMALLEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

!  Xrange  VALUE  OF  RANGE  OF  X  VALUES  FOR  VARIOUS  PLOTS. 

!  Xvar< * )  VALUE  OF  ARRAY  OF  X  VALUES  TO  BE  PLOTTED  IN  SUBPROGRAM  Pt. 

I  YlabelS  LABEL  ON  Y  AXIS  PASSED  TO  SUBPROGRAM  Pt  FOR  PLOTTING. 

I  Ynax  LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

I  Ynin  SMALLEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

!  Yrange  RANGE  OF  Y  VARIABLES  FOR  VARIOUS  PLOTS. 

!  Ze  INTERMEDIATE  VALUE  IN  CALCULATION  OF  EXCESS  ATTENUATION. 


! DECLARE  VARIABLES 
INTEGER  I 

DIM  Pou_back< 1500) ,Pou_ret< 1500) .Rangel  1500) ,CteZ< 1500) .Signal  1500) 
DIM  Cve2( 1500) ,CneZ< 1500) ,Exc_att(  1500) ,Xvar(  1500) 

DIM  Title*! 501  ,XlabelStS0]  .YlabelSUGl 
! 

PLOTTER  IS  705 , “ HPGL " 

LINE  TYPE  1 
! 

! INPUT  ATMOSPHERIC  DATA 

INPUT  "ENTER  TEMPERATURE  IN  DEGREES  CELCIUS" .Temp 
INPUT  "ENTER  ATMOSPHERIC  PRESSURE  IN  MILLIBARS" ,Aton_pres 
INPUT  "ENTER  UATER  VAPOR  PRESSURE  IN  MILLIBARS" . Uat_pres 
! 

! 

! INPUT  ECHOSOUNOER  DATA 

INPUT  "ENTER  FREQUENCY  OF  ECHOSOUNOER  IN  HERTZ", Freq 
INPUT  "ENTER  ANTENNA  DIAM  IN  METERS" ,Ant_di an 
INPUT  "ENTER  ECHOSOUNOER  PULSE  LENGTH  IN  MILLISECONDS-" , Pul se 
INPUT  "ENTER  POUER  TO  TRANSMITTER  IN  WATTS" ,Pou_trp is 
Et-.Z5  ! TRANSMIT  EFFICIENCY 

Er-.Z5  ! RECEIVER  EFFICIENCY 

G-.40  1  ANTENNA  EFFECTIVE  APERATURE  FACTOR 

! 

! 

I  SELECT  THE  PROFILE  FOR  THE  TEMPERATURE  STRUCTURE 

PRINT  "YOU  HAVE  A  CHOICE  OF  TEMPERATURE  STRUCTURE  PROFILES" 

PRINT  "PROFILE  TO  BE  USED." 

PRINT  "YOUR  SELECTIONS  ARE" 

PRINT  "  1  A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS  " 

PRINT  "  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  UHICH  " 

PRINT  "  GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 

PRINT  "  " 

PRINT  "  Z  THE  TEMPERATURE  STRUCTURE  PROFILE  ABOVE  BUT" 

PRINT  ”  WITH  A  FACTOR  OF  TWO  TO  APPROXIMATE  LOOKING" 

PRINT  "  UP  A  CONVECTIVE  PLUME" 

PRINT  "  " 
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PRINT  “  3  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 

PRINT  "  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 

PRINT  "  NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH  HEIGHT  " 

PRINT  “  IS  EXP<-.001*HEIGHT  ABOVE  BS  METERS  AND  HEIGHT" 

PRINT  "  TO  THE  -1.4S  BELOU  G5  METERS" 

PRINT  "  " 

PRINT  "  4  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 

PRINT  "  HEIGHT  TO  THE  -4/3  " 

! 

PRINT  ”  “ 

PRINT  “  " 

Mees_up- 1 
UHILE  Mees_up-1 

INPUT  "ENTER  THE  DESIRED  PROFILE  (  1  OR  2  OR  3  OR  4)". Profile 
IF  Profile-4  THEN 

INPUT  “ENTER  THE  HEIGHT  OF  THE  INVERSION  IN  METERS" . Inver 
END  IF 

IF  Prof  lie- 1  THEN 
Mese_up-0 

ELSE 

IF  Profile-2  THEN 
Mess_up-0 

ELSE 

IF  Profile-3  THEN 
Mess_up“0 

ELSE 

IF  Profile-4  THEN 
Mess_up-0 

ELSE 

PRINT  Profile,”  UAS  NOT  ONE  OF  THE  OPTIONS!!!! 
Mess_up- 1 
END  IF 
END  IF 
END  IF 
END  IF 
END  UHILE 
OUTPUT  KBD; "K " ; 


Again- 1 

! INITIALIZE  ARRAYS  FOR  SUCCESSIVE  RUNS 
UHILE  Again- 1 
FOR  J-1  TO  I 

Pou_back( J)-0 
Pou_ret<  J)-0 
Range <  J>-0 
CteZ< J>-0 
Signa( J)-0 
CveZ( J  >-0 


Exc_att(  J)-0 
NEXT  J 

FOR  J-1  TO  15 
PRINT  "  " 

NEXT  J 

PRINT  "DO  YOU  UISH  TO  SEE  A  PLOT  OF  ATTENUATION  VERSUS  HUMIDITY" 

PRINT  "FOR  VARIOUS  FREQUENCIES  AROUND  THE  FREQUENCY  YOU  INPUT?" 

LINPUT  "IF  YES  ENTER  Y  .  IF  NO  ENTER  N",X$ 

Att_f req-FNYes< X$) 

IF  Att_freq-1  THEN 
OUTPUT  KBD!“K"j 

CALL  Att_f req<  At on_pr os , Freq , T enp  > 

END  IF 

OUTPUT  KBO; "K" i 
! 

! 

FOR  J-1  TO  15 
PRINT  “  " 

NEXT  J 

PRINT  "DO  YOU  UISH  TO  SEE  A  PLOT  OF  ATTENUATION  VERSUS  HUMIDITY" 

PRINT  "FOR  VARIOUS  TEMPERATURES  AROUND  THE  TEMPERATURE  YOU  INPUT?" 

LINPUT  "IF  YES  ENTER  Y  ,  IF  NO  ENTER  N‘*,XS 
Att_tenp-FNYes( X$) 

IF  Att_tenp-1  THEN 
OUTPUT  KBDi "K" s 

CALL  Att_tenp<  Ato« _pres .Freq.Tenp) 

END  IF 

OUTPUT  KBD; “K" 5 
! 

! CONVERT  TEMPERATURE  TO  KELVIN 
T-Tenp+Z73 

! 

! CALCULATE  SPEED  OF  SOUND  BASED  ON  INPUT  TEMPERATURE( CELCIUS ) 
Speed_sound-Z0.05«<  T)“ .5 

! CALCULATE  THE  SPEED  OF  SOUND  AT  0  DEGREES  CELCIUS 
C-Z0. 05*Z73“ . 5 
( 

•CALCULATE  THE  ATTENUATION  COEFFICIENT  BASED  ON  ATMOSPHERIC  DATA 
•EQUATION  IS  FROM  NEFF  1975 

CALL  Rttenuat lon( Ato«_pre5 , At ten.Freq , Temp , Uat_pre® ) 

1-1  ! INDEX  FOR  LOOP 

Pou_back(0)-0  UNITIALIZE  VARIABLE  FOR  POWER  BACKSCATERRED 
K-Z*PI»Freq/Speed_sound  • WAVENUMBER 
Ant_area-PI»( Ant_dian/Z)*Z  ! ANTENNA  AREA 
Interval -( Speed_sound«Pul se« 1 . E- 3 ) / Z 

Bn-Z*Freq*< 1  -  1 /( 6/Speed_oound+ 1 ) ) ! BANDWIDTH  FOR  3M/S  VERTICAL  VELOCITIES 
Noiae-1 . 38E-Z3*Bn»< T1+Z.E- 14! MINIMUM  DETECTABLE  SIGNAL  AS 
!  JOHNSON  NOISE( NEGLIGABLE )+  ESTIMATED  BACKGROUND 

! 

! CALCULATE  RETURNED  SIGNAL  POWER  UNTIL  IT  IS  LESS  THAN  BACKGROUND. 
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T 


• 1  •'  ■  ».» ■  ■  v  i;  i;,v.ri.T 


y  .■*  - 


If 


Range  ( 0)-0 
Sumpou_back “0 
REPEAT 

Range < I ) -Ranget I- 1 )+Z 
SELECT  Profile 
CASE  1 

I  THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
IAS  PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398. 

! WAS  FOR  MIDDAY  ANO  CLEAR  WEATHER  ABOVE  THE 
I TULAROSA  BASIN,  NEW  MEXICO. 

! AN  ADDITONAL  FACTOR  OF  1 . 8  WAS  INCLUDED  AS  PER  DIFFERENT 
! AVERAGING  TIME. 

CteZ< I )«2. 12«Range( I >A<-1 . IB) 

CASE  Z 

•THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
! AS  PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398. 

•WAS  FOR  MIDDAY  AND  CLEAR  WEATHER  ABOVE  THE 
! TULAROSA  8ASIN,  NEW  MEXICO. 

! AN  ADDITONAL  FACTOR  OF  1 . 8  WAS  INCLUDED  AS  PER  DIFFERENT 
! AVERAGING  TIME  AND  AN  ADDITIONAL  FACTOR  OF  TWO  WAS 
•INCLUDED  TO  APPROXIMATE  LOOKING  UP  A  CONVECTIVE  PLUME 
CteZd >-Z«Z. 1 Z*Range< I >*(-1 . 1G> 

CASE  3 

•THIS  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 

! PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398  FOR" 

•NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH  HEIGHT  " 

•IS  EXP( - . 001 *H£IGHT  ABOVE  65  METERS  AND  HEIGHT" 

•TO  THE  -1.46  BELOW  65  METERS" 

IF  Ranged) <65  THEN 

CteZ(  I ) -75. 5*Rar»ge<  I ) * ( - 1 . 46 ) 

ELSE 

CteZ( I >-3.66E-Z*EXP< - .001«Range< I ) ) 

END  IF 
CASE  4 

•CALCULATE  THE  TEMPERATURE  STRUCTURE  FACTOR  BASED  ON  ASSUMPTION 
!  THAT  IT’S  PROPORTIONAL  TO  HEIGHT  TO  THE  -4/3  POWER 

•  THAT  THE  SURFACE  VERTICAL  HEAT  FLUX  IS  .095  C.n./sec 

•THE  DEPENDENCE  WITH  HEIGHT  ABOVE  THE  INVERSION  I  PULLED  FROM 
•MY  LEFT  EAR 

•EQUATION  FROM  NEFF,  1975 
SELECT  Range< I ) /Inver 
CASE  <.9 

C3-<< . 0Z4>«<  T)“( .667)) 

CteZ < I )-C3*< Rangel  I ) )*< - 1 .33) 

CASE  .9  TO  1 

CteZ< I >-CteZd-l ) 

R1 -Range< I > 

Cl-CteZd) 

CASE  1  TO  1.3 

Cte2< I )-10"<  <LGT(C3*Range< I ) )-LGT( C3«R1 ) )+LGT(C1 ) ) 


> 

t 

) 


R2«Range(  I ) 

C2«CteZ(I> 

CASE  ELSE 

CteZ< I )-C3*<  Range< I >*< - 1 . 33>-R2*< - 1.33) )+CZ 
END  SELECT 
ENO  SELECT 

! CALCULATE  VELOCITY  STRUCTURE  FACTOR.  FORMULA  FROM  6AYN0R  77 
(AND  TETARSKI 

Epsilon-< .2 856/ Range ( I ) )»<  1 +.07*< Range( I >  )*< .5) >*< 1.5) 

Cve2< I )“2*Epsilon“< . 557) 

! CALCULATE  ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  FACTOR. 

! FORMULA  FROM  TETARSKI  1961 
Cne2(I)-(Cte2(I)/(2.98E+5))+<CveZ(I)/<C*C>) 

! CALCULATE  THE  FRACTION  OF  POWER  BACKSCATTERED  FROM  INTERVAL. 

1  EQUATION  FROM  NEFF.  1975 
Signal I )■( .0039*1 K*< 1/3) )*CieZ< I )  )/<  T) *2 
! 

(CALCULATE  THE  EXCESS  ATTENUATION.  EQUATION  FROM  CLIFFORD,  1980 
(THE  EXCESS  ATTENUATION  IS  Exc_a-ti(I) 

(THE  MODULO  STRUCTURE  IS  TO  SKIP  SOME  OF  THE  INTEGRALS  ONCE  THE 
(RESOLUTION  IS  LESS  IMPORTANT. 

IF  I <50  THEN 

Renainder-0 

ELSE 

IF  I < 1 50  THEN 

Renainder-Range< I )  MODULO  10 

ELSE 

Renainder«Range( I >  MODULO  20 
END  IF 
END  IF 

IF  Renainder-0  THEN 
Rho-0 
L“0 
H-0 

Rge-Range(I)  (CONSTANT  IN  INTEGRAL 
R-0 

FOR  J-0  TO  2*1 

F-Cne2(INT(J/Z+J)> 

F-F»<  1  -R/(  Rge  )  *(  1 . 57  >+<  R/<  Rge ) )  *  <  1 . 57  > ) 

IF  J>0  THEN 

IF  J<Z*I  THEN 

IF  INT(J/Z)-J/2  THEN 
L-L+F 
F=0 

ELSE 

H-H+F 
F-0 
END  IF 
END  IF 
END  IF 
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Rho-Rho+F 

R-R+1 


NEXT  J 

Rho-Rho+4*L+Z*H 

Rho-(<<Rho».33)*K*K*1 .4G)A(-.B)) 

N-<  Ant_diam/Rho)A2 
IF  N<-1  THEN 

Ze-1/< 1+N) 

ELSE 

Ze-1 .5/< 1+N) 

ISTEP  OF  1.5-->SEE  CLIFFORD  1980 
END  IF 

Exc^attl  I  >-Ze*Ze 

ELSE 

Exc_att< I )-Exc_att< 1-1 ) 

END  IF 

| 

! CALCULATE  THE  POUER  BACK SCATTERED 

Pou_back( I )-( Poui_trans*Et-Sumpou_back )*EXP( -Atten*Range( I ) ) 

Pou_back< I )-Pou_back< I )*Interval*Exc_att( I )*Signa< I ) 
Sunpou_back-Sunpou_back+Pou_back( I > 

! CALCULATE  THE  POUER  BACK SCATTERED  TO  THE  ANTENNA 
! TUO  WAY  PATH  ATTENUATION  ACCOUNTED  FOR  ABOVE. 

Pou_ret( I )“Pou_back< I )*EXP( - At ten* Range < I ) >*Ant_area*6*Er/Range(  I )AZ 
PRINT  "RANGE*5"  .Range <  I ) 

PRINT  "POUER  RETURNED-" fPou_ret< I ) 

I-I  +  1 

UNTIL  Pou_ret(  1-1 > >0  AND  Pou_ret< I- 1 X-Noise 
! 

OUTPUT  KBD; "K" ; 

! 

PRINT  "INPUT  CONDITIONS" 

| 

PRINT  H  " 

PRINT  USING  "K"; "TEMPERATURE-  " , Tenp , "CELCIUS" 

PRINT  USING  "K"s  "ATMOSPHERIC  PRESSURE-  ".Atom  pres,  "rib" 

PRINT  USING  "K"  s  "UATER  VAPOR  PRESSURE-  "  ,Ua-t_pres ,  "mb" 

PRINT  USING  "K"; "PULSE  LENGTH-  " .Pulse, "ns" 

PRINT  USING  "K"; "TRANSMITTED  FREQUENCY-  ",Freq,"  Hz" 

PRINT  USING  "K" i "ANTENNA  DIAMETER-  " ,Ant_dian,"  n." 

PRINT  USING  "K" s "POUER  TRANSMITTED-  " ,Pou_trans,"  UATTS" 

PRINT  "TEMPERATURE  STRUCTURE  PROFILE  USED  ".Profile 
IF  Profile-4  THEN 

PRINT  "INVERSION  HEIGHT  ".Inver 
END  IF 

SELECT  Profile 
CASE  1 

PRINT  "  1  A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA" 

PRINT  "  AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  " 


PRINT  ”  WHICH  GIVES  ft  HEIGHT  TO  THE  -1.16  PROFILE" 

CASE  Z 

PRINT  "  2  A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA 

PRINT  "  AS  PRESENTED  IN  WALTERS/ KUNDEL  1981  PAGE  398 

PRINT  **  WHICH  GIVES  A  HEI6HT  TO  THE  -1.16  PROFILE" 

PRINT  "  BUT  WITH  A  FACTOR  OF  TWO  TO  APPROXIMATE" 

PRINT  “  LOOKING  IF  A  CONVECTIVE  PLUME" 

CASE  3 

PRINT  “  3  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA 

PRINT  ”  AS  PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398" 

PRINT  "  FOR  NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH" 

PRINT  "  HEIGHT  IS  EXP< -.001 ‘HEIGHT  ABOVE  65  METERS" 

PRINT  "  AND  HEIGHT  TO  THE  -1.46  BELOW  65  METERS" 

CASE  4 

PRINT  "  4  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 

PRINT  “  HEIGHT  TO  THE  -4/3  " 

END  SELECT 
I 

PRINT  USING  “K" » "THE  MTNIMUM  DETECTABLE  SIGNAL  WAS  SET  AT" .Noise, "WATTS 
! 

PRINT  “  " 

PRINT  "OUTPUT  CONDITIONS" 

! 

PRINT  USING  "K" ; "RANGE-" ,Range< I-Z ) n." 

! 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

OUTPUT  KBD; "K" ; 

PRINT  "THE  FOLLOWING  GRAPHS  WILL  NOW  BE  PLOTTED  WITH  RANGE  VERSUS" 

PRINT  "EXCESS  ATTENUATION" 

PRINT  "CTEZ  <  TEMPERATURE  STRUCTURE  PARAMETER)" 

PRINT  "CVEZ  <  VELOCITY  STRUCTURE  PARAMETER)" 

PRINT  "CNEZ  (ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  PARAMETER)" 

PRINT  "POWER  RETURNED" 

PRINT  "WHEN  YOU  ARE  READY  FOR  THE  FIRST  GRAPH  HIT  CONTINUE" 

PAUSE 

OUTPUT  KBD; "K" i 
! 

GRAPHICS  ON 
VIEWPORT  15,120,10,70 


! PLOT  RANGE  VERSUS  EXCESS  ATTENUATION 

Xmin-0 

Xmax-1 

Ynin»0 

IF  Range< I-Z) <-1000  THEN 
Ymax-3 
ELSE 


END  IF 

Dx-.l 

Dy-1 

XranQe-Xnax-Xnin 

Yrange-Ynax-Ynin  r> 

UINOOU  Xnln.Xnax ,Ynin, Ynax 
AXES  Dx ,Dy .Xnin, Ynin, 1 . 1 

CLIP  OFF  , 

! 

! LABEL  PLOT 
CSIZE  4, .6 
LOIR  0 
LORG  5 

MOVE  .5, Ynax* 1.1 

LABEL  "RANGE  VERSUS  EXCESS  ATTENUATION" 

1  LABEL  HORIZONTAL  AXES 
LDIR  0 
LORG  5 

FOR  J-0  TO  Xnax  STEP  Dx 
MOVE  J,-. 1*Dy 
LABEL  J 
NEXT  J 

MOVE  . 5*Xrange,- . 3*Oy 
CSIZE  4 , . G 

LABEL  "EXCESS  ATTENUATION" 

(LABEL  VERTICAL  AXES 
LORG  8 

FOR  J-0  TO  Ynax  STEP  Dy 
CSIZE  4 , . G 
MOVE  - . 3*Dx , J 
LABEL  "10“ 

CSIZE  Z 

MOVE  1*Dx,J+.05*Dy 
LABEL  USING  "K"|J 
NEXT  J 
LDIR  PI/Z 
LORG  G 

MOVE  -Dx , ,S*Yrange 
CSIZE  4, .6 

LABEL  "RANGE  (METERS)” 

CLIP  ON 

FOR  J-1  TO  1-1 

PLOT  Exc_att< J) ,LGT< Range< J ) ) 

NEXT  J 

PRINT  "HIT  CONTINUE  FOR  NEXT  PLOT" 

! 

I 

(PLOT  RANGE  VERSUS  TEMPERATURE  STRUCTURE  PARAMETER 

PAUSE 

GCLEAR 
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IF  Range!  I-ZX-1000  THEN 
Xnin— 4 
Ynax-3 
ELSE 

Xnin— 5 
Ynax-4 
ENO  IF 
Xnax-0 
Ynin-0 

Title*- "RANGE  vs  TEMPERATURE  STRUCTURE  PARAMETER" 

Xlabel*-“TEMPERATURE  STRUCTURE  PARAMETER" 

Ylabel*-" RANGE  < METERS >“ 

CALL  Log  loo< Xnin.Xnax .Ynin, Ynax , CteZ( • > ,  I- 1 , Range ( t ) , I- \ .Title*, Xlabel* 
i) 
j 
! 

I  PLOT  RANGE  VERSUS  VELOCITY  STRUCTURE  PARAMETER 

PAUSE 

GCLEAR 

I 

Xmin—  2 

Xnax-0 

Ynin-0 

IF  Rangel  1-2) <-1000  THEN 
Ynax-3 
ELSE 

Ynax-4 
END  IF 

T it le*-" RANGE  VERSUS  VELOCITY  STRUCTURE  PARAMETER" 

Xlabel*-" VELOCITY  STRUCTURE  PARAMETER" 

CALL  Log_log( Xnin.Xnax .Ynin, Ynax ,Cve2(  * ) , I- 1 .Range! * ) , I- 1 .Tit let, Xlabel* 
I) 

! 

! 

(PLOT  RANGE  VERSUS  ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  PARAMETER 

PAUSE 

GCLEAR 

Xnin--7 

Xnax— 5 

Ynin-0 

IF  Rangel  I -2) <-1000  THEN 
Ynax-3 
ELSE 

Ynax-4 
END  IF 

Title*-"RANGE  vs  ACOUSTIC  INDEX  STRUCTURE  PARAMETER” 

Xlabel*- "ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  PARAMETER" 

CALL  Loo  log! Xnin.Xnax , Ynin, Ynax .CneZl  • ) ,1- 1 .Rangel  • ) , I- 1 .Title*, Xlabel* 


§ 


I 


. 

>»• 


I 

I 

IPLOT  RANGE  VERSUS  POWER  RETURNED  TO  ANTENNA 

PAUSE 

GCLEAR 

IF  Ranged -ZX- 1000  THEN 
Xnin*-15 
Xnax“-5 

Ynax“3 

ELSE 

Xnin*-17 
Xnax“-4 
Ynaxm4 
END  IF 
Ynin*0 

Title*-"RANGE  VERSUS  POWER  RETURNED  TO  ANTENNA" 

Xlabel«-"POWER  RETURNED  TO  ANTENNA" 

Y label ••"RANGE  (METERS)" 

CALL  Log_log(  Xnin.Xnax ,Ynin, Ymax ,Pou_ret( *  > ,1-1 ,Range( * ) ,1-1 ,Title$,Xlab 
YlabelS) 

PAUSE 

GCLEAR 

! 

I 

! 

LINPUT  "WOULD  YOU  LIKE  TO  MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?",X» 
Again-FNYes( X*) 

Meas_up“1 
WHILE  Mess_up»1 
SELECT  Again 
CASE  1 

LINPUT  "WOULD  YOU  LIKE  TO  CHANGE  A  VARIABLE  (ENTER  Y  OR  N)?",X$ 
Neu_va“FNYes(  XS) 


Ano_change» 1 
WHILE  Ano_change“l 
SELECT  New  va 


CASE  1 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


USING 

USING 

USING 

il  ii 

USING 

USING 

USING 

USING 

USING 


VARIABLE  CURRENT  VALUE" 

"K'VI  TEMPERATURE  "  ,Tenp, "CELSIUS" 

"K'VZ  ATOMOSPHERIC  PRESSURE  "  ,  Atom_pres .  "mb 
"K" i "3  WATER  VAPOR  PRESSURE  " ,Wat_pres , "mb" 

"K" i "4  FREQUENCY  OF  ECHOSOUNDER  ",Freq,"  Hz" 
"K" ; "5  ANTENNA  DIAMETER  " ,Ant  dia«,"  n.  " 
"K"«"G  PULSE  LENGTH  ".Pulse,"  ns" 

"K" i "7  POWER  TRANSMITTED  ",Pou  trans,"  WATTS 
" K " i " 8  ATMOSPHERIC  PROFILE  ".Profile 


"ENTER  THE  NUMBER  OF  THE  VARIABLE  YOU  WISH  TO" 
"CHANGE" 
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INPUT  Var 
SELECT  Var 

CASE  1 

INPUT  "TEfFERATURE-" .Temp 
CASE  2 

INPUT  “ATMOSPHERIC  PRESSURE  IN  mb-” , Atom_pres 
CASE  3 

INPUT  "WATER  VAPOR  PRESSURE  IN  mb-" .Wat.pree 
CASE  4 

INPUT  "FREQUENCY  IN  Hz-",Freq 
CASE  S 

INPUT  "ANTENNA  DIAMETER  IN  m. , Ant_diam 
CASE  E 

INPUT  "PULSE  LEN6TH  IN  ma-", Pulse 
CASE  7 

INPUT  "POWER  TRANSMITTED  IN  WATTS-" ,Pou_trana 
CASE  8 

PRINT  "  1  -->  Ct*Z  PROFILE  OF  Z"<-1.16>" 

PRINT  "  FROM  WALTERS/KUNDEL  1981" 

PRINT  "2  — >  SAME  AS  ONE  BUT  WITH  FACTOR" 

PRINT  "  OF  TWO  TO  APPROXIMATE  LOOKING" 

PRINT  "  UP  A  CONVECTIVE  PLUME" 

PRINT  "  3  — >  Ct *2  PROFILE  OF  EXP(-.0O1*Z)" 

PRINT  M  FROM  WALTERS/KUNDEL  1981." 

PRINT  "  4  — >  CT*Z  PROFILE  OF  Z*C-4/3)" 

INPUT  "ENTER  NUMBER  OF  DESIRED  PROFILE" .Prof ile 
IF  Profile-4  THEN 

INPUT  "HEIGHT  OF  INVERSION  IN  METERS-" . Inver 
END  IF 
CASE  ELSE 

PRINT  Var. "IS  NOT  ONE  OF  THE  OPTIONS" 

END  SELECT 

LINPUT  "MADE  ANOTHER  CHANGE1  ENTER  Y  OR  N)?",X* 
Neu_va-FNYee(X*> 

Mesa_up-Z 
CASE  2 

Meaa_up-2 
Ano_change-Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  WHERE  Y  OR  N  I  111 ! " 

Ano_change- 1 

LINPUT  "CHANGE  A  VARIABLE  (ENTER  Y  OR  N)?",X* 
Neu_va-FNYea<  X$) 

END  SELECT 
END  WHILE 
CASE  2 

Mesa_up-Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  WHERE  Y  OR  N  !!!!!!” 

Mess_up- 1 
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LINPUT  "MAKE  ANOTHER  RUN  (ENTER  Y  OR  N>?" ,X* 
Again-FNYee<X$> 

ENO  SELECT 
ENO  WHILE 
ENO  WHILE 

PRINT  "THAT’S  ALL,  FOLKS" 
end 


ICALCULTE  THE  ATTENUATION 

eub  Attenuat ion<  Ato«_pree , Atten , Freq , Temp , Wat_pree ) 
(THIS  SUBPROGRAM  CALCULATES  THE  ATTENUATION  OF  SOUND 
!  IN  AIR  BASED  UPON  EQUATIONS  IN  NEFF  IS75 


INPUT 


ATMOSPHERIC  PRESSURE  IN  MILLIBARS 
FREQUENCY  OF  SOUND  WAVE  IN  HERTZ 
TEMPERATURE  IN  DEGREES  CELCIUS 
WATER-VAPOR  PRESSURE  IN  MILLIBARS 


OUTPUT  ATTENUATION  IN  1 /METERS 


(VARIABLES 

!  Atom_prefi  INPUT  OF  ATMOSPHERIC  PRESSURE  IN  nb. 

I  Atten  ATTENUATION  OF  ACOUSTIC  WAVE.  CALCULATED  IN 

!  SUBPROGRAM  ATTENUATION. 

I  Att_max  VARIABLE  IN  SUBPROGRAM  ATTENUATION.  IT  IS  THE 
!  ATTENUATION  AT  THE  FREQUENCY  OF  THE  MAXIMUM 

!  ATTENUATION  FOR  THE  INPUT  CONDITIONS. 

!  F  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.  IS  THE 

!  RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 

!  Fnax  FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 

I  ATTENUATION. 

!  Freq  INPUT  FREQUENCY  OF  ECHOSOUNDER. 

!  H  VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

!  Petar  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 

!  CALCULATION. 

!  Temp  INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

!  Tetar  INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

I  IN  SUBPROGRAM  ATTENUATION. 

!  Wat_pree  ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.  INPUT  BY 

!  OPERATOR. 

H- 1 00» Wat _pr es / At on_pr e s 
Tstar-< 1 .5»Temp+49z7/519 
Pstar»Aton_pres/ 10)4 

Fnax-< 10+Gl00*H+44400*H«H>*Pstar/Tstar‘\8 
Att_max«.0078*Fnax» Tatar* < -Z.5>*EXP<7.77*< 1-1/Tetar ) > 

F"Freq/Fmax 

Atten*(  Att_nax/304. 8 )•< < . 1 8*F )"Z  +  <  Z»F*F/( 1 +F*F ) >*Z  )* . S 

Atten-<  Atten+1 .74E- 10»Freq*Freq)/4. 35 

SUBEND 
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- \V-.\V  -.-.V'.- 


.'.v.v.'V. 


def  FNYes(Xt) 

!  THIS  FUNCTION  INTERPRETS  THE  OPERATORS  RESPONSES  TO  YES  NO  QUESTIONS 


I  INPUT  X* 

I  OUTPUT  FNYes 

t VARIABLES 

!  Tempt  VARIABLE  STRIN6  USED  IN  FUNCTION  YES. 

I  X*  STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

!  RESPONSE  TO  YES  OR  NO  QUESTION. 

DIM  Tempt! 1 ] 

Tempt! 1 ,11-TRIMt! Xt) 

SELECT  Tempt 
CASE  "Y” , "y" 

RETURN  I 
CASE  "N" .“n" 

RETURN  2 
CASE  "  " 

RETURN  I 
CASE  ELSE 

RETURN  -2 
END  SELECT 

FNEND 

! 

sub  Att_freq( A-tom_pres.Freq.Temp) 

I 

!  THIS  SUBPROGRAM  PLOTS  THE  ATTENUATION  OF  SOUND  VERSUS  UATER- 
! VAPOR  PRESSURE  FOR  FIVE  DIFFERENT  FREQUENCIES  AT  1/3  OCTIVE 
! INTERVALS  AROUND  THE  INPUT  FREQUENCY. 

! 

!  INPUT  ATOMOSPHERIC  PRESSURE  IN  MILLIBARS 
I  FREQUENCY  IN  HERTZ 

!  TEMPERATURE  IN  CELCIUS 

I 

!  OUTPUT  PLOT  OF  ATTENUATION  VERSUS  WATER- VAPOR  PRESSURE 
! 

! VARIABLES 

!  Atom_pres  INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 

!  A-tten  ATTENUATION  OF  ACOUSTIC  WAVE.  CALCULATED  IN 

!  SUBPROGRAM  ATTENUATION. 

!  Att.freq  VARIABLE  USED  TO  DETERMINE  IF  OPERATOR 
!  WANTS  TO  PLOT  ATTENUATION  VERSUS  WATER_VAPOR 

f  PRESSURE  FOR  VARIOUS  FREQUENCIES  AROUND  THE  INPUT 

!  FREQUENCY.  IF  SO  THEN  PLOT  IS  DONE  IN  SUBPROGRAM 

!  Att_freq. 

!  Es  SATURATION  VAPOR  PRESSURE  AT  GIVEN  TEMPERATURE. 

!  Freq  INPUT  FREQUENCY  OF  ECHOSOUNDER. 

I  Freq_con  USED  TO  INSURE  VALUE  OF  FREQUENCY  PASSED  TO  SUBPROGRAM 
!  Attfreq  WAS  NOT  CHANGED. 


!  J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

I  L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

!  Tenp  INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

I  Tenp_con  USED  TO  INSURE  VALUE  OF  TEMPERATURE  PASSED  TO  SUBPROGRAM 

!  Wat_pree  ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.  INPUT  BY 
!  OPERATOR. 

I  X  THIRD  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 

!  X*  STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

I  RESPONSE  TO  YES  OR  NO  QUESTION. 

!  Ynax  LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

! 

Freq_con-Fr«q 
T  enp_con«T enp 
Ynax-0 

Freq**DROUND<  Fraq*ZA(  2/3) ,  3  > 

FOR  J=0  TO  4  STEP  . t 
Wat  _pres-J 

CALL  AttenuationC Atom _pres,Atten,Freq,Tenp,Wat_pres) 

IF  Atten>Ynax  THEN 
Ynax-Atten 
END  IF 
NEXT  J 

Ynax -PROUND< Ynax +. 005. -Z> 

GRAPHICS  ON 
VIEWPORT  15,120,10,70 
WINDOW  0,15,0, Ynax 
AXES  .5, .005,0,0,2,2 
CLIP  OFF 
! LABEL  PLOT 
CSIZE  4,  .5 
LDIR  0 

MOVE  7.5,Ynax*1 . 1 
LORG  5 

LABEL  "ATTENUATION  VERSUS  WATER-VAPOR  PRESSURE” 

CSIZE  4, .5 

f  LABEL  HORIZONTAL  AXES 
LDIR  0 
LORG  6 

FOR  J-0  TO  15 
MOVE  J.0 
LABEL  J 
NEXT  J 

MOVE  7. 5, -Ynax*. 1 
CSIZE  4..G 

LABEL  "WATER- VAPOR  PRESSURE  mb" 

(LABEL  VERTICAL  AXES 
LORG  8 
CSIZE  4. .5 

FOR  J-0  TO  Ynax  STEP  .01 
MOVE  .25, J 


76 


LABEL  J 
NEXT  J 
LOIR  PI/Z 
L0R6  6 

MOVE  -Z,Y«ax».5 
CSIZE  4,. 6 

LABEL  '•ATTENUATION  I/m" 

! PRINT  INPUT  CONDITIONS 
LOIR  0 
L0R6  1 

MOVE  5.5, Ynax 
CSIZE  4, . G 

LABEL  USING  "K" ? "ATMOSPHERIC  PRESSURE-" , Atom_pree , "  mb" 

MOVE  5.5, Ynax* .95 
CSIZE  4, .6 

LABEL  USING  MK" ; "TEMPERATURE-" .Temp , “  CELSIUS" 

J-0 

FOR  L-1  TO  5 

FOR  X-0.  TO  G.5+J  STEP  .05 
Wat _pres-X 

CALL  Attenuat ion< Atom_pres , Atten , Freq .Temp , Wat_pres ) 
IF  INK  L/Z )-L/Z  THEN 
LINE  TYPE  Z 

ELSE 

LINE  TYPE  1 
END  IF 

PLOT  X, Atten 
NEXT  X 
LDIR  0 
LORG  Z 
CSIZE  3. .6 
MOVE  X. At ten 
LINE  TYPE  1 

LABEL  USING  "K";Freq,”Hz" 

Freq-DROUND<  Freq«Z“< - 1/3) ,3) 

J-J+1.5 
NEXT  . 

I 

PRIN'  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

GCLEAR 

FOR  J-0  TO  15 
PRINT  "  " 

NEXT  J 

PRINT  "WOULD  YOU  LIKE  TO  SEE  THE  SAME  PLOT  ONLY" 

PRINT  "WITH  ATTENUATION  VERSUS  RELATIVE  HUMIDITY?" 

LINPUT  "IF  YES  ENTER  Y,  ELSE  ENTER  N",X* 

Att_freq-FNYes<X«) 

IF  Att_freq-1  THEN 
OUTPUT  KBD» "K" ; 


WINDOW  0,100,0, Ynax 
CLIP  ON 

AXES  5, .005, 0,0, 2, 2 
CLIP  OFF 
1  LABEL  PLOT 
CSIZE  4, .5 
LOIR  0 

MOVE  50, Ynax* 1 . 1 
L0R6  5 

LABEL  "ATTENUATION  VERSUS  RELATIVE  HUMIDITY” 

CSIZE  4,.E 

'LABEL  HORIZONTAL  AXES 
LDIR  0 
LORG  6 

FOR  J-0  TO  100  STEP  10 
MOVE  J.0 
LABEL  J 
NEXT  J 

MOVE  50, -Ynax* . 1 
CSIZE  4 , . 6 

LABEL  "RELATIVE  HUMIDITY  Z" 

! LABEL  VERTICAL  AXES 
LORG  8 
CSIZE  4,.E 

FOR  J-0  TO  Ynax  STEP  .01 
MOVE  1  ,J 
LABEL  J 
NEXT  J 
LDIR  PI/2 
LORG  G 

MOVE  -15 .Ynax*. 5 
CSIZE  4  , . G 

LABEL  "ATTENUATION  1/n" 

'PRINT  INPUT  CONDITIONS 
LDIR  0 
LORG  1 

MOVE  40, Ynax 
CSIZE  4 , . G 

LABEL  USING  "K" 5 "ATMOSPHERIC  PRESSURE-" , Aton_pree . "  nb" 

MOVE  40, Ynax *.95 
CSIZE  4. .6 

LABEL  USING  "K" j "TEMPERATURE-" .Temp,"  CELSIUS" 

CLIP  ON 

T-Tenp+273 

Es- 10* ( 9 . 4-2353/T ) 

J-0 

Freq-OROUND<  Freq_con*2*< -2/3 ), 3) 

FOR  L-1  TO  5 

FOR  X-0.  TO  40+ J  STEP  .5 

Ua't_pres-£5*Aton_pres*X/< <  Aton_pres-Es*( 1 +X/ 1 00) )* 100) 


CALL  Att enuat  ion<  Atom_pree , Att en ,  Fr eq ,  T emp ,  Wat_pres ) 

IF  INK  L/Z  )“L/Z  THEN 
LINE  TYPE  2 

ELSE 

LINE  TYPE  1 
END  IF 

PLOT  X.Atten 
NEXT  X 
LOIR  0 
LORG  2 
CSI2E  3..E 
MOVE  X.Atten 
LINE  TYPE  1 

LABEL  USING  "K“ ; Freq , "Hz" 

Freq«DROUND<Freq*ZK 1/3) ,3) 

J- J+ 1 1 
NEXT  L 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT  KBD; "K" s 
ELSE 

GCLEAR 

OUTPUT  KBDj “K" s 
END  IF 

T  enp-T  enp_con 
Freq“FreqL_con 
SUBEND 

! 

! 

sub  Att_te«p< Atom _presfFreq,Temp) 

! 

! 

!  THIS  SUBPROGRAM  PLOTS  THE  ATTENUATION  OF  SOUND  VERSUS  UATER 
! VAPOR  PRESSURE  FOR  FIVE  DIFFERENT  TEMPERATURES  AT  1®  DEGREE 
(INTERVALS  AROUND  THE  INPUT  TEMPERATURE. 

! 

! 

!  INPUT  ATOMOSPHERIC  PRESSURE  IN  MILLIBARS 
!  FREQUENCY  OF  ACOUSTIC  ENERGY  IN  HERTZ 

!  TEMPERATURE  IN  CELSIUS 

( 

!  OUTPUT  PLOTS  OF  ATTENUATION  VERSUS  UATER- VAPOR  PRESSURE 
! 

(VARIABLES 

( 

!  Atom _pres  INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 

!  Atten  ATTENUATION  OF  ACOUSTIC  WAVE.  CALCULATED  IN 

!  SUBPROGRAM  ATTENUATION. 

!  Att_temp  VARIABLE  USED  TO  DETERMINE  IF  OPERATOR  WANTS 


TO  PLOT  ATTENUATION  VERSUS  WATERJ/APOR  PRESSURE 
FOR  VARIOUS  TEMPERATURE  AROUND  THE  INPUT  TEMPERATURE. 
THE  PLOT  IS  DONE  IN  THE  SUBPROGRAM  Att_temp. 

Es  SATURATION  VAPOR  PRESSURE  AT  GIVEN  TEMPERATURE. 

Freq  INPUT  FREQUENCY  OF  ECHQSOUNDER. 

I  MAIN  INDEX  FOR  LOOP  TO  CHECK  RETURNED  POWER  UNTIL 

J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

Temp  INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

Tenp_con  USED  TO  INSURE  VALUE  OF  TEMPERATURE  PASSED  TO  SUBPROGRAM 
Att_temp  UAS  NOT  CHANGED. 

Wat_pres  ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.  INPUT  BY 
OPERATOR. 

X  THIRD  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 

X*  STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

RESPONSE  TO  YES  OR  NO  QUESTION. 

Ynax  LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 


Ynax*0 

T  emp_con*  Tenp 
Tenp«Tenp+20 
FOR  J=0  TO  4  STEP  . 1 
Wat_pres“J 

CALL  Attenuat ion< Aton_pres , Atten . Fr eq , Temp . Wat_pres ) 
IF  Atten>Ynax  THEN 
Ynax -At ten 
END  IF 
NEXT  J 

Ynax«PROUND<  Ynax+ . 005 , -2 ) 

GRAPHICS  ON 
VIEWPORT  15,120.10.70 
WINDOW  0,15.0. Ynax 
AXES  .5, .005, 0,0, 2, 2 
CLIP  OFF 
I  LABEL  PLOT 
CSIZE  4, .6 
LOIR  0 

MOVE  7. 5, Ynax* 1.1 
LORG  5 

LABEL  "ATTENUATION  VERSUS  WATER-VAPOR  PRESSURE" 

CSIZE  4, .5 

1  LABEL  HORIZONTAL  AXES 
LDIR  0 
LORG  G 

FOR  J-0  TO  IS 
MOVE  J.0 
LABEL  J 
NEXT  J 

MOVE  7. 5,- Ynax*. 1 


'FREQUENCY-", Freq,"  HZ" 


CSIZE  4 . . E 
LABEL  " WATER- VAPOR  PRESSURE  nb" 

(LABEL  VERTICAL  AXES 
LORG  8 
CSIZE  4,.B 

FOR  J-0  TO  Ymax  STEP  .01 
MOVE  .ZS,J 
LABEL  J 
NEXT  J 
LDIR  PI/Z 
LORE  6 

MOVE  -Z,Ymax«.5 
CSIZE  4, .6 

LABEL  "ATTENUATION  1/m" 

(PRINT  INPUT  CONDITIONS 

LDIR  0 

LORG  1 

MOVE  6, Ymax 

CSIZE  4..G 

LABEL  USING  "K" ; "ATMOSPHERIC  PRESSURE-" . Atom_pres . "  mb" 

MOVE  6, Ymax*. 95 
CSIZE  4. .6 

LABEL  USING  "K" i "FREQUENCY-" .Freq,"  HZ" 

CLIP  ON 
J-0 

FOR  L-1  TO  S 

FOR  X-0.  TO  G+J  STEP  .05 
Uat_pres»X 

IF  INT( L/Z  >*L/Z  THEN 
LINE  TYPE  3 

ELSE 

LINE  TYPE  1 
END  IF 

CALL  Attenuation<  Atom _pres, Atten, Freq, Temp, Uat_pres) 
PLOT  X, Atten 
NEXT  X 
LDIR  0 
LORG  Z 
CSIZE  3. .6 
MOVE  X, Atten 
LINE  TYPE  1 

LABEL  USING  "K"sTemp,"C" 

Temp- Temp- 10 
J-J+l .6 
NEXT  L 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT  KBO; "K" ; 

FOR  J-0  TO  15 


PRINT  "  “ 

NEXT  J 

! 

PRINT  "WOULD  YOU  ARE  LIKE  TO  SEE  THE  SAME  PLOT  ONLY" 
PRINT  "WITH  ATTENUATION  VERSUS  RELATIVE  HUMIDITY" 
LINPUT  "IF  SO  ENTER  Y,  IF  NOT  ENTER  N",X* 
At4_temp-FNYes<  X$ ) 

IF  Att_temp-  1  THEN 
OUTPUT  KBDi "K" ; 

T  ernp- T enp_con 
WINDOW  0,100,0, Ymax 
AXES  1C, .005,0,0,2,2 
CLIP  OFF 
! LABEL  PLOT 
CSIZE  4, .6 
LOIR  0 

MOVE  50, Ymax* t . 1 
LORG  5 

LABEL  "ATTENUATION  VERSUS  RELATIVE  HUMIDITY" 

! LABEL  HORIZONTAL  AXIS 
LDIR  0 
LORG  G 

FOR  J-0  TO  100  STEP  10 
MOVE  J,0 
LABEL  J 
NEXT  J 

MOVE  50, -Ymax* . 1 
CSIZE  4 , . G 

LABEL  "RELATIVE  HUMIDITY" 

(LABEL  VERTICAL  AXES 
LORG  8 
CSIZE  4,. 6 

FOR  J-0  TO  Ymax  STEP  .01 
MOVE  Z,J 
LABEL  J 
NEXT  J 
LDIR  PI/Z 
LORG  G 

MOVE  -  10, Ymax* . 5 
CSIZE  4, .6 

LABEL  "ATTENUATION  1/m" 

(PRINT  INPUT  CONDITIONS 

LDIR  0 

LORG  6 

MOVE  50, Ymax 

CSIZE  4,.G 

LABEL  USING  "K" j "ATMOSPHERIC  PRESSURE-” ,Atom_pres , 
MOVE  50, Ymax *.9 5 
CSIZE  4,.G 

LABEL  USING  "K" ; "FREQUENCY-" ,Freq,"  Hz" 


CLIP  ON 
J-0 

Temp“Te«p-20 
FOR  L-1  TO  B 
T“To«p+273 
Es* 1 0A ( 9 . 4-2  353/ T ) 

FOR  X-0  TO  90- J  STEP  .5 

Uat_pree*Ea*Ato«_pre8*X/(  ( Atoro_pree-Ea*( 1 +X/ 100) )* 100) 
IF  INT*L/2)«L/2  THEN 
LINE  TYPE  3 

ELSE 

LINE  TYPE  1 
END  IF 

CALL  Attenuation*  Aton_prea , Atten,Freq,Temp,Uat_prea ) 
PLOT  X, Atten 
NEXT  X 
LDIR  0 
LORG  2 
CSIZE  3..S 
MOVE  X, At ten 
LINE  TYPE  1 

LABEL  USING  MK"sTerlp,‘•  C" 

Tenp-Tenp+10 
J-J+1S 
NEXT  L 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT  KBDj "K" i 
ELSE 

OUTPUT  KBDs "X" ; 

END  IF 

T  enp«T  ertp_con 
SUBEND 


aub  Log_log*  Xnin.Xnax , Ynin, Ynax ,Xvar( • ) ,J .Range*  * ) ,L , Title*, XI abelS, Ylabel 


THIS  SUBROUTINE  MAKES  A  LOG-LOG  PLOT  OF  DATA  PASSED  FROM  THE 
MAIN  PROGRAM. 

INPUT  MINIMUM  VALUES  OF  X  AND  Y  FOR  PLOT 
MAXIMUM  VALUES  OF  X  AND  Y  FOR  PLOT 
X  AND  Y  VALUES  TO  BE  PLOTTED 
TITLE  OF  PLOT 
LABELS  FOR  X  AND  Y  AXIS 
OUTPUT  LOG- LOG  PLOT 
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I  VARIABLES 

!  Dx  STEP  SIZE  FOR  X  AXIS  FOR  VARIOUS  PLOTS. 

!  Dy  STEP  SIZE  FOR  Y  AXIS  FOR  VARIOUS  PLOTS. 

!  J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

I  L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

!  Titles  STRING  PASSED  TO  SUBPROGRAM  PT  FOR  TITLE  OF  PLOT. 

!  XlabelS  LABEL  ON  X  AXIS  PASSED  TO  SUBPROGRAM  Pt  FOR  PLOTTING. 

I  Xnax  VALUE  OF  LARGEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

I  Xnin  VALUE  OF  SMALLEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

!  Xrange  VALUE  OF  RANGE  OF  X  VALUES  FOR  VARIOUS  PLOTS. 

!  Xvar<*>  VALUE  OF  ARRAY  OF  X  VALUES  TO  BE  PLOTTED  IN  SUBPROGRAM  Pt . 

!  YlabelS  LABEL  ON  Y  AXIS  PASSED  TO  SUBPROGRAM  Pt  FOR  PLOTTING. 

!  Ymax  LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

!  Ynin  SMALLEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

!  Yrange  RANGE  OF  Y  VARIABLES  FOR  VARIOUS  PLOTS. 

) 

Ox-1 

Dy-t 

Xrange-ABS<  Xnax- Xnin) 

Yrange«ABS<  Ynax-Ynin) 

UINDOU  Xnin, Xnax , Ynin, Ymax 
AXES  Dx ,Dy, Xnin, Ynin, 1 , 1 
CLIP  OFF 
I 

! LABEL  PLOT 
CSIZE  4..G 
LDIR  © 

LORG  5 

MOVE  Xnin+ .5* Xrange , 1 . 1*Yrange 
TitleS-TRIMS(TitleS) 

LABEL  Titles 
! LABEL  HORIZONTAL  AXES 
LDIR  0 
LORG  5 

FOR  J-Xnin  TO  Xnax  STEP  Dx 
CSIZE  4, .6 

MOVE  J- .01 3«Xrange ,- .0G7*Yrange 
LABEL  "10" 

MOVE  J+ .01 3*Xrange,- .033» Yrange 
CSIZE  Z 
LABEL  J 
NEXT  J 

MOVE  Xnin+. 5* Xrange,-. 1Z*Yrange 
LORG  5 
CSIZE  4 , . G 

XlabelS- TRIMS< XlabelS) 

LABEI  XlabelS 
! LABEL  VERTICAL  AXES 
LORG  8 

FOR  J-0  TO  Ynax  STEP  Dy 


CSIZE  4..G 

MOVE  Xnin- ,©Z5*XranQef J 
LABEL  "I©" 

CSIZE  Z 

MOVE  Xmin- .0©Z5«Xrange, J+.©3*Yrange 
LABEL  J 
NEXT  J 
LOIR  PI/Z 
LORG  G 

MOVE  X«in-. 1*Xrange, .5»Yranoe 
CSIZE  4, .6 

Ylabel*«TRIM*<  Ylabel*) 

LABEL  Ylabel « 

CLIP  ON 
FOR  J-1  TO  L 

PLOT  LGT<  Xvar<  J) ) ,LGT<  Range(  J) ) 

NEXT  J 
SUBENO 


APPENDIX  B 
COMPUTER  PROGRAM  2 


FULLER,  ROBERT  PROG_Z  FREQUENCY 

9  SEP  85 


*******PURPOSE«**»*** 

THIS  PROGRAM  WILL  TAKE  PARAMETERS  OF  AN  ECHOSOUNDER  AND 
ESTIMATE  THE  RANGE  AS  A  FUNCTION  OF  FREQUENCY.  THE  FOLLOWING  INPUTS 
ARE  REQUIRED: 

ATMOSPHERIC  DATA 

1 ) ATOMOSPHERIC  PRESSURE  in  millibars 

Z ) SELECT  ONE  OF  FOUR  ATMOSPHERIC  PROFILES  FOR  THE  TEMPERATURE 
STRUCTURE  PROFILE. 

a)FOR  ONE  PROFILE  THE  HEIGHT  OF  THE  INVERSION  LAYER  IS  ALSO  INPUT 

3)  TEMPERATURE  IN  DEGREES  CELSIUS 

4)  WATER  VAPOR  PRESSURE  IN  millibars 

ECHOSOUNDER  DATA 

5)  ANTENNA  DIAMETER  IN  METERS 

G) POWER  TRANSMITTED  BY  ECHOSOUNDER  IN  WATTS 
7)PULSE  LENGTH  OF  THE  TRANSMITTED  ACOUSTIC  ENERGY 

THE  PROGRAM  OUTPUTS  THE  FOLLOWING  GRAPHS  TO  AN  EXTERNAL  PLOTTER: 
t )RANGE( m. )  VERSUS  FREQUENCY  OF  ECHOSOUNDER 

2 ) RANGE( m. )  VERSUS  EXCESS  ATTENUATION  FOR  VARIOUS  FREQUENCUIES 


VARIABLES 


!  Again  TEST  VALUE  IN  JHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

!  WISHES  TO  MAKE  ANOUTHER  RUN. 

!  Ano_change  TEST  VALUE  IN  WHILE  LOOP  TO  DETERMINE  IF  OPERATOR 
!  WISHES  TO  CHANGE  THE  VALUE  OF  INPUT  DATA  BEFORE  MAKING 

!  ANOUTHER  RUN. 

!  Ant_area  ECHOSOUNDER  ANTENNA  AREA  CALCULATED  FROM  INPUT  OF 
!  ANTENNA  DIAMETER. 

!  Ani.diam  INPUT  OF  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS. 

!  Atom_pres  INPUT  Or  ATMOSPHERIC  PRESSURE  IN  mb. 

I  Atten  ATTENUATION  OF  ACOUSTIC  WAVE.  CALCULATED  IN 

!  SUBPROGRAM  ATTENUATION. 

!  Bn  BANDWIDTH  OF  RECEIVER  NEEDED  TO  RECEIVE  SIGNALS 

!  DOPPLER  SHIFFED  BY  THREE  METER  PER  SECOND 

!  VERTICAL  VELOCITIES. 

!  C  VELOCITY  OF  SOUND  CALCULATION  AT  ZERO  CELSIUS. 

!  Cl  VARIABLE  USED  TO  SCALE  INVERSION  HEIGHT  IN 

I  TEMPERATURE  STRUCTURE  PROFILE  THREE. 

!  C3  VARIABLE  USED  IN  CALCULATING  THE  FOURTH  CteZ  PROFILE. 

I  CneZ(  • )  ARRAY  OF  VALUES  OF  THE  ACOUSTIC  REFRACTIVE  INDEX 


Dx 

Dy 

Eps i 1 on 

Et 

Er 

Exc_att( * 
F 


IND 

Interval 

Inver 

J 

K 

L 

Meso_up 


Neu_va 

Noise 
Pou_back  ( * 
Pou_ret(  * ) 
Pow_trans 
Profile 
Pstar 

Pulse 

R 

Rt 

Range< • ) 
Remainder 


PARAMETER.  CALCULATED  BASED  ON  SELECTION  OF 
PROFILE  FOR  TEMPERATURE  STRUCTURE  PARAMETER 
AND  ASSUMED  PROFILE  FOR  THE  VELOCITY  STRUCTURE 
PARAMETER. 

VALUE  OF  THE  TEMPERATURE  STRUCTURE  PARAMETER. 

THE  OPERATOR  SELECTS  A  PROFILE  FROM  SEVERAL  PROVIDED. 

VALUE  OF  THE  VELOCITY  STRUCTURE  PARAMETER. 

VALUES  BASED  ON  CALCULATION  USING  ASSUMED  DISSAPATION 
RATE. 

STEP  SIZE  FOR  X  AXIS  FOR  VARIOUS  PLOTS. 

STEP  SIZE  FOR  Y  AXIS  FOR  VARIOUS  PLOTS. 

DISSAPATION  RATE  USED  IN  CALCULATION  OF  CveZ. 

TRANSMISSION  EFFICIENCY  OF  ECHOSOUNDER. 

EFFICIENCY  OF  ECHOSOUNDER  WHEN  RECEIVING  ACOUSTIC 
BACKSCATTER. 

*)  EXCESS  “ ATTENUATION"  AT  GIVEN  RANGE. 

VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.  IS  THE 
RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 
FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 
ATTENUATION. 

INPUT  FREQUENCY  OF  ECHOSOUNDER. 

ANTENNA  EFFECTIVE  APERATURE  FACTOR. 

VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

MAIN  INDEX  FOR  LOOP  TO  CHECK  RETURNED  POUER  UNTIL 
POUER  UAS  LESS  THAN  MINIMUM  DETECTABLE  SIGNAL. 

INDEX  USED  FOR  LOOP  FOR  DIFFERENT  FREQUENCIES 
MAXIMUM  EFFECTIVE-SCATTERING  VOLUME  THICKNESS. 

HEIGHT  OF  INVERSION  LAYER. 

FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

WAVENUMBER 

SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

TEST  VALUE  TO  CHECK  RESPONSES  WHEN  OPERATOR  IS  INPUTTING 
RESPONSES. 

TEST  VALUE  IN  CALCULATION  OF  COHERENCE  LENTH  IN  CALCULATION 
OF  EXCESS  ATTENUATION. 

TEST  VALUE  IN  CHECKING  IF  OPERATOR  WISHES  TO  CHANGE 
A  VARIABLE  BEFORE  A  NEW  RUN. 

ASSUMED  MINIMUM  DETECTABLE  SIGNAL. 

)  POUER  BACKSCATTERED  FROM  GIVEN  RANGE. 

POUER  BACKSCATTERED  TO  ECHOSOUNDER  FROM  GIVEN  RANGE. 

INPUT  OF  POUER  SUPPLIED  TO  TRANSDUCER  OF  ECHOSOUNDER. 
OPERATOR  INPUT  OF  CteZ  PROFILE  FROM  AVAILABLE  PROFILES. 
VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 
CALCULATION. 

TRANSMITTED  PULSE  LENGTH  IN  MILLISECONDS.  INPUT  BY 
OPERATOR, 

RANGE  FOR  INPUT  TO  INTEGRAL  OF  EXCESS  ATTENUATION. 

VARIABLE  USED  IN  SCALING  INVERSION  HEIGHT  FOR  PROFILE  THREE. 
ARRAY  OF  RANGE  VALUES. 

REMAINDER  OF  MODULO  FUNCTION  USED  TO  DECREES  THE 
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NUMBER  OF  PASSES  THROUGH  THE  INTEGRATION  FOR  EXCESS 
ATTENUATION. 

Rge  VALUE  OF  RANGE  IN  INTEGRAL  FOR  EXCESS  ATTENUATION. 

Rho  CORRELATION  LENGTH  USED  IN  CALCULATION  OF  EXCESS 

ATTENUATION. 

SiQma  FRACTION  OF  ENERGY  BACKSCATTERED  FROM  GIVEN  RAN6E. 

Speed_aound  SPEED  OF  SOUND  IN  AIR  AT  INPUT  TEMPERATURE. 

Sumpou_back  SUM  OF  BACKSCATTERED  ENERGY 

T  INPUT  TEMPERATURE  IN  DEGREES  KELVIN. 

Temp  INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

Tempt  VARIABLE  STRING  USED  IN  FUNCTION  YES. 

Title*  STRING  PASSED  TO  SUBPROGRAM  Semi_log  FOR  TITLE  OF  PLOT. 

Tatar  INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

IN  SUBPROGRAM  ATTENUATION. 

Var  USED  TO  SELECT  THE  VARIABLE  THE  OPERATOR  WISHED  TO 

CHANGE  BEFORE  MAKING  ANOTHER  RUN. 

Wat_prea  ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.  INPUT  BY 
OPERATOR. 

X  THIRO  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 

X*  STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

RESPONSE  TO  YES  OR  NO  QUESTION. 

Xlabel*  LABEL  ON  X  AXIS  PASSED  TO  SUBPROGRAM  Semi_log  FOR  PLOTTING. 

Xmax  VALUE  OF  LARGEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

Xmin  VALUE  OF  SMALLEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

Xran0e  VALUE  OF  RANGE  OF  X  VALUES  FOR  VARIOUS  PLOTS. 

Xvar( * )  VALUE  OF  ARRAY  OF  X  VALUES  TO  BE  PLOTTED  IN  SUBPROGRAM  Semi 

Ylabel*  LABEL  ON  Y  AXIS  PASSED  TO  SUBPROGRAM  Semi_loQ  FOR  PLOTTING. 

Ymax  LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

Ymin  SMALLEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

Yrange  RANGE  OF  Y  VARIABLES  FOR  VARIOUS  PLOTS. 

Ze  INTERMEDIATE  VALUE  IN  CALCULATION  OF  EXCESS  ATTENUATION. 


(DECLARE  VARIABLES 
INTEGER  I 

DI M  Pou_back ( 1 500  > , Pow_ret  <  1 500 ) , Ranoe(  1 500  > 

DIM  CneZ< 1500) ,Exc_att( 1500,30) ,Xvar< 1 500 ) , Ran< 30 ) 

DIM  T itletf  50]  ,Xlabel*l50]  ,Ylabel*t 1EJ 
! 

PLOTTER  IS  705 , "HPGL" 

LINE  TYPE  1 

! 

(INPUT  ATMOSPHERIC  DATA 

INPUT  "ENTER  TEMPERATURE  IN  DEGREES  CELCIUS" .Temp 
INPUT  "ENTER  ATMOSPHERIC  PRESSURE  IN  MILLIBARS" ,Atom_prea 
INPUT  "ENTER  WATER  VAPOR  PRESSURE  IN  MILLIBARS" , Wat.prea 


! INPUT  ECHOSOUNDER  DATA 

INPUT  "ENTER  ANTENNA  DIAM  IN  METERS” , Arvt_diam 
INPUT  "ENTER  ECHOSOUNDER  PULSE  LENGTH  IN  MILLISECONDS'" .Pulse 
INPUT  "ENTER  POWER  TO  TRANSMITTER  IN  WATTS", Pou.trana 
Et-.Z5  (TRANSMIT  EFFICIENCY 

EP-.25  (RECEIVER  EFFICIENCY 

G-.40  (ANTENNA  EFFECTIVE  APERATURE  FACTOR 

I 
( 

(SELECT  THE  PROFILE  FOR  THE  TEMPERATURE  STRUCTURE 
PRINT  "YOU  HAVE  A  CHOICE  OF  TEMPERATURE  STRUCTURE  PROFILES" 

PRINT  "PROFILE  TO  BE  USED." 

PRINT  "YOUR  SELECTIONS  ARE" 

PRINT  "  1 

PRINT  " 

PRINT  " 

PRINT  "  " 

PRINT  “  Z 

PRINT  " 

PRINT  " 

PRINT  "  “ 

PRINT  "  3 

PRINT  ” 

PRINT  " 

PRINT  M 
PRINT  " 

PRINT  "  “ 

PRINT  ”  4 

PRINT  " 

! 

PRINT  "  " 

PRINT  "  " 

Meaa_up« ( 

WHILE  Meaa_up-1 

INPUT  "ENTER  THE  DESIRED  PROFILE  <1  OR  Z  OR  3  OR  4>", Profile 
IF  Profile-4  THEN 

INPUT  "ENTER  THE  HEIGHT  OF  THE  INVERSION  IN  METERS" . Inver 
END  IF 

IF  Profile-1  THEN 
Meas_up-0 

ELSE 

IF  Profile-2  THEN 
Meaa_up-0 

ELSE 

IF  Profile-3  THEN 
Meas_up-0 

ELSE 

IF  Profile-4  THEN 
Meaa_up-0 


A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS  " 
PRESENTED  IN  WALTERS/KUNDEL  081  PAGE  398  WHICH  " 
GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 

THE  TEMPERATURE  STRUCTURE  PROFILE  ABOVE  BUT" 

WITH  A  FACTOR  OF  TWO  TO  APPROXIMATE  LOOKING- 
UP  A  THERMAL  PLOOM" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 
PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398  FOR" 
NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH  HEIGHT  " 

IS  EXP< - .001 ‘HEIGHT  ABOVE  G5  METERS  AND  HEIGHT" 

TO  THE  -1.46  BELOW  65  METERS" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 

HEIGHT  TO  THE  -4/3  " 


ELSE 


PRINT  Prof ile,"  UAS  NOT  ONE  OF  THE  OPTIONS! 
Mess_up*1 
END  IF 
END  IF 
END  IF 
END  IF 
END  UHILE 
OUTPUT  KBDs,,K,‘s 


Again* 1 
UHILE  Again*! 

FOR  J-l  TO  I 

Pou_back(  J)*0 
Pou_ret(  J)*0 
Range  ( J )*© 

CneZ<  J )*0 
FOR  L-l  TO  Ind 

Exc_att( J,L>*0 
NEXT  L 
NEXT  J 
I 

! CONVERT  TEMPERATURE  TO  KELVIN 
T*Tenp+Z73 
! 

(CALCULATE  SPEED  OF  SOUND  BASED  ON  INPUT  TEMPERATURE1 CELCIUS) 
Speed_sound*Z0. 05*  <  T ) * . 5 

(CALCULATE  THE  SPEED  OF  SOUND  AT  0  OEGREES  CELCIUS 
C=Z0.05*Z73'> .  5 
! 

1-1  (INDEX  FOR  LOOP 
Ind»1 

Pou_back<0>-0  (INITIALIZE  VARIABLE  FOR  POUER  BACKSCATERRED 
Ant_area-PI*(Ani_dian/Z)"Z  (ANTENNA  AREA 
Intervals  Speed_oound*Pulse* 1 . E-3)/Z 
I 

(INITIALIZE  PLOTTING  PARAMETERS  SO  CAN  COMPARE  DURING  RUN  TO  SCALE 

(PLOTS 

Ymin-1000 

Ynax*0 

(CALCULATE  RETURNED  SIGNAL  POUER  UNTIL  IT  IS  LESS  THAN  BACKGROUND 
(FOR  EACH  FREQUENCY. 

Range( 0)-0 

FOR  Freq-100  TO  3000  STEP  100 

Bn*Z*Freq*(  1  - 1 /<  6/Speed_sound+ 1 ) > ! BANDUIDTH  FOR  3M/S  VERTICAL  VELOCITY 
Noiee- 1 . 38E-Z3*Bn*<  T )+Z . E- 14 1  MINIMUM  DETECTABLE  SIGNAL  AS 
!  JOHNSON  NOISE( NEGLIGABLE )+  ESTIMATED  BACKGROUND 
(CALCULATE  THE  ATTENUATION  COEFFICIENT  BASED  ON  ATMOSPHERIC  DATA 
CALL  Attenuat ion< Atom_pres , At ten ,Freq , Temp , Uat_pres ) 
K-Z*PI*Freq/Speed_sound  ! UAVENUMBER 
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Sunpou_back “0 
REPEAT 

Range ( I )«Range( 1-1 )+2 
SELECT  Profile 
CASE  1 

•THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
IAS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 

IUAS  FOR  MIDDAY  AND  CLEAR  WEATHER  ABOVE  THE 
! TULAROSA  BASIN,  NEW  MEXICO. 

! AN  ADDITONAL  FACTOR  OF  1.8  WAS  INCLUDED  AS  PER  DIFFERENT 
•AVERAGING  TIME. 

CteZ-Z . 1 Z*Range( I ) *  < - 1 . 1 B ) 

CASE  Z 

I  THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
IAS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 

(WAS  FOR  MIDDAY  AND  CLEAR  WEATHER  ABOVE  THE 
! TULAROSA  BASIN,  NEW  MEXICO. 

! AN  ADDITONAL  FACTOR  OF  1 . 8  WAS  INCLUDED  AS  PER  DIFFERENT 
•AVERAGING  TIME  AND  AN  ADDITIONAL  FACTOR  OF  TWO  WAS 
(INCLUDED  TO  APPROXIMATE  LOOKING  UP  A  THERMAL  PLOOM 
Cte2-2*Z.  12*Ranged)A<-l.  16) 

CASE  3 

(THIS  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 

(PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 

(NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH  HEIGHT  " 

(IS  EXP(-. 001 ‘HEIGHT  ABOVE  B5  METERS  AND  HEIGHT" 

(TO  THE  - 1 .46  BELOW  65  METERS’* 

IF  Ranged) <65  THEN 

CteZ-75. 5*Range< I ) *  < - 1 . 46 ) 

ELSE 

CteZ-3.66E-Z*EXP< - .001 *Range< I ) ) 

END  IF 
CASE  4 

(CALCULATE  THE  TEMPERATURE  STRUCTURE  FACTOR  BASED  ON  ASSUMPTION 
!  THAT  IT’S  PROPORTIONAL  TO  HEIGHT  TO  THE  -4/3  POWER 

•  THAT  THE  SURFACE  VERTICAL  HEAT  FLUX  IS  .095  C.m./sec 

(THE  DEPENDENCE  WITH  HEIGHT  ABOVE  THE  INVERSION  I  PULLED  FROM 
(MY  LEFT  EAR 

(EQUATION  FROM  NEFF,  1975 
SELECT  Range(I)/ Inver 
CASE  <.9 

C3-<< . 0Z  4 ) * ( T ) “ ( .667)) 

CteZ-C3*(  Ranged  )>'*(- 1.33) 

CASE  .9  TO  1 
CteZ-CteZ 
Cl-CieZ 
R1 »Range( I ) 

CASE  1  TO  1.3 

CteZ“ 10*  <  <  LGT<  C3*Range< I ) )-LGT<  C3*R1 ))+LGT<C1 )) 


RZ-Range( I ) 

C2-CteZ 
CASE  ELSE 

Cte2”C3»( Range ( I )*(  - 1 . 33 )-RZ“< -  1 . 33 ) )+CZ 
END  SELECT 
END  SELECT 

! CALCULATE  VELOCITY  STRUCTURE  FACTOR.  FORMULA  FROM  GAYNOR  77 
Epsilon-< .Z 858/ Range ( I ) >*( 1+.07*< Range(  I ) )A( .5) >“< 1.5) 
CveZ-Z*Epsilon“< .557) 

I  CALCULATE  ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  FACTOR. 

! FORMULA  FROM  TETARCKI 
CneZ<I)-<CteZ/(Z.9BE-i5)>  +  <CveZ/(C*C>> 

(CALCULATE  THE  FRACTION  OF  POUER  BACKSCATTERED  FROM  INTERVAL. 
(EQUATION  FROM  NEFF,  1 S75 
Signa-< .0039*<K‘< 1/3) )*Cte2 )/< T>*2 
l 

(CALCULATE  THE  EXCESS  ATTENUATION.  EQUATION  FROM  CLIFFORD.  1980 
(THE  EXCESS  ATTENUATION  IS  Exc_att< I , L . X ) 

IF  I <50  THEN 

Reinainder-0 

ELSE 

IF  I <150  THEN 

Remainder-Range< I )  MODULO  10 

ELSE 

Remainder-Range ( I )  MODULO  Z0 
END  IF 
END  IF 

IF  Remainder3©  THEN 
Rho-0 
L-0 
H-0 

Rge-Range(I)  (CONSTANT  IN  INTEGRAL 
R-0 

FOR  J-0  TO  Z*I 

F”CneZ< INT< J/2+1 ) > 

F”F*< 1 -R/( Rge  )'< 1 . 57>  +  <  R/<  Rge ) > A( 1 . 57) ) 

IF  J>0  THEN 

IF  J<Z*I  THEN 

IF  INT<  J/Z )“ J/Z  THEN 
L-L+F 
F-0 

ELSE 

H-H+F 
F-0 
END  IF 
END  IF 
END  IF 
Rho-Rho+F 
R-R4  1 


NEXT  J 


Rho-Rho+4*L+Z«H 

Rho-<  <  <  Rho« . 33)*K*K»1.46)A<-.S)> 

N-<  Ant_diam/Rho>AZ 
IF  N<-1  THEN 

Ze-1/< 1+N) 

ELSE 

Ze-1 .5/< 1+N) 

'STEP  OF  1.5— >SEE  CLIFFORD  1980 
END  IF 

Exc_att< I , Ind)-Ze«Ze 

ELSE 

Exc_att<  I , Ind)-Exc_att(  1-1  ,Ind) 

END  IF 
! 

! CALCULATE  THE  POWER  BACKSCATTERED 

Pou_back(  I )-<  Pou_trans*Et-Sumpou_back )*EXP< -Atten*Range( I ) ) 

Pou_back< I )»Pou_back( I )*Interval*Exc_att(  I , Ind)*Signa 
Sumpow__back-Sunpou_back+Pou_back(  I ) 

(CALCULATE  THE  POuIr  BACKSCATTERED  TO  THE  ANTENNA 

Pou_ret( I )-Pou_back( I >*EXP< -At ten* Rangel  I ) )«Ant_area*G*Er/Range< I )  AZ 
PRINT  "RANGE-" .Range(I) 

PRINT  "POUER  RETURNED-"  ,Pou_retd) 

I-I+1 

UNTIL  Pou_ret< I- 1 ) >0  AND  Pou_ret(  I- 1 X-Noise 
Ran( Ind)-Range< I-Z ) 

IF  Range< I-Z ) >Ynax  THEN 
Ynax-Range< I-Z ) 

END  IF 

IF  Ranged -ZKYnin  THEN 
Ymin»Range( I-Z ) 

END  IF 
Ind-Ind+1 
Again- 1 
FOR  J-1  TO  I 

Pou_back( J)-0 
Pou_ret<  J )-0 
Rangel J)-0 
CneZ<  J )-0 
NEXT  J 

PRINT  "FREQUENCY-", Freq 
NEXT  Freq 
! 

I 

OUTPUT  KBD; "K" s 

PRINT  "INPUT  CONDITIONS" 

! 

PRINT  "  " 

PRINT  USING  "K"; "TEMPERATURE-  " . Tenp , "CELCIUS" 

PRINT  USING  "K"s "ATMOSPHERIC  PRESSURE-  " , Aton.prea , "mb" 

PRINT  USING  "K" ; "UATER  VAPOR  PRESSURE-  " , Ulat_pres . "mb" 


PRINT  USING  "K"  5  "PULSE  LENGTH-  "  .Pulse, “ns" 

PRINT  USING  "K“ ; "ANTENNA  DIAMETER-  " ,Ant_dia«, "  in." 

PRINT  USING  "K"; "POWER  TRANSMITTED-  “ ,Pou_trans , “  UATTS" 

PRINT  "TEMPERATURE  STRUCTURE  PROFILE  USED  ", Profile 
IF  Profile-4  THEN 

PRINT  "INVERSION  HEIGHT  ",  Inver 
END  IF 

SELECT  Profile 
CASE  1 

PRINT  M  1  A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA 

PRINT  “  AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398 

PRINT  "  WHICH  GIVES  A  HEIGHT  TO  THE  - 1 . 1 G  PROFILE" 

CASE  Z 

PRINT  “  Z  A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA 

PRINT  "  AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PA6E  398 

PRINT  "  WHICH  GIVES  A  HEIGHT  TO  THE  -1.1G  PROFILE- 

PRINT  "  BUT  WITH  A  FACTOR  OF  TWO  TO  APPROXIMATE" 

PRINT  "  LOOKING  UP  A  THERMAL  PLOOM" 

CASE  3 

PRINT  3  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA 

PRINT  "  AS  PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398" 

PRINT  "  FOR  NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH" 

PRINT  "  HEIGHT  IS  EXP< - .001 ‘HEIGHT  ABOVE  G5  METERS" 

PRINT  "  AND  HEIGHT  TO  THE  -1.4G  BELOW  G5  METERS" 

CASE  4 

PRINT  "  4  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 

PRINT  "  HEIGHT  TO  THE  -4/3  " 

END  SELECT 
! 

PRINT  USING  "K’VTHE  MINIMUM  DETECTABLE  SIGNAL  WAS  SET  AT"  .Noise, "UATTS 
I 

PRINT  "  " 

I 

! 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

OUTPUT  KBDi "K" ; 

I 

GRAPHICS  ON 
VIEWPORT  15. 1Z0. 10,70 
I 
I 

I  PLOT  RANGE  VERSUS  FREQUENCY 

Xnin-0 

Xmax-3000 

Ymin-PROUND< Ymin- 50 , Z ) 

Ywax-PROUND<  Y«ax+50,Z  > 

Ox -500 
Dy-100 

Xrange-Xnax-Xnin 
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Yrange-Ymax-Ynin 
UINDOU  Xnin.Xnax , Ynin, Ymax 
AXES  Ox ,Qy .Xnin.Ynin.Z , 1 
CLIP  OFF 
! 

! LABEL  PLOT 
CSIZE  4 , . G 
LDIR  0 
L0R6  5 

MOVE  . 5*Xrange,Ynax*1 .05 
LABEL  "RANGE  VERSUS  FREQUENCY" 

! LABEL  HORIZONTAL  AXES 
LOIR  0 
LORG  5 

FOR  J-0  TO  Xmax  STEP  Ox 
C SIZE  4,.B 

MOVE  J,Ymin-.05*Yrange 
LABEL  J 
NEXT  J 

MOVE  . 5*Xrange,Ynin- . 1 »Yrange 

CSIZE  4..B 

LABEL  "FREQUENCY" 

(LABEL  VERTICAL  AXES 
LORG  8 

FOR  J-Ymin  TO  Ymax  STEP  Dy 
CSIZE  4, . G 

MOVE  Xmin- . 00Z5*Xrange , J 
LABEL  USING  "K"»J 
NEXT  J 
LDIR  PI/Z 
LORG  G 

MOVE  Xnin- . 1 5*Xrange , . 5* Yrange+Ymin 
CSIZE  4 , . G 

LABEL  "RANGE  (METERS)" 

CLIP  ON 
Ind»1 

FOR  Freq-100  TO  3000  STEP  100 
PLOT  Freq,Ran( Ind) 

Ind«Ind+ 1 
NEXT  Freq 

! 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT  KBO: "K" j 


PLOT  RANGE  VERSUS  EXCESS  ATTENUATION 
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UNCLASSIFIED 


F/G  28/1 
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Y«ax“4 

Xnax-1 

Y*in-0 

Range<0>«0 

Title*-"RAN6E  vs  EXCESS  ATTENUATION" 
Xlabel*«"EXCESS  ATTENUATION” 
Ylabel*-”RAN6E  < METERS >“ 

Dx-.l 

Dy*1 

Xrange”ABS( Xmax-Xnin ) 

Yranoe-ABS( Y«ax-Y«in) 

UINOOU  Xnin.Xaax .Ynin.Ymax 
AXES  Ox .Dy.Xain, Yain, t , I 
CLIP  OFF 
I 

(LABEL  PLOT 
CSI2E  4..B 
LOIR  0 
LORG  5 

MOVE  X«in+ . 5*Xrange ,1.1 • Yrange 
Title»-TRIM*< Title*) 

LABEL  Title* 

(LABEL  HORIZONTAL  AXES 
LOIR  0 
LORG  5 

FOR  J-X«in  TO  Xnax  STEP  Ox 
MOVE  J,- .033* Yrange 
CSIZE  4, .6 
LABEL  J 
NEXT  J 

MOVE  X«in+.5*Xrange,-. 12*Yrange 
LORG  5 
CSIZE  4..G 

XI abel *“TRIM»< XI abel *  > 

LABEL  XI abel* 

(LABEL  VERTICAL  AXES 
LORG  8 

FOR  J-0  TO  Y*ax  STEP  Dy 
CSIZE  4..G 

MOVE  Xnin-.0Z5*Xrange, J 
LABEL  "10" 

CSIZE  2 

MOVE  Xnin- .00ZS*Xrange ,J+. 03* Yrange 
LABEL  J 
NEXT  J 
LOIR  PI/2 
L0R6  6 

MOVE  Xnin- . ( *Xrange, .5* Yrange 
CSIZE  4,.G 

Ylabel*-TRIMS<  Ylabel*) 


LABEL  Ylabel* 

FOR  L-5  TO  Ind-1  STEP  5 
X-1 

WHILE  Exc_att(X.L)>0 

Xvar<X)-Exc_att<X,L> 

Range  ( X)-Range( X-1 >+Z 
X-X+1 
END  WHILE 

IF  INT(L/Z)-L/Z  THEN 
LINE  TYPE  1 

ELSE 

LINE  TYPE  3 
END  IF 
Freq-L« 100 

CALL  Se«i__log(  Freq,Xvar<  *  )  ,X-1 , Rangel  * )  ,X- 1 ) 

NEXT  L 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

! 

6CLEAR 

OUTPUT  KBOs“K"i 

LINPUT  "WOULD  YOU  LIKE  TO  MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?",X$ 
Again-FNYes(XS) 

Meas_up-1 
WHILE  Meas_up- 1 
SELECT  Again 
CASE  1 

LINPUT  "WOULD  YOU  LIKE  TO  CHAN6E  A  VARIABLE  (ENTER  Y  OR  N>?" ,X* 
Neu_va«FNYee( X$> 

Ano_change- 1 
WHILE  Ano_ change- 1 
SELECT  Neu_va 
CASE  1 

PRINT  "  VARIABLE  CURRENT  VALUE" 

PRINT  USING  "K" I " 1  TEMPERATURE  ", Temp, "CELSIUS" 

PRINT  USING  "K“ < "Z  ATOMOSPHERIC  PRESSURE  " .Aton_pres,"nt 
PRINT  USING  "K" » ”3  WATER  VAPOR  PRESSURE  " .Wat_pres."nb" 
PRINT  "  " 

PRINT  USING  "K"  <  "4  ANTENNA  DIAMETER  " ,Ant_diam,"  m.‘ 

PRINT  USING  "K'V'5  PULSE  LENGTH  ".Pulse,"  ns" 

PRINT  USING  “K" < “G  POWER  TRANSMITTED  " ,Pou_trana ,"  WATT* 
PRINT  USING  "K” i "7  ATMOSPHERIC  PROFILE  ".Profile 
PRINT  "  " 

PRINT  "ENTER  THE  NUMBER  OF  THE  VARIABLE  YOU  WISH  TO" 
PRINT  "CHANGE" 

INPUT  Var 
SELECT  Var 

CASE  1 

INPUT  "TEMPERATURE-" .Temp 
CASE  Z 


CURRENT  VALUE" 


'K"i 

"1 

TEMPERATURE 

Temp. "CELSIUS" 

K"« 

"Z 

ATOMOSPHERIC  PRESSURE 

" , Aton_pres . " mb" 

K"S 

"3 

WATER  VAPOR  PRESSURE 

" , Wat_pres , "mb" 

’K”  i 

"4 

ANTENNA  DIAMETER 

" .Ant.diam,"  m." 

•K"» 

"S 

PULSE  LENGTH 

" .Pulse,"  ma" 

'KM? 

"G 

POWER  TRANSMITTED  " ,Pou_trane ,"  WATTS" 

'K"» 

"7 

ATMOSPHERIC  PROFILE 

".Profile 

.Temp 
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INPUT  “ATMOSPHERIC  PRESSURE  IN  mb-" , Atomj»ree 
CASE  3 

INPUT  "WATER  VAPOR  PRESSURE  IN  mb-" ,Wat_pree 
CASE  4 

INPUT  "ANTENNA  DIAMETER  IN  m.-",Ant  dia* 

CASE  5 

INPUT  “PULSE  LENGTH  IN  ma«". Pulse 
CASE  6 

INPUT  "POUER  TRANSMITTED  IN  WATTS-" ,Pou_trana 
CASE  7 

PRINT  “  I  — >  Ct*2  PROFILE  OF  Z*(-1.1G>" 

PRINT  “  FROM  UALTERS/KUNDEL  198!" 

PRINT  "  2  — >  SAME  AS  ONE  BUT  WITH  FACTOR" 

PRINT  “  OF  TWO  TO  APPROXIMATE  LOOKING" 

PRINT  "  UP  A  THERMAL  PLOOM" 

PRINT  "  3  — >  Ct*2  PROFILE  OF  EXP( - .001*2 >" 

PRINT  "  FROM  UALTERS/KUNDEL  1981." 

PRINT  "  4  -->  CT*2  PROFILE  OF  ZA(-4/3>“ 

INPUT  "ENTER  NUMBER  OF  DESIRED  PROFILE" .Prof il 
IF  Profile-4  THEN 

INPUT  “HEIGHT  OF  INVERSION  IN  METERS-" . Inver 
ENO  IF 
CASE  ELSE 

PRINT  Var , " I S  NOT  ONE  OF  THE  OPTIONS" 

END  SELECT 

L INPUT  "MADE  ANOTHER  CHANGE (ENTER  Y  OR  N)?",X* 
Neu_va-FNYe»( X*> 

Meae_up-2 
CASE  2 

Meee_up-Z 
Ano_chanoe-2 
CASE  ELSE 

PRINT  “YOUR  CHOICES  WHERE  Y  OR  N  !!!!!” 

Ano_ change- 1 

LINPUT  "CHANGE  A  VARIABLE  (ENTER  Y  OR  N)?",X* 
Neu_va-FNYee( XS ) 

END  SELECT 
END  WHILE 
CASE  2 

Mess_up-2 
CASE  ELSE 

PRINT  "YOUR  CHOICES  WHERE  Y  OR  N  11!!!!" 

Meas_up- 1 

LINPUT  "MAKE  ANOTHER  RUN  (ENTER  Y  OR  N>?" ,X* 

Again-FNYes( X$) 

END  SELECT 
END  WHILE 
END  WHILE 

PRINT  "THAT’S  ALL.  FOLKS" 


I CALCULTE  THE  ATTENUATION 

sub  Attenuat ion(  Atom_pres , Attan , Freq , Temp , Uat_prea ) 

! THIS  SUBPROGRAM  CALCULATES  THE  ATTENUATION  OF  SOUND 
I  IN  AIR  BASED  UPON  EOUATIONS  IN  NEFF  1975 
I 

I  INPUT  ATMOSPHERIC  PRESSURE  IN  MILLIBARS 
!  FREQUENCY  OF  SOUND  WAVE  IN  HERTZ 

I  TEMPERATURE  IN  DEGREES  CELCIUS 

I  WATER- VAPOR  PRESSURE  IN  MILLIBARS 

! 

!  OUTPUT  ATTENUATION  IN  1 /METERS 
I 

I  VARIABLES 

!  Atom_pree  INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 

!  Atten  ATTENUATION  OF  ACOUSTIC  WAVE.  CALCULATED  IN 

I  SUBPROGRAM  ATTENUATION. 

!  Att.max  VARIABLE  IN  SUBPROGRAM  ATTENUATION.  IT  IS  THE 
!  ATTENUATION  AT  THE  FREQUENCY  OF  THE  MAXIMUM 

I  ATTENUATION  FOR  THE  INPUT  CONDITIONS. 

I  F  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.  IS  THE 

!  RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 

I  Fmax  FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 

I  ATTENUATION. 

!  Freq  INPUT  FREQUENCY  OF  ECHOSOUNOER. 

I  H  VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

!  Pstar  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 

!  CALCULATION. 

!  Temp  INPUT  TEMPERATURE  IN  OE6REES  CELSIUS. 

!  Tatar  INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

!  IN  SUBPROGRAM  ATTENUATION. 

I  Wat_prea  ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.  INPUT  BY 

1  OPERATOR. 

H-1O0«Wat_pree/Atom_prea 
Tatar-< 1 .5«Temp+49Z )/519 
Pstar-Atom_prea/ 1014 

Fmax-< 1®+66®0*H+44400«H*H)»Patar/Tetar“ .8 
Att_max-.O078*Fmax*Tstar '< -Z . 5 )«EXP( 7. 77*< 1  - 1 /Tatar ) ) 

F -Freq/ Fmax 

Atten-<  Att_max/304.B)«<  < . l8»F)*Z+<Z»F«F/< l+F«F> )*Z)*.5 

Atten-( Atten+1 .74E- 1®*Freq*Freq)/4. 35 

SUBEND 

! 

del  FNYea(XS) 

! THIS  FUNCTION  INTERPRETS  THE  OPERATORS  RESPONSES  TO  YES  NO  QUESTIONS 
! 

!  INPUT  X* 


OUTPUT 


FNYes 


>vvyvvv^»v  v  v  • /  rrm  r*.*  w 


1  VARIABLES 

I  Tempt  VARIABLE  STRING  USED  IN  FUNCTION  YES. 

I  Xt  STRIN6  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

I  RESPONSE  TO  YES  OR  NO  QUESTION. 

DIM  Tempt! 1 1 
TenptM  .n-TRIMt(Xt) 

SELECT  Tempt 
CASE  "Y","y" 

RETURN  I 
CASE  “N“ ,“n“ 

RETURN  Z 
CASE  ”  ” 

RETURN  1 
CASE  ELSE 

RETURN  -Z 
END  SELECT 

FNEND 


sub  Seni_log!  Freq.Xvar!  *•)  .L , Range!  *  > ,  J) 


THIS  SUBROUTINE  MAKES  A  SEMI -LOG  PLOT  OF  DATA  PASSED  FROM  THE 
MAIN  PROGRAM. 


INPUT 

X  AND  Y  VALUES  TO  BE  PLOTTED 
OUTPUT  SEMI -LOG  PLOT 


!  VARIABLES 
!  J 
!  L 

I  Xvar< • ) 


FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

VALUE  OF  ARRAY  OF  X  VALUES  TO  BE  PLOTTED  IN  SUBPROGRAM  Semi 


CLIP  ON 

FOR  J-1  TO  L  STEP  1 

PLOT  Xvar(J).LGT< Range! J>) 
NEXT  J 
LDIR  PI/4 
LORG  Z 
CSIZE  3..S 
LINE  TYPE  1 

MOVE  Xvar! J-1 ) ,LGT( Range! J-1) ) 
LABEL  USING  “K" !  Freq , "Hz" 

SUBEND 


I  oo 


r/.v.y.'/v 


APPENDIX  C 
COMPUTER  PROGRAM  3 


FULLER.  ROBERT  PR0G_3  EFFICIENCY 

10  SEP  8S 


•••••••PURPOSE******* 

THIS  PR06RAM  WILL  TAKE  PARAMETERS  OF  AN  ECHOSOUNOER  AND 
ESTIMATE  THE  RAN6E  AS  A  FUNCTION  OF  EFFICIENCY  OF  THE  THE  TRANSDUCER. 

THE  F0LL0UIN6  INPUTS  ARE  REQUIRED: 

ATMOSPHERIC  DATA 

1  >ATOMOSPHERIC  PRESSURE  in  millibars 

2 )  SELECT  ONE  OF  FOUR  ATMOSPHERIC  PROFILES  FOR  THE  TEMPERATURE 
STRUCTURE  PROFILE. 

a)FOR  ONE  PROFILE  THE  HEIGHT  OF  THE  INVERSION  LAYER  IS  ALSO  INPUT 

3)  TEMPERATURE  IN  DE6REES  CELSIUS 

4)  WATER  VAPOR  PRESSURE  IN  millibars 

ECHOSOUNOER  DATA 

5)  ANTENNA  DIAMETER  IN  METERS 

6)  FREQUENCY 

7) P0UER  TRANSMITTED  BY  ECHOSOUNDER  IN  WATTS 

8) PULSE  LENGTH  OF  THE  TRANSMITTED  ACOUSTIC  ENERGY 

THE  PROGRAM  OUTPUTS  THE  FOLLOWING  GRAPH  TO  AN  EXTERNAL  PLOTTER: 

1 )RAN6E( m. )  VERSUS  EFFICIENCY  OF  ECHOSOUNDER 


VARIABLES 


Again  TEST  VALUE  IN  WHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  MAKE  ANOUTHER  RUN. 

Ano_chang«  TEST  VALUE  IN  WHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  CHANGE  THE  VALUE  OF  INPUT  DATA  BEFORE  MAKING 
ANOUTHER  RUN. 

Ani.area  ECHOSOUNOER  ANTENNA  AREA  CALCULATED  FROM  INPUT  OF 
ANTENNA  DIAMETER. 

Ant_diam  INPUT  OF  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS. 

Atom_pres  INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 

Atten  ATTENUATION  OF  ACOUSTIC  WAVE.  CALCULATED  IN 

SUBPROGRAM  ATTENUATION. 

Bn  BANDWIDTH  OF  RECEIVER  NEEDED  TO  RECEIVE  SIGNALS 

DOPPLER  SHIFFED  BY  THREE  METER  PER  SECOND 
VERTICAL  VELOCITIES. 

C  VELOCITY  OF  SOUND  CALCULATION  AT  ZERO  CELSIUS. 

Ct  VARIABLE  USED  TO  SCALE  INVERSION  HEI6HT  IN 

TEMPERATURE  STRUCTURE  PROFILE  THREE. 

C3  VARIABLE  USED  IN  CALCULATING  THE  FOURTH  CteZ  PROFILE. 

CneZ( • )  ARRAY  OF  VALUES  OF  THE  ACOUSTIC  REFRACTIVE  INDEX 


PARAMETER.  CALCULATED  BASED  ON  SELECTION  OF 
PROFILE  FOR  TEMPERATURE  STRUCTURE  PARAMETER 
AND  ASSUMED  PROFILE  FOR  THE  VELOCITY  STRUCTURE 
PARAMETER. 

CteZ  VALUE  OF  THE  TEMPERATURE  STRUCTURE  PARAMETER. 

THE  OPERATOR  SELECTS  A  PROFILE  FROM  SEVERAL  PROVIDED. 

CveZ  VALUE  OF  THE  VELOCITY  STRUCTURE  PARAMETER. 

VALUES  BASED  ON  CALCULATION  USIN6  ASSUMED  DISSAPATION 
RATE. 

Dx  STEP  SIZE  FOR  X  AXIS  FOR  VARIOUS  PLOTS. 

Oy  STEP  SIZE  FOR  Y  AXIS  FOR  VARIOUS  PLOTS. 

Epsilon  DISSAPATION  RATE  USED  IN  CALCULATION  OF  CveZ. 

Eff  TRANSMISSION  EFFICIENCY  OF  ECHOSOUNOER. 

Eff  EFFICIENCY  OF  ECHOSOUNDER  WHEN  RECEIVING  ACOUSTIC 

BACK SCATTER. 

Exc.att  EXCESS  "ATTENUATION"  AT  GIVEN  RANGE. 

F  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.  IS  THE 

RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 

Fnax  FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 

ATTENUATION. 

Freq  INPUT  FREQUENCY  OF  ECHOSOUNDER. 

G  ANTENNA  EFFECTIVE  APERATURE  FACTOR. 

H  VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

I  MAIN  INDEX  FOR  LOOP  TO  CHECK  RETURNED  POWER  UNTIL 

POUER  WAS  LESS  THAN  MINIMUM  DETECTABLE  SIGNAL. 

IND  INDEX  USED  FOR  LOOP  FOR  DIFFERENT  EFFICIENCIES 

Interval  MAXIMUM  EFFECTIVE-SCATTERING  VOLUME  THICKNESS. 

Inver  HEIGHT  OF  INVERSION  LAYER. 

J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

K  WAVENUMBER 

L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

Meaa_up  TEST  VALUE  TO  CHECK  RESPONSES  WHEN  OPERATOR  IS  INPUTTING 
RESPONSES. 

N  TEST  VALUE  IN  CALCULATION  OF  COHERENCE  LENTH  IN  CALCULATION 

OF  EXCESS  ATTENUATION. 

Neu_va  TEST  VALUE  IN  CHECKING  IF  OPERATOR  WISHES  TO  CHANGE 
A  VARIABLE  BEFORE  A  NEW  RUN. 

Nolee  ASSUMED  MINIMUM  DETECTABLE  SIGNAL. 

Pow.back ( * )  POUER  BACKSCATTERED  FROM  GIVEN  RANGE. 

Pou_ret< • )  POUER  BACKSCATTERED  TO  ECHOSOUNDER  FROM  GIVEN  RANGE. 

Pou_trans  INPUT  OF  POUER  SUPPLIED  TO  TRANSDUCER  OF  ECHOSOUNDER. 

Profile  OPERATOR  INPUT  OF  CteZ  PROFILE  FROM  AVAILABLE  PROFILES. 

Pstar  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 

CALCULATION. 

Pulae  TRANSMITTED  PULSE  LENGTH  IN  MILLISECONDS.  INPUT  BY 

OPERATOR. 

R  RANGE  FOR  INPUT  TO  INTEGRAL  OF  EXCESS  ATTENUATION. 

R1  VARIABLE  USED  IN  SCALING  INVERSION  HEIGHT  FOR  PROFILE  THREE. 

Range< • )  ARRAY  OF  RANGE  VALUES. 

Remainder  REMAINDER  OF  MODULO  FUNCTION  USED  TO  DECREES  THE 


NUMBER  OF  PASSES  THROUGH  THE  INTEGRATION  FOR  EXCESS 
ATTENUATION. 

Rge  VALUE  OF  RANGE  IN  INTEGRAL  FOR  EXCESS  ATTENUATION. 

Rho  CORRELATION  LENGTH  USED  IN  CALCULATION  OF  EXCESS 

ATTENUATION. 

Signs  FRACTION  OF  ENERGY  BACK SCATTERED  FROM  GIVEN  RANGE. 

Spaed_sound  SPEED  OF  SOUND  IN  AIR  AT  INPUT  TEMPERATURE. 
Sunpou_back 

T  INPUT  TEMPERATURE  IN  DEGREES  KELVIN. 

Temp  INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

TenpC  VARIABLE  STRING  USED  IN  FUNCTION  YES. 


Tatar 


Uat  _pres 


Xnax 

Xnin 

Xrange 

YlabelS 

Ynax 

Ynin 

Yranga 

Ze 


INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 
IN  SUBPROGRAM  ATTENUATION. 

USED  TO  SELECT  THE  VARIABLE  THE  OPERATOR  UISHED  TO 
CHAN6E  BEFORE  MAKING  ANOTHER  RUN. 

ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.  INPUT  BY 
OPERATOR. 

THIRD  OROER  INDEX  USED  IN  VARIOUS  LOOPS. 

STRIN6  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 
RESPONSE  TO  YES  OR  NO  OUESTION. 

VALUE  OF  LARGEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

VALUE  OF  SMALLEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

VALUE  OF  RANGE  OF  X  VALUES  FOR  VARIOUS  PLOTS. 

LABEL  ON  Y  AXIS  PASSED  TO  SUBPROGRAM  Seni.log  FOR  PLOTTING. 
LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

SMALLEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

RANGE  OF  Y  VARIABLES  FOR  VARIOUS  PLOTS. 

INTERMEDIATE  VALUE  IN  CALCULATION  OF  EXCESS  ATTENUATION. 


I  DECLARE  VARIABLES 
INTEGER  I 

DIM  Pou_back! 1500) ,Pou_ret< 1500) .Range!  1500) 

DIM  CneZ< 1500) .Ran!  30) 

! 

PLOTTER  IS  705,"HPGL" 

LINE  TYPE  1 

! 

(INPUT  ATMOSPHERIC  DATA 

INPUT  "ENTER  TEMPERATURE  IN  DEGREES  CELCIUS" , Temp 
INPUT  "ENTER  ATMOSPHERIC  PRESSURE  IN  MILLIBARS” ,Ato«_preo 
INPUT  "ENTER  WATER  VAPOR  PRESSURE  IN  MILLIBARS" ,Wat_pres 


(INPUT  ECHOSOUNDER  DATA 

INPUT  "ENTER  ANTENNA  DIAM  IN  METERS" , Ant.di an 
INPUT  "ENTER  ECHOSOUNDER  FREQUENCY  IN  HERTZ", Freq 
INPUT  "ENTER  ECHOSOUNDER  PULSE  LENGTH  IN  MILLISECONDS-" .Pulse 
INPUT  "ENTER  POWER  TO  TRANSMITTER  IN  WATTS" ,Pou  trana 


I  03 
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G-. 40  ! ANTENNA  EFFECTIVE  APERATURE  FACTOR 

I 
t 

(SELECT  THE  PROFILE  FOR  THE  TEMPERATURE  STRUCTURE 
PRINT  "YOU  HAVE  A  CHOICE  OF  TEMPERATURE  STRUCTURE  PROFILES" 

PRINT  "PROFILE  TO  BE  USED." 

PRINT  "YOUR  SELECTIONS  ARE" 

PRINT  "  1 

PRINT  " 

PRINT  " 

PRINT  "  “ 

PRINT  "  Z 

PRINT  " 

PRINT  “ 

PRINT  "  “ 

PRINT  "  3 

PRINT  " 

PRINT  " 

PRINT  " 

PRINT  “ 

PRINT  "  " 

PRINT  "  4 

PRINT  “ 

! 

PRINT  "  " 

PRINT  “  “ 

Meea_up- 1 
UHILE  Meaa_up-1 

INPUT  "ENTER  THE  DESIRED  PROFILE  <1  OR  Z  OR  3  OR  4)", Profile 

IF  Profile-4  THEN 

INPUT  "ENTER  THE  HEIGHT  OF  THE  INVERSION  IN  METERS" . Inver 

END  IF 

IF  Profile-!  THEN 
Meaa_up-0 

ELSE 

IF  Profile-Z  THEN 
Meaa_up-0 

ELSE 

IF  Profile-3  THEN 
Meaa  up-0 

ELSE 

IF  Profile-4  THEN 

Meaa_up-0 

ELSE 

PRINT  Profile,"  WAS  NOT  ONE  OF  THE  OPTIONS!!!!" 

Meaa_up- 1 
END  IF 
END  IF 
END  IF 
END  IF 


A  TEWERARURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS 
PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  WHICH  " 
GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 

THE  TEMPERATURE  STRUCTURE  PROFILE  ABOVE  BUT" 

WITH  A  FACTOR  OF  TWO  TO  APPROXIMATE  LOOKING" 

UP  A  THERMAL  PLOOH" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 
PRESENTED  IN  WALTERS/ KUNDEL  1381  PAGE  398  FOR" 
NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH  HEIGHT  " 

IS  EXP< - . 001 *HEI6HT  ABOVE  65  METERS  AND  HEIGHT" 

TO  THE  -1.46  BELOW  65  METERS" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 

HEIGHT  TO  THE  -4/3  “ 
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END  WHILE 
OUTPUT  KBDi "K"  i 

I 

I 

Aoain-1 
WHILE  Again- 1 
FOR  J-1  TO  I 

Pou_back  <  J  )-0 
Pou_ret( J)-0 
Range< J)-0 
CneZ( J>-0 
NEXT  J 
! 

! CONVERT  TEMPERATURE  TO  KELVIN 
T-Tenp+Z73 
! 

! CALCULATE  SPEED  OF  SOUND  BASED  ON  INPUT  TEMPERATURE( CELCIUS ) 
Speed_sound-Z0.05*(  T)*.5 

! CALCULATE  THE  SPEED  OF  SOUND  AT  0  DEGREES  CELCIUS 
C-Z0.05*Z73*.S 
! 

1-1  ! INDEX  FOR  INNER  LOOP  (POWER  RETURNED  TO  NOISE) 

Ind-1  (INDEX  FOR  OUTER  LOOP  (EFFICIENCY) 

Pou_back ( 0 ) -0  (INITIALIZE  VARIABLE  FOR  POWER  BACKSCATERRED 

Ant_area-PI*( Ant_diam/Z )*Z  (ANTENNA  AREA 
Interval-!  Speed_sound*Pulse*  1  .  E-.D/2 
1 

UNITIALIZE  PLOTTING  PARAMETERS  SO  CAN  COMPARE  DURING  RUt'.  TO  SCALE 

(PLOTS 

Ynin-1000 

Ynax-0 

(CALCULATE  RETURNED  SIGNAL  POWER  UNTIL  IT  IS  LESS  THAN  BACKGROUND 
(FOR  EACH  FREOUENCY. 

Range( 0)»0 

FOR  Eff-.05  TO  .5  STEP  .0S 

Bn-Z»Freq*( 1 - 1 /( G/Speed_sound+ 1 ) ) (BANDWIDTH  FOR  3M/S  VERTICAL  VELOCITY 
Noise- 1 . 38E-Z3*Bn*( T )+Z . E- 14 ! MINIMUM  DETECTABLE  SIGNAL  AS 
I  JOHNSON  NOISE( NEGLIGABLE )+  ESTIMATED  BACKGROUND 
(CALCULATE  THE  ATTENUATION  COEFFICIENT  BASED  ON  ATMOSPHERIC  DATA 
CALL  Attenuation( Atom_pres , Atten.Freq.Tenp ,Wat_pres ) 
K-Z*PI»Freq/Speed_sound  ( WAVENUMBER 
1-1 

Sumpou_back -0 
REPEAT- 

Range(  I  )-Range(  I  - 1  )+Z 
SELECT  Profile 
CASE  1 

(THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
(AS  PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398. 

(WAS  FOR  MIDDAY  AND  CLEAR  WEATHER  ABOVE  THE 
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! TULAROSA  BASIN,  NEU  MEXICO. 

! AN  ADDITONAL  FACTOR  OF  1.8  UAS  INCLUDED  AS  PER  DIFFERENT 
(AVERAGING  TIME. 

CteZ“Z. 12*Range( I ) A< - 1 . 1 6) 

CASE  Z 

(THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
(AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  388. 

I UAS  FOR  MIDDAY  ANO  CLEAR  WEATHER  ABOVE  THE 
(TULAROSA  BASIN,  NEU  MEXICO. 

(AN  ADOI TONAL  FACTOR  OF  1.8  UAS  INCLUDED  AS  PER  DIFFERENT 
(AVERAGING  TIME  AND  AN  ADDITIONAL  FACTOR  OF  TWO  UAS 
(INCLUDED  TO  APPROXIMATE  LOOKING  UP  A  THERMAL  PLOOM 
CteZ-2‘2. 12*Range<  I )“(  -1 .  1G> 

CASE  3 

(THIS  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 

(PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 

(NIGHT  CONDITIONS.  THE  DEPENDENCE  UITH  HEIGHT  “ 

(IS  EXP< - .001 ‘HEIGHT  ABOVE  B5  METERS  AND  HEIGHT" 

(TO  THE  -1.4G  BELOU  GS  METERS" 

IF  Range< I ) <E5  THEN 

CteZ*75. 5»Range( I )"( -1 .46) 

ELSE 

Cte2-3.GGE-Z*EXP(-.001*Range< I ) ) 

END  IF 
CASE  4 

(CALCULATE  THE  TEMPERATURE  STRUCTURE  FACTOR  BASED  ON  ASSUMPTION 
!  THAT  IT’S  PROPORTIONAL  TO  HEIGHT  TO  THE  -4/3  POUER 

!  THAT  THE  SURFACE  VERTICAL  HEAT  FLUX  IS  .095  C.m./sec 

(THE  DEPENDENCE  UITH  HEIGHT  ABOVE  THE  INVERSION  I  PULLED  FROM 
(MY  LEFT  EAR 

(EQUATION  FROM  NEFF,  197S 
SELECT  Range< I )/Inver 
CASE  <.9 

C3-<< . 0Z4)*<  T)*( .B67) ) 

Ct eZ-C3* ( Range < I ) ) * ( -  1 . 33 ) 

CASE  .9  TO  1 
CteZ^CteZ 
Cl-CteZ 
R1 “Range< I ) 

CASE  1  TO  1.3 

CteZ" 10" ( <  LGT<  C3*Range( I ) )-LGT<  C3»R1 )>+LGT(C1  )) 
RZ»Range( I ) 

CZ-CteZ 
CASE  ELSE 

CteZ“C3*<  Range( I )"(-1 .33>-R2*<- 1 .33) )+CZ 
END  SELECT 
END  SELECT 

(CALCULATE  VELOCITY  STRUCTURE  FACTOR.  FORMULA  FROM  GAYNOR  77 
Epsilon-< . Z866/Range< I ) >•< 1 +.07«( Range( I ) ) *( .G))“( 1.5) 
CveZ“Z»EpsiIon*( .SG7) 


! CALCULATE  ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  FACTOR. 

(FORMULA  FROM  TETARSKI 

CneZ( I >-<  CteZ/<  Z . 98E+5 )  )+<  Cve2/( C«C  > ) 

(CALCULATE  THE  FRACTION  OF  POWER  BACKSCATTERED  FROM  INTERVAL. 

I  EQUATION  FROM  NEFF,  (975 
Signa-<  .0039«<K*<  1/3)  )*C-teZ)/<T>*2 
! 

(CALCULATE  THE  EXCESS  ATTENUATION.  EQUATION  FROM  CLIFFORD,  198® 
I  THE  EXCESS  ATTENUATION  IS  Exc.att 
IF  I <50  THEN 

Renal nder-0 

ELSE 

IF  I<150  THEN 

Renaindep-Range(  I )  MODULO  t® 

ELSE 

Rena inder- Range ( I )  MODULO  20 
END  IF 
END  IF 

IF  Renai nder-0  THEN 
Rho-C 
L-0 
H-0 

Rge-Range<I)  (CONSTANT  IN  INTEGRAL 
R-0 

FOR  J-0  TO  Z*I 

F-CneZ<INTU/2  +  1  )) 

F«F*< ( - R/ <  Rge  >*(  1 . 57 ) + ( R/ ( Rge ) ) * ( 1.E7)) 

IF  J>®  THEN 

IF  J<2*I  THEN 

IF  INT(J/Z)-J/Z  THEN 
L-L+F 
F-0 

ELSE 

H-H+F 
F-0 
END  IF 
END  IF 
END  IF 
Rho-Rho+F 
R-R+l 

NEXT  J 

Rho-Rho+4«L+Z*H 

Rho- ( ( <  Rho« . 33)»K*K*(.4G)*(-.G)) 

N*(  Ant_dian/Rho ) A 2 
IF  N<- I  THEN 

Ze-I /( 1+N) 

ELSE 

Ze-1 . 5/( 1+N) 

(STEP  OF  1.5— >SEE  CLIFFORD  1980 
END  IF 


Exc_att-Ze*Ze 


ELSE 

Exc_at+-Exc_att 
END  IF 

! 

I  CALCULATE  THE  POWER  BACK SCATTERED 

Pou_back(  I )-<Pou_trans*Eff-Sunpou_back  >*EXP( -Atten*Range( I ) ) 
Pou_back<  I )-Pou_back( I )*Interval*Exc_att«Signa 
Sunpou_back-Sunpou_back+Pou_back< I ) 

(CALCULATE  THE  POWER  BACKSCATTERED  TO  THE  ANTENNA 
Pou_ret(  I  )-Pou_back<  I  )*EXP(  -A't-ten*Range<  I )  )*Ant_area*G*Ef  f  / Range ( I )  *2 
PRINT  "RANGE-"  .Ranged) 

PRINT  "POWER  RETURNED-", Pou.retd) 

I-I  +  1 

UNTIL  Pou_reid-1  >>0  AND  Pou_retd- 1  X-Noise 
Ran( Ind)-Range< I-Z ) 

IF  Range( I-Z  >>Ynax  THEN 
Ynax-Range( I-Z ) 

END  IF 

IF  Range<  I-Z  XYinin  THEN 
Ynin-Range< I-Z ) 

END  IF 
Ind-Ind+1 
Again- I 
FOR  J-(  TO  I 

Pou_back( J)-0 
Pou_ret<  J)-0 
Range<  J )-0 
CneZ< J)-0 
NEXT  J 
NEXT  Eff 
I 
! 

OUTPUT  KBO; "K" t 
PRINT  "INPUT  CONDITIONS" 

! 

PRINT  "  " 

PRINT  USING  "K"; "TEMPERATURE- 
PRINT  USING  "K"j "ATMOSPHERIC  PRESSURE- 
PRINT  USING  "K'VUATER  VAPOR  PRESSURE- 
PRINT  USING  "K" i "PULSE  LENGTH- 
PRINT  USING  "K"; "ANTENNA  DIAMETER- 
PRINT  USING  "K“ 5 "ECHOSOUNDER  FREQUENCY- 
PRINT  USING  "K"j "POWER  TRANSMITTED- 
PRINT  "TEMPERATURE  STRUCTURE  PROFILE  USED  ".Profile 
IF  Profile-4  THEN 

PRINT  "INVERSION  HEIGHT  ".Inver 
END  IF 

SELECT  Profile 
CASE  ( 


",Tenp,"CELCIUS" 

" , Aton_pres , "mb" 

"  ,Wa't_pres,"mb" 

" .Pulse, "ns" 
",Ant_dian,"  n." 

" ,Freq,"Hz” 
",Pou_trans,"  WATTS" 


PRINT  "  1  ft  TEMPERftRURE  STRUCTURE  PROFILE  BftSED  ON  DftTft 

PRINT  “  AS  PRESENTED  IN  WALTERS/ KUNOEL  1981  PAGE  398 

PRINT  "  WHICH  GIVES  ft  HEIGHT  TO  THE  -1.16  PROFILE" 

CASE  Z 

PRINT  "2ft  TEMPERftRURE  STRUCTURE  PROFILE  BASED  ON  DATA 
PRINT  "  AS  PRESENTED  IN  WfiLTERS/KUNDEL  1981  PftGE  398 

PRINT  “  WHICH  GIVES  ft  HEIGHT  TO  THE  -1.16  PROFILE- 

PRINT  "  BUT  WITH  ft  FACTOR  OF  TWO  TO  APPROXIMATE" 

PRINT  “  LOOKING  UP  ft  THERMAL  PLOOM" 

CASE  3 

PRINT  "  3  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA 

PRINT  "  AS  PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398“ 

PRINT  "  FOR  NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH" 

PRINT  "  HEIGHT  IS  EXP<-.001*HEIGHT  ABOVE  65  METERS" 

PRINT  "  AND  HEIGHT  TO  THE  -1.46  BELOW  65  METERS" 

CASE  4 

PRINT  M  4  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 

PRINT  "  HEIGHT  TO  THE  -4/3  “ 

END  SELECT 
! 

PRINT  USING  "K" s "THE  MINIMUM  DETECTABLE  SIGNAL  WAS  SET  AT" .Noise, "WATTS 
! 

PRINT  "  " 

! 

! 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

OUTPUT  KBO; "K“ ; 

I 

GRAPHICS  ON 
VIEWPORT  15,120,10.70 
I 
! 

I  PLOT  RANGE  VERSUS  FREQUENCY 
XKiin"® 

Xaax* . 5 

Y*i n-PROUND< Ymi  n-  50 ,  Z ) 

Ywax“PROUND< Ynax+50 , 2 ) 

Ox-. 05 
Dy- 100 

XranQe*Xwax-X«in 
Yrange-Ywax-Ynin 
WINDOW  Xnin.Xnax ,Ynin, Ynax 
AXES  Dx .Dy.Xnin, Ynin.Z , 1 
CLIP  OFF 

I 

! LABEL  PLOT 
CSIZE  4,. 6 
LOIR  0 
LORG  5 


MOVE  . 5»Xrange , Ynax* 1 . 05 

LABEL  “RANGE  VERSUS  TRANSDUCER  EFFICIENCY" 

! LABEL  HORIZONTAL  AXES 
LOIR  0 
LORG  5 

FOR  J-0  TO  Xnax  STEP  Dx 
CSIZE  4, .5 

MOVE  J.Ynin- .05*Yrange 
LABEL  J 
NEXT  J 

MOVE  .5»Xrange,Ynin-. 1*Yrange 

CSIZE  4..G 

LABEL  "EFFICIENCY" 

! LABEL  VERTICAL  AXES 
LORG  6 

FOR  J-Ynin  TO  Yixax  STEP  Dy 
CSIZE  4..G 

MOVE  Xmin-.0025*Xrange, J 
LABEL  USING  "K"jJ 
NEXT  J 
LDIR  PI/Z 
LORG  6 

MOVE  Xmin- . 1 S«Xranga , . 5*Yrange+Y«in 
CSIZE  4 , . G 

LABEL  "RANGE  (METERS)" 

CLIP  ON 
Ind-1 

FOR  Eff “.05  TO  .5  STEP  .05 
PLOT  Eff.Ran(Ind) 

Ind-Ind+1 
NEXT  Eff 
! 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT  KBD!"K"! 

I 

! 

LINPUT  "WOULD  YOU  LIKE  TO  MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?",X* 
Again-FNYea(X*) 

Mess_up*1 
WHILE  Meaa_up*l 
SELECT  Again 
CASE  I 

LINPUT  "WOULD  YOU  LIKE  TO  CHANGE  A  VARIABLE  (ENTER  Y  OR  N)?",X* 
Neu_va“FNYea( X*) 

Ano_change- 1 
WHILE  Ano_change"1 
SELECT  Neu_va 
CASE  t 


1 1  O 


PRINT  “  VARIABLE  CURRENT  VALUE" 

PRINT  USING  TEMPERATURE  " .Temp, "CELSIUS" 

PRINT  USIN6  "K“ ; "2  ATOMOSPHERIC  PRESSURE  " ,Atom_pres . "mb 
PRINT  USIN6  “K"i“3  WATER  VAPOR  PRESSURE  “  ,Wa-t_pres,"mb“ 
PRINT  “  " 

PRINT  USING  "K” j "4  ANTENNA  DIAMETER  ",Ant_diam,"  m." 

PRINT  USIN6  “K" i "S  PULSE  LENGTH  ".Pulse,"  ms" 

PRINT  USING  "K" i“G  POUER  TRANSMITTED  “ ,Pou_trans UATTS 
PRINT  USIN6  "K‘V7  ATMOSPHERIC  PROFILE  ".Profile 
PRINT  "  " 

PRINT  "ENTER  THE  NUMBER  OF  THE  VARIABLE  YOU  WISH  TO" 

PRINT  "CHANGE" 

INPUT  Var 
SELECT  Var 

CASE  1 

INPUT  “TEMPERATURE-". Temp 
CASE  Z 

INPUT  "ATMOSPHERIC  PRESSURE  IN  «ib-" .  Atom_pres 
CASE  3 

INPUT  "WATER  VAPOR  PRESSURE  IN  mb-" ,Wat_pres 
CASE  4 

INPUT  "ANTENNA  DIAMETER  IN  m. -" .Ant_diam 
CASE  S 

INPUT  "PULSE  LENGTH  IN  ms-" .Pulse 
CASE  6 

INPUT  "POWER  TRANSMITTED  IN  WATTS-", Pou  trans 
CASE  7 

PRINT  "  1  — >  Ct*Z  PROFILE  OF  Z*<-I.16>" 

PRINT  "  FROM  WALTERS/KUNDEL  1981" 

PRINT  "  Z  — >  SAME  AS  ONE  BUT  WITH  FACTOR- 
PRINT  "  OF  TWO  TO  APPROXIMATE  LOOKING" 

PRINT  "  UP  A  THERMAL  PLOOM" 

PRINT  "  3  — >  Ct*Z  PROFILE  OF  EXP< -.0®1«Z>" 

PRINT  “  FROM  WALTERS/KUNDEL  1981." 

PRINT  "  4  — >  CT'Z  PROFILE  OF  Z‘<-4/3>" 

INPUT  "ENTER  NUMBER  OF  DESIRED  PROFILE" .Prof il 
IF  Profile-4  THEN 

INPUT  "HEIGHT  OF  INVERSION  IN  METERS-" . Inver 
END  IF 
CASE  ELSE 

PRINT  Var. "IS  NOT  ONE  OF  THE  OPTIONS" 

ENO  SELECT 

LINPUT  "MADE  ANOTHER  CHANGE ( ENTER  Y  OR  N)?".X* 
Neu_va-FNYes( X$) 

Mess_up-Z 
CASE  Z 

Mess_up-Z 
Ano_change-Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  WHERE  Y  OR  N  !!!!!" 


Ano_change* 1 

LINPUT  "CHANGE  A  VARIABLE  (ENTER  Y  OR  N)?",X* 
Neu_va-FNYee(  XS  > 

END  SELECT 
ENO  WHILE 
CASE  Z 

Meee_up"2 
CASE  ELSE 

PRINT  "YOUR  CHOICES  UHERE  Y  OR  N  I !  !  !!!  “ 

Mess_up“ I 

LINPUT  "MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?“,X* 

Again-FNYeef  X$) 

ENO  SELECT 
ENO  WHILE 
END  WHILE 

PRINT  "THAT’S  ALL.  FOLKS" 
end 


(CALCULTE  THE  ATTENUATION 

sub  Attenuat ion(  At on_pres . At ten , Freq , Tenp , Wat_pres ) 

! THIS  SUBPROGRAM  CALCULATES  THE  ATTENUATION  OF  SOUND 
! IN  AIR  BASED  UPON  EQUATIONS  IN  NEFF  <975 
! 

!  INPUT  ATMOSPHERIC  PRESSURE  IN  MILLIBARS 
I  FREQUENCY  OF  SOUND  WAVE  IN  HERTZ 

!  TEMPERATURE  IN  DEGREES  CELCIUS 

!  WATER-VAPOR  PRESSURE  IN  MILLIBARS 

I 

!  OUTPUT  ATTENUATION  IN  1 /METERS 
! 

! VARIABLES 

I  Aton_pres  INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 

ATTENUATION  OF  ACOUSTIC  WAVE.  CALCULATED  IN 
SUBPROGRAM  ATTENUATION. 

VARIABLE  IN  SUBPROGRAM  ATTENUATION.  IT  IS  THE 
ATTENUATION  AT  THE  FREQUENCY  OF  THE  MAXIMUM 
ATTENUATION  FOR  THE  INPUT  CONDITIONS. 

VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.  IS  THE 
RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 
FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 
ATTENUATION. 

INPUT  FREQUENCY  OF  ECHOSOUNDER. 

VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 
VARIABLE  USEO  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 
CALCULATION. 

INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 
IN  SUBPROGRAM  ATTENUATION. 

ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.  INPUT  BY 


! 


t  Atten 
! 

!  Att_nax 
! 

I  F 

!  Fnax 

!  Freq 
!  H 

!  Pstar 

I 

I  Tenp 
!  Tatar 

!  Wat_pree 


12 


I  OPERATOR. 

H- 1 0®*Ua-t_pres/ft'to*_pre8 
Ts-tar-(  1 .8«Te«p449Z  )/5IS 
Patar-Ato«_prea/ 1 01 4 

Fnax-< 10+G6Q0*H+44400*H*H)*Patar/Tatar* . 8 
Att_*ax-.0078*F*ax#Tatar*<-Z.S>*EXP<7.77*< 1-1/Tetar>> 

F“Freq/F«ax 

Atten-(fl-tt_«ax/304.8)*<( . 18*F>*Z+(Z*F*F/<  1+F*F)  )*Z)*.S 

Atten«<  Atten-H .74E- 10*Freq*Freq)/4. 35 

SUBENO 

! 

def  FNYee<  XC) 

'THIS  FUNCTION  INTERPRETS  THE  OPERATORS  RESPONSES  TO  YES  NO  QUESTIONS 

I 

I  INPUT  X* 

! 

I  OUTPUT  FNYea 

I  VARIABLES 

!  Tenpt  VARIABLE  STRIN6  USED  IN  FUNCTION  YES. 

!  X*  STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

!  RESPONSE  TO  YES  OR  NO  QUESTION. 

DIM  Tenp«( 1 ] 

Tenp*t  1  .  H-TRIM*<  X*> 

SELECT  Tenp« 

CASE  "Y-.-y" 

RETURN  1 
CASE  “N" ,"n" 

RETURN  2 
CASE  “  " 

RETURN  1 
CASE  ELSE 

RETURN  -2 
END  SELECT 


FULLER.  ROBERT 


APPENDIX  D 

COMPUTER  PROGRAM  4 
PR0G_4  ANTENNA  SIZE 

I©  SEP  85 


•••••••PURPOSE* •*•**• 

THIS  PROGRAM  WILL  TAKE  PARAMETERS  OF  AN  ECHOSOUNOER  AND 
ESTIMATE  THE  RANGE  AS  A  FUNCTION  OF  ANTENNA  SIZE  OF  THE  THE  ECHOSOUNDER. 
THE  FOLLOUING  IW»UTS  ARE  REQUIRED: 

ATMOSPHERIC  DATA 

1 )ATOMOSPHERIC  PRESSURE  in  millibars 

Z ) SELECT  ONE  OF  FOUR  ATMOSPHERIC  PROFILES  FOR  THE  TEMPERATURE 
STRUCTURE  PROFILE. 

a)FOR  ONE  PROFILE  THE  HEIGHT  OF  THE  INVERSION  LAYER  IS  ALSO  INPUT 

3 )  TEMPERATURE  IN  DEGREES  CELSIUS 

4)  WATER  VAPOR  PRESSURE  IN  millibars 

ECHOSOUNDER  OATA 

5)  FREQUENCY 

G) POWER  TRANSMITTED  BY  ECHOSOUNDER  IN  WATTS 
7)PULSE  LENGTH  OF  THE  TRANSMITTED  ACOUSTIC  ENERGY 

THE  PROGRAM  OUTPUTS  THE  FOLLOWING  GRAPH  TO  AN  EXTERNAL  PLOTTER: 

1 ) RANGE ( m. )  VERSUS  ANTENNA  DIAMETER  FOR  VARIOUS  FREQUENCIES 


VARIABLES 


I  Again 

! 

!  Ano_change 

,  !  Ant_area 

! 

I  Ant_diam 
!  Atom _pre s 
!  Atton 

| 

!  Bn 
! 

I 

!  C 
!  Cl 
( 

!  C3 

!  CneZ(*) 

| 

i 


TEST  VALUE  IN  WHILE  LOOP  TO  DETERMINE  IF  OPERATOR 
WISHES  TO  MAKE  ANOUTHER  RUN. 

TEST  VALUE  IN  WHILE  LOOP  TO  DETERMINE  IF  OPERATOR 
WISHES  TO  CHANGE  THE  VALUE  OF  INPUT  DATA  BEFORE  MAKING 
ANOUTHER  RUN. 

ECHOSOUNDER  ANTENNA  AREA  CALCULATED  FROM  INPUT  OF 
ANTENNA  DIAMETER. 

INPUT  OF  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS. 

INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 

ATTENUATION  OF  ACOUSTIC  WAVE.  CALCULATED  IN 
SUBPROGRAM  ATTENUATION. 

BANOWIOTH  OF  RECEIVER  NEEDED  TO  RECEIVE  SIGNALS 
DOPPLER  SHIFFED  BY  THREE  METER  PER  SECOND 
VERTICAL  VELOCITIES. 

VELOCITY  OF  SOUND  CALCULATION  AT  ZERO  CELSIUS. 

VARIABLE  USED  TO  SCALE  INVERSION  HEIGHT  IN 
TEMPERATURE  STRUCTURE  PROFILE  THREE. 

VARIABLE  USED  IN  CALCULATING  THE  FOURTH  CtoZ  PROFILE. 

ARRAY  OF  VALUES  OF  THE  ACOUSTIC  REFRACTIVE  INDEX 
PARAMETER.  CALCULATED  BASED  ON  SELECTION  OF 


!  PROFILE  FOR  TEMPERATURE  STRUCTURE  PARAMETER 

I  AND  ASSUMED  PROFILE  FOR  THE  VELOCITY  STRUCTURE 

I  PARAMETER. 

!  CteZ  VALUE  OF  THE  TEMPERATURE  STRUCTURE  PARAMETER. 

!  THE  OPERATOR  SELECTS  A  PROFILE  FROM  SEVERAL  PROVIDED. 

I  CveZ  VALUE  OF  THE  VELOCITY  STRUCTURE  PARAMETER. 

I  VALUES  BASED  ON  CALCULATION  USING  ASSUMED  OISSAPATION 

I  RATE. 

I  Dx  STEP  SIZE  FOR  X  AXIS  FOR  VARIOUS  PLOTS. 

'  Dy  STEP  SIZE  FOR  Y  AXIS  FOR  VARIOUS  PLOTS. 

I  Epsilon  OISSAPATION  RATE  USED  IN  CALCULATION  OF  CveZ. 

!  Eff  TRANSMISSION  EFFICIENCY  OF  ECHOSOUNDER. 

I  Eff  EFFICIENCY  OF  ECHOSOUNDER  UHEN  RECEIVING  ACOUSTIC 

I  BACK SCATTER. 

!  Exc_att  EXCESS  “ATTENUATION"  AT  GIVEN  RANGE. 

!  F  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.  IS  THE 

!  RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 

!  Fmax  FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 

!  ATTENUATION. 

!  Freq  INPUT  FREQUENCY  OF  ECHOSOUNDER. 

!  G  ANTENNA  EFFECTIVE  APERATURE  FACTOR. 

!  H  VARIA8LE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

!  I  MAIN  INDEX  FOR  LOOP  TO  CHECK  RETURNED  POWER  UNTIL 

!  POWER  WAS  LESS  THAN  MINIMUM  DETECTABLE  SICNAL. 

!  IND  INDEX  USED  FOR  LOOP  FOR  DIFFERENT  EFFICIENCIES 

!  Interval  MAXIMUM  EFFECTIVE-SCATTERING  VOLUME  THICKNESS. 

!  Inver  HEIGHT  OF  INVERSION  LAYER. 

!  J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

!  K  WAVENUMBER 

!  L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

!  Meas_up  TEST  VALUE  TO  CHECK  RESPONSES  UHEN  OPERATOR  IS  INPUTTING 
!  RESPONSES. 

!  N  TEST  VALUE  IN  CALCULATION  OF  COHERENCE  LENTH  IN  CALCULATION 

I  OF  EXCESS  ATTENUATION. 

!  Neu_va  TEST  VALUE  IN  CHECKING  IF  OPERATOR  WISHES  TO  CHANGE 
!  A  VARIABLE  BEFORE  A  NEW  RUN. 

I  Noise  ASSUMED  MINIMUM  DETECTABLE  SIGNAL. 

!  Pow_back ( * )  POWER  BACKSCATTERED  FROM  GIVEN  RANGE. 

!  Pou_ret<*>  POWER  BACKSCATTERED  TO  ECHOSOUNDER  FROM  GIVEN  RANGE. 

!  Pou_trans  INPUT  OF  POWER  SUPPLIED  TO  TRANSDUCER  OF  ECHOSOUNDER. 

!  Profile  OPERATOR  INPUT  OF  CteZ  PROFILE  FROM  AVAILABLE  PROFILES. 

!  Pstar  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 

I  CALCULATION. 

!  Pulse  TRANSMITTED  PULSE  LENGTH  IN  MILLISECONDS.  INPUT  BY 

!  OPERATOR. 

!  R  RANGE  FOR  INPUT  TO  INTEGRAL  OF  EXCESS  ATTENUATION. 

!  R1  VARIABLE  USED  IN  SCALING  INVERSION  HEIGHT  FOR  PROFILE  THREE. 

!  Range< * )  ARRAY  OF  RANGE  VALUES. 

!  Remainder  REMAINDER  OF  MODULO  FUNCTION  USED  TO  DECREES  THE 
!  NUMBER  OF  PASSES  THROUGH  THE  INTEGRATION  FOR  EXCESS 


I  I  5 


!  ATTENUATION. 

!  Rge  VALUE  OT  RANGE  IN  INTEGRAL  FOR  EXCESS  ATTENUATION. 

!  Rho  CORRELATION  LENGTH  USED  IN  CALCULATION  OF  EXCESS 

!  ATTENUATION. 

!  Signs  FRACTION  OF  ENERGY  BACK SCATTERED  FROM  GIVEN  RANGE. 

!  Speed_BOund  SPEED  OF  SOUNO  IN  AIR  AT  INPUT  TEMPERATURE. 

I  Sunpow_back 

!  T  INPUT  TEMPERATURE  IN  DEGREES  KELVIN. 

!  Tenp  INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

!  Tenp*  VARIABLE  STRING  USED  IN  FUNCTION  YES. 

!  Tatar  INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

!  IN  SUBPROGRAM  ATTENUATION. 

!  Var  USED  TO  SELECT  THE  VARIABLE  THE  OPERATOR  UISHED  TO 

!  CHANGE  BEFORE  MAKING  ANOTHER  RUN. 

!  Wat_pres  ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.  INPUT  BY 
!  OPERATOR. 

!  X  THIRD  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 

!  X$  STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

!  RESPONSE  TO  YES  OR  NO  QUESTION. 

!  Xnax  VALUE  OF  LARGEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

!  Xnin  VALUE  OF  SMALLEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

!  Xrange  VALUE  OF  RANGE  OF  X  VALUES  FOR  VARIOUS  PLOTS. 

!  Ylabel*  LABEL  ON  Y  AXIS  PASSED  TO  SUBPROGRAM  Se«i_log  FOR  PLOTTING. 

!  Ynax  LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

!  Ynin  SMALLEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

!  Yrange  RANGE  OF  Y  VARIABLES  FOR  VARIOUS  PLOTS. 

!  Ze  INTERMEDIATE  VALUE  IN  CALCULATION  OF  EXCESS  ATTENUATION. 


'DECLARE  VARIABLES 
INTEGER  I 

DIM  Pou_back< 1500) ,Pou_ret< 1500) , Rangel  1500) 

DIM  CneZd 500)  ,Ran<  50,  G) 

! 

PLOTTER  IS  705, "HPGL" 

LINE  TYPE  1 
! 

! INPUT  ATMOSPHERIC  DATA 

INPUT  "ENTER  TEMPERATURE  IN  DEGREES  CELCIUS” , Tenp 
INPUT  "ENTER  ATMOSPHERIC  PRESSURE  IN  MILLIBARS" , Atom _pres 
INPUT  "ENTER  UATER  VAPOR  PRESSURE  IN  MILLIBARS" ,Wat_pres 


(INPUT  ECHOSOUNDER  DATA 

INPUT  "ENTER  ECHOSOUNDER  PULSE  LENGTH  IN  MILLISECONDS-" .Pulse 
INPUT  "ENTER  POWER  TO  TRANSMITTER  IN  WATTS" ,Pou_trans 
G-.40  (ANTENNA  EFFECTIVE  APERATURE  FACTOR 

Eff-.Z5 


! TULAROSA  BASIN,  NEU  MEXICO. 

(AN  ADDITONAL  FACTOR  OF  1.8  UAS  INCLUDED  AS  PER  DIFFERENT 
! AVERAGIN6  TIME. 

Cte2-Z . 1 Z*Range< I ) *  < -  1 . 1 6 ) 

CASE  2 

(THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
! AS  PRESENTED  IN  WALTERS/ KUNDEL  1981  PAGE  398. 

I UAS  FOR  MIDDAY  AND  CLEAR  WEATHER  ABOVE  THE 
( TULAROSA  BASIN,  NEU  MEXICO. 

! AN  ADDITONAL  FACTOR  OF  1.8  WAS  INCLUDED  AS  PER  DIFFERENT 
(AVERAGING  TIME  AND  AN  ADDITIONAL  FACTOR  OF  TWO  WAS 
(INCLUDED  TO  APPROXIMATE  LOOKING  UP  A  THERMAL  PLOOM 
CteZ*2*2 . 12*Range( I )*< - l . 1 6) 

CASE  3 

(THIS  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 

(PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398  FOR" 

(NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH  HEIGHT  ” 

(IS  EXP<-.001«HEIGHT  ABOVE  B5  METERS  AND  HEIGHT" 

(TO  THE  -1.46  BELOW  GG  METERS" 

IF  Range<  I  KS5  THEN 

Cte2-7S.5*Range< I)*<-1 .46) 

ELSE 

Cte2-3.66E-2*EXP<-.001*Ranged>> 

END  IF 
CASE  4 

(CALCULATE  THE  TEMPERATURE  STRUCTURE  FACTOR  BASED  ON  ASSUMPTION 
!  THAT  IT’S  PROPORTIONAL  TO  HEIGHT  TO  THE  -4/3  POWER 

!  THAT  THE  SURFACE  VERTICAL  HEAT  FLUX  IS  .095  C.n./sec 

f THE  DEPENDENCE  WITH  HEIGHT  ABOVE  THE  INVERSION  I  PULLED  FROM 
(MY  LEFT  EAR 

! EOUATION  FROM  NEFF,  1975 
SELECT  Range< I )/Inver 
CASE  <.9 

C3-( < -0Z4)*( T)“< .667) > 

CteZ-C3*<  Rangel  I >  >*< -1 . 33 > 

CASE  .9  TO  1 
CteZ-CteZ 
Cl-CteZ 
R1 “Range< I ) 

CASE  1  TO  1.3 

CteZ* J0*<  <  LGT<  C3* Rangel  I ) >-LGT( C3»R1 ) )+LGT( Cl )  ) 
RZ“Range( I ) 

CZ-C-teZ 
CASE  ELSE 

CteZ-C3*<  Ranged )*(-! . 33>-RZ*<- 1 .33) >+CZ 
END  SELECT 
END  SELECT 

(CALCULATE  VELOCITY  STRUCTURE  FACTOR.  FORMULA  FROM  GAYNOR  77 
Epeilon»( .2B66/Range( I ) )»< l  +  .07«<  Range(  I ) )*(  .6) )"(  1.5) 

C  veZ  »Z  *Eps i 1  on" ( .667) 


I  CALCULATE  ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  FACTOR. 

I  FORMULA  FROM  TETARSKI 

CneZ( I )-<  CteZ/<  2.98E+5) >+< CveZ/( C*C> ) 

I  CALCULATE  THE  FRACTION  OF  POUER  BACKSCATTERED  FROM  INTERVAL, 
t EQUATION  FROM  NEFF,  1975 
Sigma-(  .0O39»<ir<  1/3) >*CteZ )/< T> A2 
! 

I  CALCULATE  THE  EXCESS  ATTENUATION.  EQUATION  FROM  CLIFFORD.  1980 
(THE  EXCESS  ATTENUATION  IS  Exc_att 
IF  I <50  THEN 

Remainder -0 

ELSE 

IF  I <150  THEN 

Remainder-Range( I )  MODULO  10 

ELSE 

Remainder-Range< I )  MODULO  20 
END  IF 
END  IF 

IF  Remainder-0  THEN 
Rho-0 
L-0 
H-0 

Roe-Ranged)  !  CONSTANT  IN  INTEGRAL 
R-0 

FOR  J-0  TO  2*1 

F-CneZ< INT( J/2+1 ) ) 

F-F*< 1 -R/( Rge) A< 1 ,E7)+( R/( Rge ) )  A( 1.G7)) 

IF  J>0  THEN 

IF  J<*!*I  THEN 

IF  INT( J/2  >«J/Z  THEN 
L-L+F 
F-0 

ELSE 

H-H+F 
F-0 
END  IF 
END  IF 
END  IF 
Rho-Rho+F 
R-R+1 

NEXT  J 

Rho-Rho+4*L4Z*H 

Rho-<  <  (Rho*.33)*K*K*1 .45>*{-.B) > 

N-<  Ant_diam/Rho)‘Z 
IF  N<-7  THEN 

Ze-1 /( 1+N) 

ELSE 

Ze-1 .5/( 1+N) 

!  STEP  OF  1.5— >SEE  CLIFFORD  1980 
END  IF 
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Exc_att-Ze*Ze 

ELSE 

Exc_ait-Exc_att 
END  IF 

I  CALCULATE  THE  POWER  BACKSCATTERED 

Pou_back(  I >-<Pou_trane*£f f-Sunpou_back )*EXP( -Atten*Range( I ) ) 
Pou_back(  I  )-Pou_back<  I  )*In-terval*Exc_att*SiQna 
Sunpou_back-Suflpou_back+Pou_back(  I ) 

(CALCULATE  THE  POWER  BACKSCATTERED  TO  THE  ANTENNA 
Pou_ret(  I  )-Pou_back( I  >*EXP(  -Atten«Range(  I )  )*Ant_area*6*Ef  f / Range ( I )  "Z 
PRINT  "RANGE*" , Rango< I ) 

PRINT  "POWER  RETURNED-". Pou_ret( I > 

PRINT  "FREQUENCY-" .Freq 
I-I  +  1 

UNTIL  Pou_ret(  1-1 )>0  AND  Pou_ret< I- 1 ) <*Noise 
Ran< Ind, Ifreq)-Range< I-Z ) 

IF  Ran< Ind,Ifreq)>Ymax  THEN 
Ymax-Ran( Ind, If roq) 

END  IF 

IF  Ran(  Ind.IfreqXYnin  THEN 
Y«in-Ran( Ind.Ifreq) 

END  IF 
Ind-Ind+1 
Aoain-1 
FOR  J-0  TO  I 

Pou_back<  J)-0 
Pou_ret<  J>«0 
Range<  J ) -0 
CneZ( J)-0 
NEXT  J 

NEXT  Ant_dian 
Ifreq-Ifreq+1 
Ind-1 
NEXT  Froq 

! 

j 

OUTPUT  KBD;"K"s 

PRINT  "INPUT  CONDITIONS" 

! 

PRINT  "  “ 

PRINT  USING  "K" ; "TEMPERATURE-  " ,Temp,"CELCIUS" 

PRINT  USING  "K"; "ATMOSPHERIC  PRESSURE-  " , Aton.pree . "nb" 

PRINT  USING  "K”; "WATER  VAPOR  PRESSURE-  " , Wat_pres , "«b" 

PRINT  USING  "K"! "PULSE  LENGTH-  " .Pul se . "ns" 

PRINT  USING  "K"! "TRANSDUCER  EFFICIENCY-  ",Eff 

PRINT  USING  "K" ; "POWER  TRANSMITTED-  ”  ,Pow  trane,"  WATTS" 

PRINT  "TEMPERATURE  STRUCTURE  PROFILE  USED  ".Profile 
IF  Profile-4  THEN 

PRINT  "INVERSION  HEIGHT  ".Inver 


END  IF 

SELECT  Profile 
CASE  1 

PRINT  "  1 

PRINT  “ 

PRINT  “ 

CASE  2 

PRINT  “  2 

PRINT  “ 

PRINT  “ 

PRINT  “ 

PRINT  " 

CASE  3 

PRINT  "  3 

PRINT  “ 

PRINT  " 

PRINT  ” 

PRINT  “ 

CASE  4 

PRINT  "  4 

PRINT  M 
END  SELECT 
! 

PRINT  USING  "K’VTHE  MINIMUM  DETECTABLE  SIGNAL  WAS  SET  AT" .Noise , "WATTS 
I 

PRINT  “  " 

! 

! 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

OUTPUT  KBD«"K"; 

I 

GRAPHICS  ON 
VIEWPORT  15. 120. 10,70 
! 

! 

(PLOT  RANGE  VERSUS  ANTENNA  DIAMETER  FOR  VARIOUS  FREQUENCIES 

Xnin-0 

Xnax«6 

Ynin«PROUND< Ynin-50,2 ) 

Ynax -PROUND1 Ynax +50 , Z ) 

Dx-1 
Dy- 1 00 

Xr  ange*  Xnax  -  Xn  i  n 
Yrange"Ynax-Ynin 
WINDOW  Xnin.Xnax .Ynin.Ynax 
AXES  Dx ,Dy .Xnin.Ynin, 1 , 1 
CLIP  OFF 
! 

! LABEL  PLOT 


A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA 
AS  PRESENTED  IN  WALTERS/ KUNDEL  1981  PAGE  3S8 
WHICH  GIVES  A  HEIGHT  TO  THE  -1.1G  PROFILE" 

A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA 
AS  PRESENTED  IN  WALTERS/ KUNDEL  1981  PAGE  398 
WHICH  GIVES  A  HEI6HT  TO  THE  -1.16  PROFILE" 

BUT  WITH  A  FACTOR  OF  TWO  TO  APPROXIMATE" 
LOOKING  UP  A  THERMAL  PLOOM" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA 
AS  PRESENTED  IN  WALTERS/ KUNDEL  1981  PAGE  398" 
FOR  NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH" 
HEIGHT  IS  EXP( - . 001 + HEIGHT  ABOVE  G5  METERS" 
AND  HEIGHT  TO  THE  -1.4G  BELOW  G5  METERS" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 
HEI6HT  TO  THE  -4/3  " 
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CSIZE  4..S 
LOIR  © 

LORG  5 

MOVE  . 5*Xrange,Ynax*1 .0S 

LABEL  "RANGE  VERSUS  ANTENNA  DIAMETER" 

I  LABEL  HORIZONTAL  AXES 
LDIR  0 
LORG  5 

FOR  J-0  TO  Xnax  STEP  Dx 
CSIZE  4, .6 

MOVE  J , Ynin- .05*Yrange 
LABEL  J 
NEXT  J 

MOVE  . 5*Xrange,Ynin-. 1*Yrange 
CSIZE  4.. 6 

LABEL  "ANTENNA  DIAMETER  <n.  >" 

! LABEL  VERTICAL  AXES 
LORG  8 

FOR  J“Ynin  TO  Ymax  STEP  Dy 
CSIZE  4 , . 6 

MOVE  Xmin- .00Z5«Xrange , J 
LABEL  USING  "K";J 
NEXT  J 
LDIR  PI/Z 
LORG  6 

MOVE  Xnin- . 1 S*Xrange , . S« Yrange+Ynin 
CSIZE  4, .6 

LABEL  "RANGE  (METERS)" 

CLIP  ON 

Ind-1 

Ifraq-1 

FOR  Freq-500  TO  1500  STEP  500 

FOR  Ant_dian-.5  TO  6.Z5-Ifreq  STEP  .Z5 
IF  INT< Ifreq/Z  >«Ifraq/Z  THEN 
LINE  TYPE  3 
ELSE 

LINE  TYPE  1 
END  IF 

PLOT  Ant_dian, Ran( Ind.Ifreq) 
Ind»Ind+1 
NEXT  Ant_dian 
LOIR  0 
L0R6  Z 
CSIZE  3,.G 

MOVE  An-l_dian-  .  1  ,  Ran<  Ind-  1  , 1  f  req  > 

LINE  TYPE  1 

LABEL  USING  "K" jFreq,"Hz" 

Ifreq-Ifreq+1 
Ind»1 
NEXT  Freq 


PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT  KBO; "K"; 

1 

! 

LINPUT  " UOULD  YOU  LIKE  TO  MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?",X* 
Again-FNYes(  X$> 

Mess_up- 1 
WHILE  Mess_up-I 
SELECT  Again 
CASE  1 

LINPUT  "WOULD  YOU  LIKE  TO  CHANGE  A  VARIABLE  (ENTER  Y  OR  N)?" ,XS 
Neu_va»FNYes<  X$> 

Ano_change- 1 
WHILE  Ano_ change- 1 
SELECT  Neu.va 
CASE  1 

PRINT  "  VARIABLE  CURRENT  VALUE" 

PRINT  USING  “K” ; " 1  TEMPERATURE  "1Te»p,“ CELSIUS" 

PRINT  USING  "K“ s "2  ATOMOSPHERIC  PRESSURE  " . Aton_pres , "nb‘ 
PRINT  USING  "K"j"3  WATER  VAPOR  PRESSURE  " . Wat_pres , "nb" 
PRINT  "  ” 

PRINT  USING  “K" s ”4  ANTENNA  DIAMETER  ",Ant_dian,“  n." 

PRINT  USING  "K" j "5  PULSE  LEN6TH  ".Pulse,"  ns" 

PRINT  USING  "K" ; “6  POWER  TRANSMITTED  “ ,Pou_trans WATTS 
PRINT  USIN6  "K'V'7  ATMOSPHERIC  PROFILE  ", Profile 
PRINT  "  " 

PRINT  "ENTER  THE  NUMBER  OF  THE  VARIABLE  YOU  WISH  TO” 

PRINT  "CHANGE" 

INPUT  Var 
SELECT  Var 

CASE  1 

IMPUT  "TEMPERATURE-", Temp 
CASE  Z 

INPUT  "ATMOSPHERIC  PRESSURE  IN  nb-" ,Atom_pres 
CASE  3 

INPUT  "WATER  VAPOR  PRESSURE  IN  nb-" ,Wat_pres 
CASE  4 

INPUT  "ANTENNA  DIAMETER  IN  n.-" ,Ant_dian 
CASE  5 

INPUT  "PULSE  LENGTH  IN  ns-", Pulse 
CASE  G 

INPUT  "POWER  TRANSMITTED  IN  WATTS-" ,Pou_trans 
CASE  7 

PRINT  "  1  — >  Ct"Z  PROFILE  OF  Z*(-1.16>" 

PRINT  "  FROM  WALTERS/KUNOEL  1981" 

PRINT  "  Z  — >  SAME  AS  ONE  BUT  WITH  FACTOR- 
PRINT  "  OF  TWO  TO  APPROXIMATE  LOOKING" 


PRINT  “  UP  ft  THERMAL  PLOOM" 

PRINT  *'  3  — >  Ct*Z  PROFILE  OF  EXP( -.0B1*Z>“ 

PRINT  "  FROM  UflLTERS/ K UNOEL  1981." 

PRINT  “  4  -->  CT-Z  PROFILE  OF  Z*(-4/3>“ 

INPUT  "ENTER  NUMBER  OF  DESIRED  PROFILE” .Prof ile 
IF  Profile-4  THEN 

INPUT  “HEIGHT  OF  INVERSION  IN  METERS-” . Inver 
END  IF 
CASE  ELSE 

PRINT  Var.”IS  NOT  ONE  OF  THE  OPTIONS" 

END  SELECT 

LINPUT  "MADE  ANOTHER  CHANGE< ENTER  Y  OR  N>?",X$ 
Neu_va-FNYes( XS ) 

Mees_up-Z 
CASE  Z 

Mese_up-Z 
Ano_chanQe-Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  WHERE  Y  OR  N  ! !!! ! ” 

Ano_chanQe-1 

LINPUT  "CHANGE  A  VARIABLE  (ENTER  Y  OR  N)?”,XS 
Neu_va-FNYee(X$> 

END  SELECT 
END  WHILE 
CASE  Z 

Mesa_up-Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  WHERE  Y  OR  N  !!!!!!" 

Meaa_up**1 

LINPUT  "MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?",XS 
AQain-FNYee<  X*) 

END  SELECT 
END  WHILE 
END  WHILE 

PRINT  "THAT’S  ALL.  FOLKS" 

end 

I 


(CALCULTE  THE  ATTENUATION 

sub  At t enuat i on<  At on_pr es , At t en , Fr eq , T enp , Wat_pree  > 

! THIS  SUBPROGRAM  CALCULATES  THE  ATTENUATION  OF  SOUND 
! IN  AIR  BASED  UPON  EOUATIONS  IN  NEFF  1975 

I 

!  INPUT  ATMOSPHERIC  PRESSURE  IN  MILLIBARS 
!  FREQUENCY  OF  SOUND  WAVE  IN  HERTZ 

!  TEMPERATURE  IN  DEGREES  CELCIUS 

!  WATER-VAPOR  PRESSURE  IN  MILLIBARS 

1 

!  OUTPUT  ATTENUATION  IN  1 /METERS 

I 
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!  variables 

!  Atom _prea  INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 

I  Atten  ATTENUATION  OF  ACOUSTIC  UAVE.  CALCULATED  IN 

!  SUBPROGRAM  ATTENUATION. 

!  Att_max  VARIABLE  IN  SUBPROGRAM  ATTENUATION.  IT  IS  THE 
!  ~  ATTENUATION  AT  THE  FREQUENCY  OF  THE  MAXIMUM 

!  ATTENUATION  FOR  THE  INPUT  CONDITIONS. 

!  F  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.  IS  THE 

I  RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 

!  Fmax  FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 

!  ATTENUATION. 

!  Freq  INPUT  FREQUENCY  OF  ECHOSOUNDER. 

I  H  VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

!  Pstar  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 

!  CALCULATION. 

I  Temp  INPUT  TEMPERATURE  IN  DEGREES  CELSIUS, 

j  Tatar  INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

!  IN  SUBPROGRAM  ATTENUATION. 

|  Uat_prea  ATMOSPHERIC  UATER  PRESSURE  IN  MILLIBARS.  INPUT  BY 

I  ~  OPERATOR. 

H” I 00*Uat_prea/Atom_prea 
Tatar«<  1 . 8*Temp+49Z )/519 
Patar-Atom_prea/ 1014 

Fmax-<  10+6600*H+44400»H»H>*Pstar/TstarA .  B 

Att_max« .0078* Fmax* Tatar *< -Z . 5 )*EXP<  7. 77»( 1  -  1 /Tstar ) ) 

F«Freq/Fmax 

At ten«( Att_max/304. 8 ) • < ( . 1 8»F ) * Z  +  < 2*F*F/( 1 +F»F )  ) AZ >* . 5 
Atten«<  Atten+  t .  74E- 10»Freq*Freq)/4. 35 
SUBEND 
l 

def  FNYes(XS) 

! THIS  FUNCTION  INTERPRETS  THE  OPERATORS  RESPONSES  TO  YES  NO  QUESTIONS 
! 

!  INPUT  XS 

! 

!  OUTPUT  FNYea 

(VARIABLES 

!  Temp*  VARIABLE  STRING  USED  IN  FUNCTION  YES. 

!  X*  STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

t  RESPONSE  TO  YES  OR  NO  QUESTION. 

DIM  TempSl I ] 

TempSl t , t]«TRIMS<  XS) 

SELECT  Temp* 

CASE  "Y","y" 

RETURN  I 
CASE  "N" , "n" 

RETURN  Z 
CASE  "  " 

RETURN  1 
CASE  ELSE 


APPENDIX  E 

COMPUTER  PROGRAM  5 

I  FULLER,  ROBERT  PR0G_5  POWER 

I 

1 10  SEP  85 


.♦.••••PURPOSE******* 

THIS  PROGRAM  WILL  TAKE  PARAMETERS  OF  AN  ECHOSOUNOER  AND 
ESTIMATE  THE  RANGE  AS  A  FUNCTION  OF  POWER  TRANSMITTED  BY  THE  ECHOSOUNDER 
FOR  A  RANGE  OF  FREQUENCIES 
THE  FOLLOWING  INPUTS  ARE  REQUIRED: 

ATMOSPHERIC  DATA 

1  )ATOMOSPHERIC  PRESSURE  in  millibar* 

Z ) SELECT  ONE  OF  FOUR  ATMOSPHERIC  PROFILES  FOR  THE  TEMPERATURE 
STRUCTURE  PROFILE. 

a)FOR  ONE  PROFILE  THE  HEIGHT  OF  THE  INVERSION  LAYER  IS  ALSO  INPUT 

3 )  TEMPERATURE  IN  DEGREES  CELSIUS 

4)  WATER  VAPOR  PRESSURE  IN  millibar* 

ECHOSOUNOER  DATA 

5)  FREQUENCY 

G) ANTENNA  DIAMETER  IN  METERS 

71PULSE  LENGTH  OF  THE  TRANSMITTED  ACOUSTIC  ENERGY 

THE  PROGRAM  OUTPUTS  THE  FOLLOWING  GRAPH  TO  AN  EXTERNAL  PLOTTER: 

1 )RAN6E( m. )  VERSUS  POWER  TRANSMITTED  IN  WATTS  FOR  VARIOUS  FREQUENCIES 


VARIABLES 


Again  TEST  VALUE  IN  WHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  MAKE  ANOUTHER  RUN. 

Ano_change  TEST  VALUE  IN  WHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  CHANGE  THE  VALUE  OF  INPUT  DATA  BEFORE  MAKING 
ANOUTHER  RUN. 

Ant.araa  ECHOSOUNDER  ANTENNA  AREA  CALCULATED  FROM  INPUT  OF 
ANTENNA  DIAMETER. 

Ant_diam  INPUT  OF  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS. 

Atom_prea  INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 

Atten  ATTENUATION  OF  ACOUSTIC  WAVE.  CALCULATED  IN 

SUBPROGRAM  ATTENUATION. 

Bn  BANDWIDTH  OF  RECEIVER  NEEDED  TO  RECEIVE  SIGNALS 

DOPPLER  SHIFFED  BY  THREE  METER  PER  SECOND 
VERTICAL  VELOCITIES. 

C  VELOCITY  OF  SOUND  CALCULATION  AT  ZERO  CELSIUS. 

Cl  VARIABLE  USED  TO  SCALE  INVERSION  HEI6HT  IN 

TEMPERATURE  STRUCTURE  PROFILE  THREE. 

C3  VARIABLE  USED  IN  CALCULATING  THE  FOURTH  Cte2  PROFILE. 

CneZU)  ARRAY  OF  VALUES  OF  THE  ACOUSTIC  REFRACTIVE  INDEX 


I  CteZ 


I  CveZ 


I  Ox 
I  Dy 

!  Epsilon 
!  Eff 
I  Eff 
! 

!  Exc_att 
I  F 
I 

!  Fnax 


!  Ifreq 
!  IND 
!  Interval 
!  Inver 
!  J 
I  K 
!  L 

!  Hesa_up 


!  Neu_va 

! 

!  Noise 
!  Pou_back  < • 
!  Pou_ret(*) 
!  Pou_trans 
!  Profile 
!  Pstar 


Pulse 


!  R 
!  R1 

!  Range( • ) 


PARAMETER.  CALCULATED  BASED  ON  SELECTION  OF 
PROFILE  FOR  TEMPERATURE  STRUCTURE  PARAMETER 
AND  ASSUMED  PROFILE  FOR  THE  VELOCITY  STRUCTURE 
PARAMETER. 

VALUE  OF  THE  TEMPERATURE  STRUCTURE  PARAMETER. 

THE  OPERATOR  SELECTS  A  PROFILE  FROM  SEVERAL  PROVIDED. 

VALUE  OF  THE  VELOCITY  STRUCTURE  PARAMETER. 

VALUES  BASED  ON  CALCULATION  USING  ASSUMED  OISSAPATION 
RATE. 

STEP  SIZE  FOR  X  AXIS  FOR  VARIOUS  PLOTS. 

STEP  SIZE  FOR  Y  AXIS  FOR  VARIOUS  PLOTS. 

DISSAPATION  RATE  USED  IN  CALCULATION  OF  CveZ. 

TRANSMISSION  EFFICIENCY  OF  ECHOSOUNDER. 

EFFICIENCY  OF  ECHOSOUNDER  WHEN  RECEIVING  ACOUSTIC 
BACK SCATTER. 

EXCESS  "ATTENUATION"  AT  GIVEN  RANGE. 

VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.  IS  THE 
RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 

FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 
ATTENUATION. 

INPUT  FREQUENCY  OF  ECHOSOUNDER. 

ANTENNA  EFFECTIVE  APERATURE  FACTOR. 

VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

MAIN  INDEX  FOR  LOOP  TO  CHECK  RETURNED  POUER  UNTIL 
POWER  WAS  LESS  THAN  MINIMUM  DETECTABLE  SIGNAL. 

INDEX  FOR  INCREMENTS  OF  FREQUENCY 

INDEX  USED  FOR  LOOP  FOR  DIFFERENT  TRANSMITTED  POUER 

MAXIMUM  EFFECTIVE-SCATTERING  VOLUME  THICKNESS. 

HEIGHT  OF  INVERSION  LAYER. 

FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

UAVENUMBER 

SECOND  ORDER  INOEX  FOR  ASSORTED  LOOPS. 

TEST  VALUE  TO  CHECK  RESPONSES  WHEN  OPERATOR  IS  INPUTTING 
RESPONSES. 

TEST  VALUE  IN  CALCULATION  OF  COHERENCE  LENTH  IN  CALCULATION 
OF  EXCESS  ATTENUATION. 

TEST  VALUE  IN  CHECKING  IF  OPERATOR  WISHES  TO  CHANGE 
A  VARIABLE  BEFORE  A  NEW  RUN. 

ASSUMEO  MINIMUM  DETECTABLE  SIGNAL. 

)  POUER  BACKSCATTERED  FROM  6IVEN  RANGE. 

POWER  BACKSCATTERED  TO  ECHOSOUNDER  FROM  GIVEN  RANGE. 

INPUT  OF  POWER  SUPPLIED  TO  TRANSDUCER  OF  ECHOSOUNDER. 

OPERATOR  INPUT  OF  CteZ  PROFILE  FROM  AVAILABLE  PROFILES. 

VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 
CALCULATION. 

TRANSMITTED  PULSE  LENGTH  IN  MILLISECONDS.  INPUT  BY 
OPERATOR. 

RANGE  FOR  INPUT  TO  INTEGRAL  OF  EXCESS  ATTENUATION. 

VARIABLE  USED  IN  SCALING  INVERSION  HEIGHT  FOR  PROFILE  THREE. 

ARRAY  OF  RANGE  VALUES. 
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!  Remainder  REMAINDER  OF  MODULO  FUNCTION  USED  TO  DECREES  THE 

NUMBER  OF  PASSES  THR0U6H  THE  INTEGRATION  FOR  EXCESS 
ATTENUATION. 

Rge  VALUE  OF  RANGE  IN  INTESRAL  FOR  EXCESS  ATTENUATION. 

Rho  CORRELATION  LENGTH  USED  IN  CALCULATION  OF  EXCESS 

ATTENUATION. 

Sigma  FRACTION  OF  ENERGY  BACKSCATTERED  FROM  GIVEN  RANGE. 
Speed. sound  SPEED  OF  SOUND  IN  AIR  AT  INPUT  TEMPERATURE. 
Sumpou.back 

T  INPUT  TEMPERATURE  IN  DEGREES  KELVIN. 

Temp  INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

Temp*  VARIABLE  STRING  USED  IN  FUNCTION  YES. 


Tstar  INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

IN  SUBPROGRAM  ATTENUATION. 

Var  USED  TO  SELECT  THE  VARIABLE  THE  OPERATOR  UISHED  TO 

CHANGE  BEFORE  MAKING  ANOTHER  RUN. 

Uat_pres  ATMOSPHERIC  UATER  PRESSURE  IN  MILLIBARS.  INPUT  BY 
OPERATOR. 

X  THIRD  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 

X*  STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

RESPONSE  TO  YES  OR  NO  QUESTION. 

Xnax  VALUE  OF  LARGEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

Xmin  VALUE  OF  SMALLEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

Xrange  VALUE  OF  RANGE  OF  X  VALUES  FOR  VARIOUS  PLOTS. 

Ylabel*  LABEL  ON  Y  AXIS  PASSED  TO  SUBPROGRAM  Semi_log  FOR  PLOTTING. 
Ymax  LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

Ymin  SMALLEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

Yrange  RANGE  OF  Y  VARIABLES  FOR  VARIOUS  PLOTS. 

Ze  INTERMEDIATE  VALUE  IN  CALCULATION  OF  EXCESS  ATTENUATION. 


! DECLARE  VARIABLES 
INTEGER  I 

DIM  Pou_back < 1500) ,Pou_ret( 1500) ,Range( 1500) 

DIM  Cne2 < 1 500 ) . Ran(  50,5) 

! 

PLOTTER  IS  705 , " HP6L” 

LINE  TYPE  1 
I 

! INPUT  ATMOSPHERIC  DATA 

INPUT  "ENTER  TEMPERATURE  IN  DEGREES  CELCIUS" .Temp 
INPUT  "ENTER  ATMOSPHERIC  PRESSURE  IN  MILLIBARS" .Atom jsrea 
INPUT  "ENTER  WATER  VAPOR  PRESSURE  IN  MILLIBARS" , Wat.pree 


! INPUT  ECHOSOUNDER  DATA 

INPUT  "ENTER  ECHOSOUNDER  PULSE  LENGTH  IN  MILLISECONDS-" .Pulse 
INPUT  "ENTER  THE  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS" .Ant.di an 
G- .40  ! ANTENNA  EFFECTIVE  APERATURE  FACTOR 
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Eff-.Z5 

! 

I 

I  SELECT  THE  PROFILE  FOR  THE  TEMPERATURE  STRUCTURE 

PRINT  "YOU  HAVE  A  CHOICE  OF  TEMPERATURE  STRUCTURE  PROFILES" 

PRINT  "PROFILE  TO  BE  USEO." 

PRINT  "YOUR  SELECTIONS  ARE" 

PRINT  "  1  A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS  " 

PRINT  “  PRESENTED  IN  UALTERS/KUNDEL  1981  PA6E  398  WHICH  " 

PRINT  "  GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 

PRINT  “  " 

PRINT  “  Z  THE  TEMPERATURE  STRUCTURE  PROFILE  ABOVE  BUT" 

PRINT  "  WITH  A  FACTOR  OF  TWO  TO  APPROXIMATE  LOOKING" 

PRINT  "  UP  A  THERMAL  PLOOM" 

PRINT  "  " 

PRINT  "  3  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 

PRINT  “  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 

PRINT  "  NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH  HEIGHT  " 

PRINT  "  IS  EXP(-.001»HEIGHT  ABOVE  65  METERS  AND  HEIGHT" 

PRINT  "  TO  THE  -1.46  BELOU  65  METERS" 

PRINT  "  " 

PRINT  “  4  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 

PRINT  "  HEIGHT  TO  THE  -4/3  " 

! 

PRINT  M  " 

PRINT  "  " 

Meas_up- 1 
WHILE  Mea6_up-1 

INPUT  "ENTER  THE  DESIRED  PROFILE  <1  OR  Z  OR  3  OR  4)", Profile 
IF  Profile-4  THEN 

INPUT  "ENTER  THE  HEIGHT  OF  THE  INVERSION  IN  METERS" , Inver 
END  IF 

IF  Profile-1  THEN 
Meaa_up-0 

ELSE 

IF  Profile-Z  THEN 
Mess_up-0 

ELSE 

IF  Profile- 3  THEN 
Mess_up-0 

ELSE 

IF  Profile-4  THEN 
Mes6_up-0 

ELSE 

PRINT  Profile,"  WAS  NOT  ONE  OF  THE  OPTIONS!!!!" 
Mess_up- ! 

END  IF 
END  IF 
ENO  IF 
END  IF 
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END  WHILE 
OUTPUT  KBDj "K" i 
! 
t 

Again" 1 
WHILE  Again- t 
FOR  J-1  TO  I 
Pow_back  ( J ) -0 
Pow_ret< J)-0 
Range ( J)-0 
Cne2< J>-0 
NEXT  J 
! 

! CONVERT  TEMPERATURE  TO  KELVIN 

T-Tenp+Z73 

I 

! CALCULATE  SPEED  OF  SOUND  BASED  ON  INPUT  TEMPERATURE< CELCIUS ) 
Speed_aound-Z0.05«(  T)“ .5 

(CALCULATE  THE  SPEED  OF  SOUND  AT  0  DEGREES  CELCIUS 
C-Z0. 05*273*. 5 
! 

I-l  I  INDEX  FOR  INNER  LOOP  (POWER  RETURNED  TO  NOISE) 

Ind-I  (INDEX  FOR  SECOND  LOOP  (EFFICIENCY) 

Ifreq-1  (INDEX  FOR  THIRD  LOOP  (FREQUENCY) 

Pou_back(0)-0  (INITIALIZE  VARIABLE  FOR  POWER  BACKSCATERRED 
Intervals  Speed_aound*Pulae« 1 .E-3)/Z 
Ant_area-PI»<  Ant_dia«/2  )*Z 
! 

(INITIALIZE  PLOTTING  PARAMETERS  SO  CAN  COMPARE  DURING  RUN  TO  SCALE 

(PLOTS 

Ynin-1000 

Ynax-0 

(CALCULATE  RETURNED  SIGNAL  POWER  UNTIL  IT  IS  LESS  THAN  BACKGROUND 
IFOR  EACH  FREQUENCY. 

Range( 0)-0 

FOR  Freq-500  TO  1500  STEP  500 
FOR  Pou_trana-50  TO  500  STEP  50 

Bn-Z*Freq*( 1  - 1 /( 6/Speed_aound+1 >) (BANDWIDTH  FOR  3M/S  VERTICAL  VELOCITY 
Noiee-1 . 38E-Z3*Bn*< T )+Z . E- 14! MINIMUM  DETECTABLE  SIGNAL  AS 
!  JOHNSON  NOISE( NEGLIGABLE  >4  ESTIMATED  BACKGROUND 
(CALCULATE  THE  ATTENUATION  COEFFICIENT  BASED  ON  ATMOSPHERIC  DATA 
CALL  At t enuat i on( At on_pr ea , At  t  en , Freq , T enp , Wat_prea  > 
K-2*PI»Freq/Speed_aound  ! WAVENUMBER 
I-l 

Sunpou_back -0 
REPEAT' 

Range( I )-Range(  I- I )+2 
SELECT  Profile 
CASE  I 

I  THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
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IAS  PRESENTED  IN  WALTERS/ KUNDEL  1981  PAGE  398. 

I  WAS  FOR  MIDOAY  AND  CLEAR  WEATHER  ABOVE  THE 
! TULAROSA  BASIN,  NEW  MEXICO. 

IAN  ADDITONAL  FACTOR  OF  1.8  WAS  INCLUDED  AS  PER  DIFFERENT 
I  AVERAGING  TIME. 

CteZ-Z. 1Z#Range<I>*<-1. 16) 

CASE  Z 

I  THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
IAS  PRESENTED  IN  WALTERS/KUNOEL  1981  PAGE  398. 

I  WAS  FOR  MIDOAY  AND  CLEAR  WEATHER  ABOVE  THE 
I TULAROSA  BASIN,  NEW  MEXICO. 

IAN  ADDITONAL  FACTOR  OF  1 . 8  WAS  INCLUDED  AS  PER  DIFFERENT 
I  AVERAGING  TIME  AND  AN  ADDITIONAL  FACTOR  OF  TWO  WAS 
I  INCLUDED  TO  APPROXIMATE  LOOKING  UP  A  THERMAL  PLOOM 
CteZ*Z»Z. 1Z*Range< I >“( -  1 . IB) 

CASE  3 

I  THIS  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 

(PRESENTED  IN  WALTERS/KUNDEL  1981  PA6E  398  FOR" 

I  NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH  HEIGHT  " 

IIS  EXP( - .001 «HEI6HT  ABOVE  GS  METERS  AND  HEIGHT" 

I  TO  THE  -1.4G  BELOW  65  METERS" 

IF  Rangel  I ) <65  THEN 

Cte2«75.5«Range< I )*(-! .4S> 

ELSE 

CteZ“3.GEE-2*EXP< -.0O1*Range<  I ) > 

END  IF 
CASE  4 

I  CALCULATE  THE  TEMPERATURE  STRUCTURE  FACTOR  BASED  ON  ASSUMPTION 
I  THAT  IT'S  PROPORTIONAL  TO  HEIGHT  TO  THE  -4/3  POWER 

I  THAT  THE  SURFACE  VERTICAL  HEAT  FLUX  IS  .095  C.w./eec 

I  THE  DEPENDENCE  WITH  HEIGHT  ABOVE  THE  INVERSION  I  PULLED  FROM 
I  MY  LEFT  EAR 

I EOUATION  FROM  NEFF,  1975 
SELECT  Rangel  I ) / 1 nver 
CASE  <.9 

C3-<<  024 ) • ( T ) *  < . GG7  >  > 

CteZ-C3«< Ranged  >)*<-!.  33) 

CASE  .9  TO  1 
CteZ-CteZ 
Cl-CteZ 
R1 “Rangel I ) 

CASE  1  TO  1.3 

Cte2“10*( < LGTCC3*Range( I ) )-LGT(C3»R1 ) >+LGT(CI  ) ) 

RZ “Rangel I ) 

CZ-CteZ 
CASE  ELSE 

C t eZ“C 3* <  Rangel  I>'>(-1.33)-RZ“(-I.33))+CZ 
END  SELECT 
END  SELECT 

!  CALCULATE  VELOCITY  STRUCTURE  FACTOR.  FORMULA  FROM  GAYNOR  77 


Epsilon-< . Z866/Range( I >  )•< 1+.07*< Range< I ) )*< .6)>*< 1.5) 
Cve2*Z*Ep5ilon“( .BB7> 

I  CALCULATE  ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  FACTOR, 

(FORMULA  FROM  TETARSKI 

CneZ( I )■<  CteZ/( Z. 98E+5) )+< CveZ/< C»C > > 

(CALCULATE  THE  FRACTION  OF  POWER  BACKSCATTEREO  FROM  INTERVAL. 
(EQUATION  FROM  NEFF,  1975 
Signa-( ,0039*< K* < 1 /3) )*CteZ >/< T)*Z 
i 

(CALCULATE  THE  EXCESS  ATTENUATION.  EQUATION  FROM  CLIFFORD,  1980 
(THE  EXCESS  ATTENUATION  IS  Exc.att 
IF  I <50  THEN 

Reminder  *0 

ELSE 

IF  I < 1 50  THEN 

Remainder-Range<  I )  MODULO  10 

ELSE 

Renainder-Range(  I )  MODULO  2® 

END  IF 
END  IF 

IF  Remainder-®  THEN 
Rho-0 
L-0 
H-0 

Rge-Range< I >  I  CONSTANT  IN  INTEGRAL 
R-0 

FOR  J-0  TO  2*1 

F-CneZ< INK  J/Z+1 ) ) 

F-F*( 1 -R/<  Rge  )*<  1 . B7)  +  (  R/(  Rge ) ) * ( 1.B7>> 

IF  J>0  THEN 

IF  J<Z*I  THEN 

IF  INK  J/Z  >-J/Z  THEN 
L-L+F 
F-0 

ELSE 

H-H+F 
F-0 
END  IF 
END  IF 
END  IF 
Rho-Rho+F 
R-R+1 

NEXT  J 

Rho-RHo+4*L+Z*H 

Rho-< ( ( Rho*. 33  )*«•«* 1 .4B)*(  - . B) ) 

N-( Ant_diam/Rho ) " Z 
IF  N<-1  THEN 

Ze-1/< 1+N) 

ELSE 

Ze-1 .5/( 1+N) 
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(STEP  OF  1.5—  >SEE  CLIFFORD  1980 
END  IF 

Exc_att-Ze*Ze 

ELSE 

Exc_att-Exc_att 
END  IF 
I 

! CALCULATE  THE  POWER  BACKSCATTERED 

Pou_back( I )-<Pou_trana*Eff-Suiiipou_back >«EXP(  -At ten* Rangel  I > > 
Pou_back( I )“Pou_back  < I )*Interval*Exc_att*5ig«a 
Su«pou_back-Sumpou_back4Pou_back(  I ) 

(CALCULATE  THE  POWER  BACKSCATTERED  TO  THE  ANTENNA 
Pou_ret< I )-Pou_back( I )*EXP< -At ten* Rangel  I ) )*Ant_area*G*Ef  f /Range( I ) A2 
PRINT  "RANGE-". Ranged) 

PRINT  "POWER  RETURNED-"  ,Pou_retd> 

PRINT  "FREQUENCY-", Freq 
I-I  +  1 

UNTIL  Pou_ret(I-1  )>0  AND  Pou_ret(  1-1 ) <-Noi.se 
Ran< Ind.Ifreq)-Range( 1-2 ) 

IF  Randnd,Ifreq)>Ymax  THEN 
Ynax-Ran< Ind.Ifreq) 

END  IF 

IF  Ran< I nd , I f req ) < Ynin  THEN 
Ymin-Ran<  Ind.Ifreq) 

END  IF 
Ind-Ind+ I 
Again- 1 
FOR  J-0  TO  I 

Pou_back<  J)-0 
Pou_ret( J)-0 
Range < J)-0 
Cne2( J)-0 
NEXT  J 

NEXT  Pou_trana 
Ifreq-Ifreq+1 
Ind-1 
NEXT  Freq 

! 

I 

OUTPUT  KBD; "K" i 

PRINT  "INPUT  CONDITIONS" 

PRINT  M  " 

PRINT  USING  "K"; "TEMPERATURE-  " , Temp, "CELCIUS" 

PRINT  USING  "K" 5 "ATMOSPHERIC  PRESSURE-  " , Aton.pres . "nb" 

PRINT  USING  "K"s "WATER  VAPOR  PRESSURE-  " ,Wat_prea . "mb" 

PRINT  USING  "K"; "PULSE  LENGTH-  " .Pulae, "na” 

PRINT  USING  "K"j” TRANSDUCER  EFFICIENCY-  ”,Eff 

PRINT  USING  "K"; "POWER  TRANSMITTED-  " ,Pou_trana . "  WATTS" 

PRINT  "TEMPERATURE  STRUCTURE  PROFILE  USED  ".Profile 
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IF  Profile-4  THEN 

PRINT  “INVERSION  HEI6HT  ".Inver 
END  IF 

SELECT  Profile 
CASE  1 

PRINT  "  1  A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA 

PRINT  "  AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398 

PRINT  “  WHICH  GIVES  A  HEIGHT  TO  THE  -1.IE  PROFILE" 

CASE  Z 

PRINT  M  2  A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA 

PRINT  "  AS  PRESENTED  IN  WA.TERS/KUNDEL  1981  PAGE  398 

PRINT  "  UHICH  GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 

PRINT  "  BUT  WITH  A  FACTOR  OF  TWO  TO  APPROXIMATE" 

PRINT  "  LOOKING  UP  A  THERMAL  PLOOM" 

CASE  3 

PRINT  "  3  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA 

PRINT  "  AS  PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398" 

PRINT  "  FOR  NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH" 

PRINT  "  HEIGHT  IS  EXP( 00! * HEIGHT  ABOVE  G5  METERS" 

PRINT  "  AND  HEIGHT  TO  THE  -1.46  BELOU  65  METERS" 

CASE  4 

PRINT  "  4  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 

PRINT  "  HEIGHT  TO  THE  -4/3  " 

END  SELECT 

I 

PRINT  USING  "K" s "THE  MINIMUM  DETECTABLE  SIGNAL  WAS  SET  AT" .Noise . "UATTS 
! 

PRINT  "  " 

! 

! 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

OUTPUT  KBD; "K" ; 

I 

GRAPHICS  ON 
VIEWPORT  15,120,10,7© 


! PLOT  RANGE  VERSUS  POWER  TRANSMITTED  FOR  VARIOUS  FREQUENCIES 

Xnin-0 

Xmax-600 

Yniin*PROUND(  Ynin-50,2 ) 

Ymax-PROUND<  Ynax+50,2 ) 

Dx-100 
Dy- 100 

Xrange-Xnax-Xnin 
Yrange-Ymax-Ynin 
WINDOW  Xnin.Xnax .Ynin.Ymax 
AXES  Dx .Oy.Xmin, Ynin, 1 , 1 
CLIP  OFF 


! 

I  LABEL  PLOT 
CSIZE  4..S 
LDIR  0 
LOR6  5 

MOVE  . 5*Xrange,Ynax*1 .05 

LABEL  "RANGE  VERSUS  POUER  TRANSMITTED" 

t LABEL  HORIZONTAL  AXES 

LDIR  0 

LORG  5 

FOR  J-0  TO  Xnax  STEP  Dx 
CSIZE  4, .6 

MOVE  J, Ynin- .05*Yrange 
LABEL  J 
NEXT  J 

MOVE  . 5*Xrange, Ymin-. 1*Yrange 
CSIZE  4, .6 

LABEL  "POUER  TRANSMITTED  ( UATTS)" 

! LABEL  VERTICAL  AXES 
LORG  8 

FOR  J-Ynin  TO  Ymax  STEP  Dy 
CSIZE  4, .6 

MOVE  Xi<iin- .  00Z5*Xrange ,  J 
LABEL  USING  "K"sJ 
NEXT  J 
LDIR  PI/2 
LORG  G 

MOVE  Xnin-. 15*Xrange, . 5* Yrange+Ymin 
CSIZE  4 , .6 

LABEL  "RANGE  (METERS)" 

CLIP  ON 

Ind“t 

Ifreq"! 

FOR  Freq-500  TO  1500  STEP  500 
FOR  Pou_trans-50  TO  500  STEP  50 
IF  INT( Ifreq/Z  >"Ifreq/Z  THEN 
LINE  TYPE  3 
ELSE 

LINE  TYPE  1 
END  IF 

PLOT  Pou_trans ,Ran< Ind, Ifreq) 
Ind«Ind+ 1 
NEXT  Pou_trans 
LDIR  0 
LORG  Z 
CSIZE  3..B 

MOVE  Pow_trans-50 ,Ran( Ind- 1 .Ifreq > 
LINE  TYPE  1 

LABEL  USING  MK" i Freq . "Hz" 
Ifreq-Ifreq+1 
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v'vA.W  R,*'. 


T 


k 


Ind-1 
NEXT  Freq 


1 


PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT  KBD; "K“ j 


! 


I 


LINPUT  "UOULD  YOU  LIKE  TO  MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?",X$ 
Again-FNYes( XS ) 

Mess_up- 1 
UHILE  Mes6_up-1 
SELECT  Again 
CASE  1 

LINPUT  "UOULD  YOU  LIKE  TO  CHAN6E  A  VARIABLE  (ENTER  Y  OR  N>?",X* 
Neu_va«FNYes(X$> 

Ano_chang«“ 1 
UHILE  Ano_change-1 
SELECT  Neu_va 
CASE  1 

PRINT  “  VARIABLE  CURRENT  VALUE" 

PRINT 
PRINT 
PRINT 
PRINT 


USING 

"K“ 

J"1 

TEMPERATURE 

Temp. "CELSIUS" 

USING 

“K" 

i  "  2 

ATOMOSPHERIC  PRESSURE 

" .Atom _pres , "mb” 

USING 

M  M 

"K" 

!"3 

WATER  VAPOR  PRESSURE 

",Uat_pre8,"mb" 

USING 

"K" 

;  “4 

ANTENNA  DIAMETER 

" , Ant_diam,"  n." 

USIN6 

"K" 

s"S 

PULSE  LENGTH 

".Pulse,"  ns" 

USING 

"K" 

!"G 

POWER  TRANSMITTED  ",Pou  trans,"  WATTS" 

USING 

"K" 

5  "7 

ATMOSPHERIC  PROFILE 

" .Prof ile 

PRINT  "  " 

PRINT  "ENTER  THE  NUMBER  OF  THE  VARIABLE  YOU  UISH  TO" 

PRINT  "CHAN6E" 

INPUT  Var 
SELECT  Var 

CASE  1 

INPUT  "TEMPERATURE*" ,Tenp 
CASE  2 

INPUT  "ATMOSPHERIC  PRESSURE  IN  mb-” . Atom.pres 
CASE  3 

INPUT  "WATER  VAPOR  PRESSURE  IN  mb-" ,Uat_pres 
CASE  4 

INPUT  "ANTENNA  DIAMETER  IN  m. -" , Ant.dian 
CASE  S 

INPUT  "PULSE  LENGTH  IN  ms-". Pulse 
CASE  E 

INPUT  "POWER  TRANSMITTED  IN  WATTS-" ,Pou_trans 
CASE  7 

PRINT  "  1  -->  Ct*Z  PROFILE  OF  Z*<-1.16>" 

PRINT  "  FROM  UALTERS/KUNOEL  1981" 
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PRINT  "  Z  — >  SAME  AS  ONE  BUT  UITH  FACTOR- 
PRINT  “  OF  TWO  TO  APPROXIMATE  LOOKING- 

PRINT  "  UP  A  THERMAL  PLOOM" 

PRINT  "  3  -- >  Ct*Z  PROFILE  OF  EXP(-.001«Z>" 

PRINT  “  FROM  UAL TERS / K UNOEL  1981." 

PRINT  “  4  -- >  CT*Z  PROFILE  OF  Z*<-4/3>“ 

INPUT  "ENTER  NUMBER  OF  DESIRED  PROFILE" .Prof ile 
IF  Profile-4  THEN 

INPUT  "HEIGHT  OF  INVERSION  IN  METERS-" . Inver 
ENO  IF 
CASE  ELSE 

PRINT  Var , "IS  NOT  ONE  OF  THE  OPTIONS" 

END  SELECT 

LINPUT  "MADE  ANOTHER  CHAN6E( ENTER  Y  OR  N)?",X$ 
Neu_va-FNYes<X*> 

Mess_up-Z 
CASE  Z 

Mess_up-Z 
Ano_chanoe-Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  UHERE  Y  OR  N  !!!  !  ! " 

Ano_change- 1 

LINPUT  "CHANGE  A  VARIABLE  <  ENTER  Y  OR  N)?",XS 
Neu_va«FNYes(X*> 

END  SELECT 
END  WHILE 
CASE  Z 

Mess_up-Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  UHERE  Y  OR  N  III!!!" 

Meae_up-t 

LINPUT  "MAKE  ANOTHER  RUN  <  ENTER  Y  OR  N)?",Xt 
Aoain-FNYee(Xt) 

END  SELECT 
END  UHILE 
END  UHILE 

PRINT  "THAT’S  ALL.  FOLKS" 
end 
i 
! 

ICALCULTE  THE  ATTENUATION 

sub  Attenuat  ion<  Ato«_pres ,  Aiten ,  Freq ,  Temp  ,Ua,t_pres ) 

! THIS  SUBPROGRAM  CALCULATES  THE  ATTENUATION  OF  SOUND 
! IN  AIR  BASED  UPON  EQUATIONS  IN  NEFF  1975 
! 

!  INPUT  ATMOSPHERIC  PRESSURE  IN  MILLIBARS 
I  FREQUENCY  OF  SOUND  WAVE  IN  HERTZ 

!  TEMPERATURE  IN  DEGREES  CELCIUS 

!  WATER- VAPOR  PRESSURE  IN  MILLIBARS 
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OUTPUT  ATTENUATION  IN  1 /METERS 
VARIABLES 

Atom _preB  INPUT  OF  ATMOSPHERIC  PRESSURE  IN  nb. 

Atten  ATTENUATION  OF  ACOUSTIC  WAVE.  CALCULATED  IN 

SUBPROGRAM  ATTENUATION. 

Att  nax  VARIABLE  IN  SUBPROGRAM  ATTENUATION.  IT  IS  THE 
ATTENUATION  AT  THE  FREQUENCY  OF  THE  MAXIMUM 
ATTENUATION  FOR  THE  INPUT  CONDITIONS. 

F  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.  IS  THE 

RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 
Fnax  FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 

ATTENUATION. 

Freq  INPUT  FREQUENCY  OF  ECHOSOUNDER. 

H  VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

Pstar  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 

CALCULATION. 

Temp  INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

Tstar  INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

IN  SUBPROGRAM  ATTENUATION. 

Uat.pres  ATMOSPHERIC  UATER  PRESSURE  IN  MILLIBARS.  INPUT  BY 
i  “  OPERATOR. 

H« 100*Uat_pres/Aton_pres 
Tstar*( 1 . 8»Tenp+49Z )/51 9 
Pstar»Aton_pres/  1014 

Fnax*(  10+6600*H+44400*H*H)*Pstar/TstarA .  8 

Att_nax“ .  0078* Fnax* Tatar*  <  -  2.5 )  *EXP<  7.77*(  1-1/Tstar)  ) 

F-Freq/Fnax 

Atten-<  Att_nax/304.8)*<  (  .  1 8*F  )  *?-f  (  2«F«F/<  I +F*F > ) AZ > *. 5 

Atten«<  Atten+1 .74E- 10*Freq*Freq)/4. 35 

SUBEND 

i 

def  FNYes(XS) 

(THIS  FUNCTION  INTERPRETS  THE  OPERATORS  RESPONSES  TO  YES  NO  QUESTIONS 
! 

!  INPUT  X$ 

! 

!  OUTPUT  FNYes 

! VARIABLES 

!  TenpS  VARIABLE  STRING  USED  IN  FUNCTION  YES. 

!  X*  STRING  PASSED  TO  FUNCTION  YE5  TO  DETERMINE  OPERATOR 

!  RESPONSE  TO  YES  OR  NO  QUESTION. 

DIM  TenpS! 1 ] 

TenpSl 1 , 11-TRIMS<XS) 

SELECT  TenpS 
CASE  "Y" . "y” 

RETURN  I 
CASE  •'N'V'n" 

RETURN  Z 
CASE  "  “ 
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APPENDIX  F 
COMPUTER  PROGRAM  6 


FULLER,  ROBERT  PROG_B  NOISE 

10  SEP  85 


!  ***#***PURPOSE*«***** 

!  THIS  PROGRAM  MILL  TAKE  PARAMETERS  OF  AN  ECHOSOUNDER  AND 
(ESTIMATE  THE  RANGE  AS  A  FUNCTION  OF  BACKGROUND  NOISE 
(THE  F0LL0UIN6  INPUTS  ARE  REQUIRED: 

!  ATMOSPHERIC  DATA 

I  1 )ATOMOSPHERIC  PRESSURE  in  millibars 

(  Z ) SELECT  ONE  OF  FOUR  ATMOSPHERIC  PROFILES  FOR  THE  TEMPERATURE 
I  STRUCTURE  PROFILE. 

(  a)FOR  ONE  PROFILE  THE  HEIGHT  OF  THE  INVERSION  LAYER  IS  ALSO  INPUT 

I  3 ) TEMPERATURE  IN  DEGREES  CELSIUS 

!  4) WATER  VAPOR  PRESSURE  IN  millibars 

I  ECHOSOUNDER  DATA 

I  S) FREQUENCY 

I  G ) ANTENNA  DIAMETER  IN  METERS 

I  71PULSE  LENGTH  OF  THE  TRANSMITTED  ACOUSTIC  ENERGY 
!  8) POWER  TRANSMITTED 

(THE  PROGRAM  OUTPUTS  THE  FOLLOWING  GRAPH  TO  AN  EXTERNAL  PLOTTER: 

!  1 (RANGE ( m. )  VERSUS  BACKGROUND  NOISE  IN  WATTS 


VARIABLES 


Again  TEST  VALUE  IN  WHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  MAKE  ANOUTHER  RUN. 

Ano.change  TEST  VALUE  IN  WHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  CHANGE  THE  VALUE  OF  INPUT  DATA  BEFORE  MAKIN6 
ANOUTHER  RUN. 

Ant_area  ECHOSOUNDER  ANTENNA  AREA  CALCULATED  FROM  INPUT  OF 
ANTENNA  DIAMETER. 

Ant.diam  INPUT  OF  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS. 

Atom_pres  INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 

At-ton  ATTENUATION  OF  ACOUSTIC  WAVE.  CALCULATED  IN 

SUBPROGRAM  ATTENUATION. 

C  VELOCITY  OF  SOUND  CALCULATION  AT  ZERO  CELSIUS. 

Cl  VARIABLE  USED  TQ  SCALE  INVERSION  HEIGHT  IN 

TEMPERATURE  STRUCTURE  PROFILE  THREE. 

C3  VARIABLE  USED  IN  CALCULATING  THE  FOURTH  CteZ  PROFILE. 

CneZ<  * )  ARRAY  OF  VALUES  OF  THE  ACOUSTIC  REFRACTIVE  INDEX 
PARAMETER.  CALCULATED  BASED  ON  SELECTION  OF 
PROFILE  FOR  TEMPERATURE  STRUCTURE  PARAMETER 
AND  ASSUMED  PROFILE  FOR  THE  VELOCITY  STRUCTURE 


Ox 

Oy 

Epsilon 

Eff 

Eff 

Exc_att 

F 


1NO 

Interval 

Inver 

J 

K 

L 

Mess_up 


Neu_va 

Noise 
Pou_back( * 
Pou_ret( • ) 
Pou_trans 
Prof  lie 
Pstar 

Pulse 

R 

R1 

Range ( • ) 
Remainder 


PARAMETER. 

VALUE  OF  THE  TEMPERATURE  STRUCTURE  PARAMETER. 

THE  OPERATOR  SELECTS  A  PROFILE  FROM  SEVERAL  PROVIDED. 

VALUE  OF  THE  VELOCITY  STRUCTURE  PARAMETER. 

VALUES  BASED  ON  CALCULATION  USIN6  ASSUMED  OISSAPATION 
RATE. 

STEP  SIZE  FOR  X  AXIS  FOR  VARIOUS  PLOTS. 

STEP  SIZE  FOR  Y  AXIS  FOR  VARIOUS  PLOTS. 

OISSAPATION  RATE  USED  IN  CALCULATION  OF  CveZ. 

TRANSMISSION  EFFICIENCY  OF  ECHOSOUNDER. 

EFFICIENCY  OF  ECHOSOUNDER  WHEN  RECEIVING  ACOUSTIC 
BACKSCATTER. 

EXCESS  "ATTENUATION"  AT  GIVEN  RAN6E. 

VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.  IS  THE 
RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 

FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 
ATTENUATION. 

INPUT  FREQUENCY  OF  ECHOSOUNDER. 

ANTENNA  EFFECTIVE  APERATURE  FACTOR. 

VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

MAIN  INDEX  FOR  LOOP  TO  CHECK  RETURNED  POWER  UNTIL 
POUER  WAS  LESS  THAN  MINIMUM  DETECTABLE  SI6NAL. 

INDEX  USED  FOR  LOOP  FOR  DIFFERENT  NOISE  CUT  OFF  LEVELS 

MAXIMUM  EFFECTIVE- SCATTERING  VOLUME  THICKNESS. 

HEIGHT  OF  INVERSION  LAYER. 

FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

WAVENUMBER 

SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

TEST  VALUE  TO  CHECK  RESPONSES  WHEN  OPERATOR  IS  INPUTTING 
RESPONSES. 

TEST  VALUE  IN  CALCULATION  OF  COHERENCE  LENTH  IN  CALCULATION 
OF  EXCESS  ATTENUATION. 

TEST  VALUE  IN  CHECKING  IF  OPERATOR  WISHES  TO  CHAN6E 
A  VARIABLE  BEFORE  A  NEW  RUN. 

ASSUMED  MINIMUM  DETECTABLE  SIGNAL. 

)  POWER  BACKSCATTERED  FROM  GIVEN  RANGE. 

POWER  BACKSCATTERED  TO  ECHOSOUNDER  FROM  GIVEN  RANGE. 

INPUT  OF  POWER  SUPPLIED  TO  TRANSDUCER  OF  ECHOSOUNDER. 

OPERATOR  INPUT  OF  CteZ  PROFILE  FROM  AVAILABLE  PROFILES. 

VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 
CALCULATION. 

TRANSMITTED  PULSE  LENGTH  IN  MILLISECONDS.  INPUT  BY 
OPERATOR. 

RANGE  FOR  INPUT  TO  INTEGRAL  OF  EXCESS  ATTENUATION. 

VARIABLE  USED  IN  SCALING  INVERSION  HEIGHT  FOR  PROFILE  THREE. 

ARRAY  OF  RANGE  VALUES. 

REMAINDER  OF  MODULO  FUNCTION  USED  TO  DECREES  THE 
NUMBER  OF  PASSES  THROUGH  THE  INTEGRATION  FOR  EXCESS 
ATTENUATION. 

VALUE  OF  RANGE  IN  INTEGRAL  FOR  EXCESS  ATTENUATION. 
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VA'N 


*  W  %  »  ’  V>  ■ 


Rho  CORRELATION  LENGTH  USED  IN  CALCULATION  OF  EXCESS 

ATTENUATION. 

Signs  FRACTION  OF  ENERGY  BACK SCATTERED  FROM  GIVEN  RANGE. 
Spsed.sound  SPEED  OF  SOUND  IN  AIR  AT  INPUT  TEMPERATURE. 

T  INPUT  TEMPERATURE  IN  DEGREES  KELVIN. 

Tenp  INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

Sunpou_back 

TsoipS  VARIABLE  STRING  USED  IN  FUNCTION  YES. 

Tatar  INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

IN  SUBPROGRAM  ATTENUATION. 

Var  USED  TO  SELECT  THE  VARIABLE  THE  OPERATOR  WISHED  TO 

CHANGE  BEFORE  MAKING  ANOTHER  RUN. 

Wat_pres  ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.  INPUT  BY 
OPERATOR. 

X  THIRD  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 

X*  STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

RESPONSE  TO  YES  OR  NO  QUESTION. 

Xnax  VALUE  OF  LARGEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

Xnin  VALUE  OF  SMALLEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

Xrange  VALUE  OF  RANGE  OF  X  VALUES  FOR  VARIOUS  PLOTS. 

YlabslC  LABEL  ON  Y  AXIS  PASSED  TO  SUBPROGRAM  Seni.log  FOR  PLOTTING. 
Ymax  LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

Yain  SMALLEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

Yrange  RANGE  OF  Y  VARIABLES  FOR  VARIOUS  PLOTS. 

Ze  INTERMEDIATE  VALUE  IN  CALCULATION  OF  EXCESS  ATTENUATION. 


! DECLARE  VARIABLES 
INTEGER  I 

DI M  Pou_back ( 5000 ) , Pou_r at  <  5000 ) , Rangex .  000 ) 

DIM  CneZ<  5000) ,Ran(  50) 

I 

PLOTTER  IS  705,"HPGL" 

LINE  TYPE  1 
! 

(INPUT  ATMOSPHERIC  DATA 

INPUT  "ENTER  TEMPERATURE  IN  DEGREES  CELCIUS" ,Te*p 
INPUT  "ENTER  ATMOSPHERIC  PRESSURE  IN  MILLIBARS" ,Aton_pres 
INPUT  "ENTER  WATER  VAPOR  PRESSURE  IN  MILLIBARS" ,Wat_pre» 


(INPUT  ECHOSOUNDER  DATA 

INPUT  "ENTER  ECHOSOUNDER  FREQUENCY  IN  HERTZ" ,Freq 
INPUT  "ENTER  ECHOSOUNDER  PULSE  LENGTH  IN  MILLISECONDS-" .Pulse 
INPUT  "ENTER  THE  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS" .Ant.dian 
INPUT  "ENTER  THE  POWER  TRANSMITTED  IN  WATTS" ,Pou_tran« 

6-. 40  (ANTENNA  EFFECTIVE  APERATURE  FACTOR 

Eff-.Z5 

i 
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f SELECT  THE  PROFILE  FOR  THE  TEMPERATURE  STRUCTURE 

PRINT  -YOU  HAVE  A  CHOICE  OF  TEMPERATURE  STRUCTURE  PROFILES" 

PRINT  “PROFILE  TO  BE  USEO." 

PRINT  “YOUR  SELECTIONS  ARE" 

PRINT  “  1  A  TEHPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS  " 

PRINT  ”  PRESENTED  IN  WALTERS/ KUNOEL  1981  PAGE  3S8  WHICH  “ 

PRINT  "  GIVES  A  HEIGHT  TO  THE  -1.1B  PROFILE" 

PRINT  "  “ 

PRINT  “  2  THE  TEMPERATURE  STRUCTURE  PROFILE  ABOVE  BUT" 

PRINT  "  WITH  A  FACTOR  OF  TWO  TO  APPROXIMATE  LOOKING" 

PRINT  “  UP  A  THERMAL  PLOOM“ 

PRINT  "  " 

PRINT  "  3  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 

PRINT  “  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 

PRINT  "  NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH  HEIGHT  " 

PRINT  "  IS  EXP< - .001 ‘HEIGHT  ABOVE  B5  METERS  AND  HEIGHT" 

PRINT  "  TO  THE  -1.4B  BELOW  B5  METERS" 

PRINT  "  " 

PRINT  "  4  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 

PRINT  "  HEIGHT  TO  THE  -4/3  " 

! 

PRINT  "  " 

PRINT  "  " 

Mes8_up“ 1 
WHILE  Meaa_up- 1 

INPUT  "ENTER  THE  DESIRED  PROFILE  (1  OR  Z  OR  3  OR  4)", Profile 
IF  Profile-4  THEN 

INPUT  "ENTER  THE  HEIGHT  OF  THE  INVERSION  IN  METERS" .Inver 
END  IF 

IF  Profile-1  THEN 
Meaa_up-0 

ELSE 

IF  Profile-Z  THEN 
Meaa_up-® 

ELSE 

IF  Profile-3  THEN 
Meaa_up-0 

ELSE 

IF  Profile-4  THEN 
Meaa_up-0 

ELSE 

PRINT  Profile,”  WAS  NOT  ONE  OF  THE  OPTIONS! It!* 
Meaa_up- 1 
ENO  IF 
END  IF 
END  IF 
END  IF 
ENO  WHILE 
OUTPUT  KBDj "K”  i 


I 

I 

Again*) 

WHILE  Again*! 

FOR  J-1  TO  I 
Pou.back ( J ) *0 
Pou_ret) J)-0 
Range ( J>-0 
CneZ) J>-0 
NEXT  J 
I 

! CONVERT  TEMPERATURE  TO  KELVIN 
T-Tenp+Z73 

j 

(CALCULATE  SPEED  OF  SOUND  BASED  ON  INPUT  TEMPERATURE) CELCIUS ) 
Speed_sound*Z0. 05* < T ) * . 5 

(CALCULATE  THE  SPEED  OF  SOUND  AT  0  DEGREES  CELCIUS 

C-Z0.05*Z73‘.S 

I 

1-1  (INDEX  FOR  INNER  LOOP  (POUER  RETURNED  TO  NOISE) 

Ind-1  (INDEX  FOR  SECOND  LOOP  (NOISE) 

Pou_back<0)-0  (INITIALIZE  VARIABLE  FOR  POUER  BACKSCATERRED 
I nt er val *( Speed_sound«Pul se«  1 . E- 3 ) / Z 
Ant_area*PI*) Ant_dia«/Z ) A2 
! 

(CALCULATE  RETURNED  SIGNAL  POUER  UNTIL  IT  IS  LESS  THAN  BACKGROUND 
(FOR  EACH  FREQUENCY. 

Range( 0)*0 
Noise- t.E- 18 
FOR  Ind-1  TO  S 

CALL  Attenuation) Aton_pres,Atten,Freq,Te«p, Wat _pres) 
K-Z*PI*Freq/Speed_sound  ! WAVENUMBER 
1-1 

REPEAT 

Range) I ) -Range) 1-1 >+2 
Sunpou_back -0 
SELECT  Profile 
CASE  1 

(THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
(AS  PRESENTED  IN  UAL TERS/K UNDE L  1S81  PAGE  398. 

!UAS  FOR  MIDDAY  AND  CLEAR  WEATHER  ABOVE  THE 
(TULAROSA  BASIN,  NEW  MEXICO. 

(AN  ADDITONAL  FACTOR  OF  1.8  WAS  INCLUDED  AS  PER  DIFFERENT 
(AVERAGING  TIME. 

CteZ-Z. 1Z ’Range) I )"(  -1 . IB) 

CASE  Z 

(THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
(AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 

(WAS  FOR  MIDOAY  AND  CLEAR  WEATHER  ABOVE  THE 
i TULAROSA  BASIN,  NEW  MEXICO. 
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! AN  AOOITONAL  FACTOR  OF  1.8  WAS  INCLUDED  AS  PER  DIFFERENT 
I  AVERAGING  TIME  AND  AN  ADDITIONAL  FACTOR  OF  TWO  WAS 
! INCLUDED  TO  APPROXIMATE  LOOKING  UP  A  THERMAL  PLOOM 
Cte2-Z*Z . 1 2* Range* I ) *  < - 1 . 1 G ) 

CASE  3 

t THIS  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 

I  PRESENTED  IN  UALTERS/KUNOEL  1981  PAGE  398  FOR" 

(NI6HT  CONDITIONS.  THE  DEPENDENCE  WITH  HEIGHT  " 

(IS  EXP* - . 001 *HEIGHT  ABOVE  65  METERS  ANO  HEIGHT" 

(TO  THE  -1.46  BELOU  65  METERS" 

IF  Range*  1X65  THEN 

CteZ-75 . 5* Range* I ) * ( - 1 . 46 ) 

ELSE 

CteZ-3. G6E-2*EXP< - .001 •Range* I ) ) 

END  IF 
CASE  4 

(CALCULATE  THE  TEMPERATURE  STRUCTURE  FACTOR  BASED  ON  ASSUMPTION 
!  THAT  IT’S  PROPORTIONAL  TO  HEIGHT  TO  THE  -4/3  POUER 

!  THAT  THE  SURFACE  VERTICAL  HEAT  FLUX  IS  .095  C.m./sec 

(THE  DEPENDENCE  UITH  HEIGHT  ABOVE  THE  INVERSION  I  PULLED  FROM 
(MY  LEFT  EAR 

(EQUATION  FROM  NEFF.  (975 
SELECT  Range* I > / I nver 
CASE  <.9 

C3-*< .024)*<T)‘< .667) > 

Cte2*C3*<  Kange* I ) )** -  1 . 33) 

CASE  .9  TO  1 
CteZ-Cte2 
C1*CteZ 
R1 -Range* I ) 

CASE  1  TO  (.3 

Cte2-10*<  < LGT* C3*Range* I >  >-LGT* C3*R1 ) )+LGT<  Cl  ) ) 

RZ “Range* I ) 

C2*CteZ 
CASE  ELSE 

CteZ“C3*<  Range* I )*( - 1 . 33>-RZ*( -  1 . 33) )+CZ 
ENO  SELECT 
END  SELECT 

(CALCULATE  VELOCITY  STRUCTURE  FACTOR.  FORMULA  FROM  GAYNOR  77 
Epsilon** .2866/Range* I ) )*( 1+.07*< Range* I > >"< .6) >*( 1 .5) 
CveZ“Z«Epsilon*< .667) 

(CALCULATE  ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  FACTOR. 

(FORMULA  FROM  TETARSKI 

CneZ* I )-( CteZ/( 2.98E+5) )  +  ( Cve2/<  C»C ) ) 

(CALCULATE  THE  FRACTION  OF  POWER  BACKSCATTERED  FROM  INTERVAL. 

I  EQUATION  FROM  NEFF.  1975 
Signa-*  .0039*<K,‘<  1/3>)*CteZ)/<T)“Z 
! 

(CALCULATE  THE  EXCESS  ATTENUATION.  EQUATION  FROM  CLIFFORD.  1980 
(THE  EXCESS  ATTENUATION  IS  Exc  att 


IF  I <5®  THEN 
Renal nder-0 

ELSE 

IF  I<1S0  THEN 

Renainder*Range( I )  MODULO  10 

ELSE 

Remainder* Range ( I )  MODULO  20 
END  IF 
END  IF 

IF  Renal nder *0  THEN 
Rho*0 
L-0 
H-0 

Rge-Range(I)  I  CONSTANT  IN  INTEGRAL 
R-0 

FOR  J-0  TO  2*1 

F-Cne2(INT< J/2+1 )) 

F*F*(  l-R/(Rge)',(  1 .67>+< R/< Rge) >*<1 .67) > 

IF  J>0  THEN 

IF  J<2*I  THEN 

IF  INT<  J/2 )*J/2  THEN 
L-L+F 
F-0 

ELSE 

H-H+F 
F-0 
END  IF 
END  IF 
END  IF 
Rho-Rho+F 
R-R+l 

NEXT  J 

Rho-Rho+4*L+Z*H 

Rho*(  ( ( Rho* . 33)*K*K* 1 . 46)*( - . B> ) 

N-<  Ant_diam/Rho)*Z 
IF  N<* 1  THEN 

Ze-1 /< 1 +N) 

ELSE 

Ze-1 .  5/( 14N) 

ISTEP  OF  1 . 5— >SEE  CLIFFORD  1980 
END  IF 

Exc_att-Ze*Ze 

ELSE 

E x c_al 1 - E x c_at 1 
END  IF 

! 

! CALCULATE  THE  POUER  BACKSCATTERED 

Pow_back< I >-(Pou_trane*Eff-Sunpog_back )*EXP( -Atten*Range( I ) ) 
Pou_back  < I >-Pou_backt I )*Interval*Exc_att*Signe 
Sunpoubaek-Sumpou  back+Pou  back< I ) 


I  CALCULATE  THE  POUER  BACK SCATTERED  TO  THE  ANTENNA 
Pow_ret( I )-Pou_back( I )*£XP< -At ten* Range ( I ) )*Ant_area*G*Eff/Range< I >*2 
PRINT  "RANGE-",  Ranged) 

PRINT  "POWER  RETURNED-".  Pou.retd) 

PRINT  “NOISE-", Noise 
I-I  +  1 

UNTIL  Pou_retd-1  )>0  AND  Pou_ret(  I- tX-Noise 
Ran< Ind)-Range< 1-2 ) 

Again- 1 
FOR  J-0  TO  I 
Pou_baek(  J)-0 
Pou_rei<  J)-0 
Range ( J )-0 
Cne2<  J)-0 
NEXT  J 

Noise-Noise*  <0 
NEXT  Ind 

) 

I 

OUTPUT  KBDj "K“ i 

PRINT  "INPUT  CONDITIONS" 

! 

PRINT  "  " 


PRINT  USING  "K"! "TEMPERATURE-  “ ,Tenp, “CELCIUS" 

PRINT  USING  "K"j "ATMOSPHERIC  PRESSURE-  " ,Aton_pres,"nb" 

PRINT  USING  “K" j "UATER  VAPOR  PRESSURE-  " .Wat_pres,"nb" 

PRINT  USING  "K"; "FREQUENCY-  " ,Freq,"  Hz" 

PRINT  USING  "K" ; "PULSE  LENGTH-  ".Pulse, "ns" 

PRINT  USING  "K" 5 "TRANSDUCER  EFFICIENCY-  ",Eff 

PRINT  USING  "K'VPOUER  TRANSMITTED-  "  ,Pou  trans,"  WATTS" 

PRINT  "TEMPERATURE  STRUCTURE  PROFILE  USED  ".Profile 
IF  Profile-4  THEN 

PRINT  "INVERSION  HEIGHT  ".Inver 
END  IF 

SELECT  Profile 
CASE  1 

PRINT  "  1  A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA" 

PRINT  "  AS  PRESENTED  IN  WALTERS/KUNDEL  1S81  PAGE  398  " 

PRINT  "  WHICH  GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 

CASE  2 

PRINT  "  2  A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA" 

PRINT  "  AS  PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398  " 

PRINT  "  WHICH  GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 

PRINT  "  BUT  WITH  A  FACTOR  OF  TWO  TO  APPROXIMATE" 

PRINT  "  LOOKING  UP  A  THERMAL  PLOOM" 

CASE  3 

PRINT  "  3  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA" 

PRINT  "  AS  PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398" 

PRINT  "  FOR  NIGHT  CONDITIONS.  THE  DEPENDENCE  WITH" 

PRINT  "  HEIGHT  IS  EXP< - .001 *HEIGHT  ABOVE  G5  METERS" 
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PRINT  “  AND  HEIGHT  TO  THE  -1.4E  BELOU  B5  METERS" 

CASE  4 

PRINT  “  4  A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 

PRINT  "  HEIGHT  TO  THE  -4/3  “ 

END  SELECT 

I 

I 

PRINT  "  " 

t 

I 

PRINT  “HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

OUTPUT  KBD; “K" ! 

! 

GRAPHICS  ON 
VIEUPORT  15.120.1®. 70 
! PLOT  RANGE  VERSUS  NOISE 
| 

Xnin— 18 

Xnax“- I  3 

Ywax«4 

Ynin«Z 

Dx-t 

Dy-1 

Yrange-2 

Xrange-5 

UINOOU  Xmin.Xnax , Ymin.Ynax 
AXES  Dx ,Dy .Xnin.Ynin, 1 , 1 
CLIP  OFF 


! LABEL  PLOT 
CSIZE  4,.G 
LDIR  0 
LORG  4 

MOVE  Xnin+.5*Xr ange ,Ynin+1*Yrange 
LABEL  "RANGE  vs  NOISE" 

! LABEL  HORIZONTAL  AXES 
LDIR  0 
LORG  5 

FOR  J-Xmin  TO  Xnax  STEP  Dx 
CSIZE  4. .6 

MOVE  J- .01 3*Xrange , Ynin- ,0E7«Yrange 
LABEL  "10" 

MOVE  J+ . 01 3*Xrange , Ymin- . 033* Yr ange 
CSIZE  2 
LABEL  J 
NEXT  J 

MOVE  Xnin+.5»Xrange,Ymin-. !Z«Yrange 
LORG  5 
CSIZE  4..G 
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LABEL  “NOISE  < WATTS)" 

! LABEL  VERTICAL  AXES 
LORG  8 

FOR  J-Ynin  TO  Yaax  STEP  Dy 
C SIZE  4, .6 

MOVE  Xnin-.01*Xrange, J 
LABEL  10 
CSIZE  Z 

MOVE  Xmin-.0OZ5*Xrange, J+.03*Yrange 
LABEL  J 
NEXT  J 
LOIR  PI/Z 
LORG  G 

MOVE  Xmin-. 1 *Xrange , Ymin+ . 5* Yrange 
CSIZE  4,.B 

LABEL  "RANGE  (METERS)" 

CLIP  ON 
Ind-1 

Noise- 1 . E- 18 
FOR  Ind-1  TO  5 

PLOT  LGT( Noise) ,L6T<  Ran< Ind) > 

Noise-Noise* 10 
NEXT  Ind 
( 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT  KBD;"K"i 
! 

( 

LINPUT  "UOULD  YOU  LIKE  TO  MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?",X* 
Aoain-FNYes( X*) 

Mess_up- 1 
WHILE  Mess_up-1 
SELECT  Again 
CASE  1 

LINPUT  "WOULD  YOU  LIKE  TO  CHANGE  A  VARIABLE  (ENTER  Y  OR  N>?" ,X* 
Neu_va-FNYes(X*> 

Ano_change- 1 
WHILE  Ano_change-1 
SELECT  Neu_va 
CASE  1 

PRINT  "  VARIABLE  CURRENT  VALUE" 

PRINT  USING  "K"s"1  TEMPERATURE  " rTenp , "CELSIUS" 

PRINT  USING  " K ”  ; " 2  ATOMOSPHERIC  PRESSURE  "  . Aton_pres  ,  "mb 
PRINT  USING  "K" t "3  WATER  VAPOR  PRESSURE  " , Wat.pres , "mb" 
PRINT  "  " 

PRINT  USING  "K" i "4  ANTENNA  DIAMETER  " ,Ant_dian,"  n." 

PRINT  USING  "K"s"5  PULSE  LENGTH  ".Pulse,"  as” 

PRINT  USING  "K" ; "G  POWER  TRANSMITTED  " ,Pou_trans ,"  WATTS 


PRINT  USING  “K" i "7  ATMOSPHERIC  PROFILE  ".Profile 
PRINT  -  - 

PRINT  "ENTER  THE  NUMBER  OF  THE  VARIABLE  YOU  UISH  TO" 

PRINT  "CHANGE" 

INPUT  Var 
SELECT  Var 

CASE  t 

IfffUT  "TEMPERATURE-", Temp 
CASE  2 

INPUT  "ATMOSPHERIC  PRESSURE  IN  mb-" , Atom_pres 
CASE  3 

INPUT  "WATER  VAPOR  PRESSURE  IN  mb-" , Wat_pres 
CASE  4 

INPUT  "ANTENNA  DIAMETER  IN  m. -“ . Ant.dian 
CASE  S 

INPUT  "PULSE  LENGTH  IN  me-" .Pulse 
CASE  G 

INPUT  "POWER  TRANSMITTED  IN  WATTS-" .Pou  trans 
CASE  7 

PRINT  “  t  — >  Ct*Z  PROFILE  OF  Z*<-1.16>" 

PRINT  "  FROM  WALTERS/KUNDEL  1981" 

PRINT  "  2  — >  SAME  AS  ONE  BUT  WITH  FACTOR" 

PRINT  "  OF  TWO  TO  APPROXIMATE  LOOKING" 

PRINT  "  UP  A  THERMAL  PLOOM" 

PRINT  "  3  — >  Ct*Z  PROFILE  OF  EXP< - .001 *Z > " 

PRINT  "  FROM  WALTERS/KUNDEL  1981." 

PRINT  "  4  — >  CT*Z  PROFILE  OF  Z*<-4/3)" 

INPUT  "ENTER  NUMBER  OF  DESIRED  PROFILE" .Prof lie 
IF  Profile-4  THEN 

INPUT  "HEIGHT  OF  INVERSION  IN  METERS-" , Inver 
END  IF 
CASE  ELSE 

PRINT  Var. "IS  NOT  ONE  OF  THE  OPTIONS" 

END  SELECT 

LINPUT  "MADE  ANOTHER  CHANGEt ENTER  Y  OR  N)7" .X* 
Neu_va-FNYes<  X$) 

Mese_up-Z 
CASE  2 

Mess_up-Z 
Ano_ change-2 
CASE  ELSE 

PRINT  "YOUR  CHOICES  WHERE  Y  OR  N  !!!  !  !  “ 

Ano_change- 1 

LINPUT  "CHANGE  A  VARIABLE  .<  ENTER  Y  OR  N)?",X$ 
Neu_va-FNYes(X*> 

END  SELECT 
END  UHILE 
CASE  2 

Mess_up-Z 
CASE  ELSE 


PRINT  "YOUR  CHOICES  WHERE  Y  OR  N  !!!!!!" 
Mess_up-1 

LINPUT  “MAKE  ANOTHER  RUN  < ENTER  Y  OR  N)?".X* 
Again-FNYea<  X#  > 

ENO  SELECT 
END  UHILE 
ENO  UHILE 

PRINT  "THAT’S  ALL,  FOLKS" 
end 


ICALCULTE  THE  ATTENUATION 

sub  Attenuat ion< At om_pres , Atten .Freq , Temp , Uat_pr es ) 

I  THIS  SUBPROGRAM  CALCULATES  THE  ATTENUATION  OF  SOUND 
! IN  AIR  BASED  UPON  EQUATIONS  IN  NEFF  1975 

I 

l  INPUT  ATMOSPHERIC  PRESSURE  IN  MILLIBARS 
|  FREQUENCY  OF  SOUND  WAVE  IN  HERTZ 

|  TEMPERATURE  IN  DEGREES  CELCIUS 

I  WATER- VAPOR  PRESSURE  IN  MILLIBARS 

! 

!  OUTPUT  ATTENUATION  IN  1 /METERS 


VARIABLES 
Atom_pres 
At  ten 

Att_max 


F 

Fmax 

! 

!  Freq 
!  H 

!  Pstar 
! 

I  Temp 
I  Tstar 
I 

!  Uat_pres 


INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 

ATTENUATION  OF  ACOUSTIC  WAVE.  CALCULATED  IN 
SUBPROGRAM  ATTENUATION. 

VARIABLE  IN  SUBPROGRAM  ATTENUATION.  IT  IS  THE 
ATTENUATION  AT  THE  FREQUENCY  OF  THE  MAXIMUM 
ATTENUATION  FOR  THE  INPUT  CONDITIONS. 

VARIABLE  USEO  IN  SUBPROGRAM  ATTENUATION.  IS  THE 
RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 

FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 
ATTENUATION. 

INPUT  FREQUENCY  OF  ECHOSOUNDER. 

VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 
CALCULATION. 

INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 
IN  SUBPROGRAM  ATTENUATION. 

ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.  INPUT  BY 
OPERATOR. 


H“ 100*Wat_pres/Atom_pres 
Tatar»< 1 .8«Temp449Z ) / S 1 9 
Pstar-Atom_pres/ 1014 

Fmax-(  10+si0C»H+4440®*H*H>  •Pstar-/ Tatar*.  8 
Att_nax-.0078*Fnax«Tstar*<-Z.S>»EXP<7.77»<  1- 1/Tstar >  > 
F*Freq/Fmax 

Atten-<Att_max/304.8>»<< . I B*F )* Z+< 2»F*F/<  » +F*F ) >*Z )* .5 
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-v, 


•-  ..v.v'. 


Atten-!  Atten+ 1 .74E- 1®*Freq*Freq>/4. 35 
SUBEND 
1 

def  FNYes(  X$ ) 

! THIS  FUNCTION  INTERPRETS  THE  OPERATORS  RESPONSES  TO  YES  NO  QUESTIONS 
! 

!  INPUT  XS 

I 

!  OUTPUT  FNYes 

(VARIABLES 

!  TenpS  VARIABLE  STRING  USED  IN  FUNCTION  YES. 

!  XS  STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

!  RESPONSE  TO  YES  OR  NO  QUESTION. 

DIM  TempSin 
TempSM  ,11 -TRIMS!  XS> 

SELECT  TempS 
CASE  "Y“ ,  "y" 

RETURN  1 
CASE  "N"  ,'*n“ 

RETURN  2 
CASE  "  M 

RETURN  1 
CASE  ELSE 

RETURN  -2 
END  SELECT 
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